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Kin Quadraturverfahren zur Berechnung der laminaren und turbulenten 
Reibungsschicht bei ebener und rotationssymmetrischer Strémung *. 


Von E, Truckenbrodt, 


1, Einleitung. Die Berechnung der laminaren und turbulenten Reibungsschicht bei Druck- 
abfall und Druckanstieg ist das entscheidende Problem bei der Bestimmung des Strémungs- 
verlustes in Kanalen und Rohren und des Strémungswiderstandes von Kérpern. Dabei mige 
im folgenden sowohl das ebene wie auch das rotationssymmetrische Problem behandelt werden. 
Bei unseren Betrachtungen wollen wir uns auf inkompressible Strémungen beschranken. 

Wahrend die theoretische Behandlung der laminaren Reibungsschicht grundsatzlich geklart 
ist, ist man bei der Behandlung der turbulenten Reibungsschicht immer noch auf halbempirische 
Zusammenhange angewiesen. Da im laminaren Fall die exakten Methoden recht langwierig sind, 
sind hierfiir, genau so wie im turbulenten Fall, verschiedene Naherungsverfahren entwickelt 
worden. Allen diesen Verfahren liegt die erstmalig von Th. von Kdérmdén! angegebene Impuls- 
gleichung zugrunde. 

Das Naherungsverfahren fiir den ebenen laminaren Fall wurde von K. Pohlhausen 2 aus- 
gearbeitet. Es wurde spater von H. Holstein und T. Bohlen® dadurch ganz wesentlich verbessert, 
daB nicht die Grenzschichtdicke wie bei Pohlhausen, sondern die Impulsverlustdicke als gesuchte 
GréBe eingefiithrt wurde. In neuerer Zeit hat nun K. Wieghardt* zusatzlich zu der Impuls- 
gleichung noch eine Energiegleichung abgeleitet, die er benutzt, um ein zweiparametriges Rechen- 
verfahren zu entwickeln. A. Walz® hat dieses Verfahren dadurch vereinfacht, daB er wieder zur 
Einparametrigkeit wie bei den Verfahren von Pohlhausen und Holstein- Bohlen zuriickkehrt. Das 
dem Pohlhausen-Verfahren analoge rotationssymmetrische Verfahren wurde von S. Tomotika ® 
angegeben. Seine Vereinfachung im Sinne des Holstein-Bohlenschen Ansatzes wurde von F. W. 
Scholkemeyer 7 vorgenommen. 

Zur Berechnung der ebenen turbulenten Reibungsschicht wurde erstmalig ein Verfahren 
von E. Gruschwitz* angegeben. Auber dem Impulssatz benutzt Gruschwitz noch eine aus gewissen 
Energiebetrachtungen gewonnene halbempirische Gleichung zur Bestimmung eines das Ge- 
schwindigkeitsprofil in der Reibungsschicht kennzeichnenden Formparameters, welcher die Ab- 
léseempfindlichkeit der Reibungsschicht charakterisiert. Die in dem Gruschwitzschen Verfahren 
auftretenden empirischen GréBen wurden von A. Kehl® erneut untersucht. Kin dem Gruschwitz- 
schen 4hnliches Verfahren wurde auch von E. Buri! entwickelt. In Amerika hat sich 
ein Verfahren von A. E. von Doenhoff und N. Tetervin™ als brauchbar erwiesen. Neben dem 
Impulssatz wird auch bei diesem Verfahren eine empirische Gleichung zur Bestimmung eines 
das Geschwindigkeitsprofil kennzeichnenden Formparameters benutzt. Ausgehend von dieser 


* Bericht des Instituts fiir Strémungsmechanik der T.H. Braunschweig. Die Untersuchung wurde 
mit Unterstiitzung der Deutschen Forschungsgemeinschaft durchgefiihrt. 

1 Th. v. Karman, Z. angew. Math. Mech. 1 (1921), S. 233. 

2 K. Pohlhausen, Z. angew. Math. Mech. 1 (1921), S. 252. 

3 H. Holstein u. T. Bohlen, Lilienthal-Bericht S 10 (1940), 5. 5. 

4 K. Wieghardt, Ing.-Arch. 16 (1948), S. 231, 

5 A, Walz, Ing.-Arch. 16 (1948), S. 243. 

6 §, Tomotika, Laminar boundary layer on the surface of asphere in a uniformstream, ARC-Rep. No, 1678 

1935). 

g e W. Scholkemeyer, Die laminare Reibungsschicht an rotationssymmetrischen Kérpern, Diss. Braun- 
schweig (1943). 

8 EH. Gruschwitz, Ing.-Arch. 2 (1931), S. 321. 

9 4. Kehl, Ing.-Arch. 12 (1943), S. 293. 

10 E, Buri, Eine Berechnungsgrundlage fiir die turbulente Grenzschicht bei beschleunigter und ver- 
zogerter Strémung, Diss. Ziirich (1931). 

11 4. E. v. Doenhoff u. N. Tetervin, Determination of general relations for the behaviour of turbulent 
boundary layers, NACA Rep. No, 772 (1943). 
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Arbeit hat H. C. Garner! in England eine Methode entwickelt, die in ihrer numerischen Aus- 
wertung derjenigen von v. Doenhoff- Tetervin itberlegen ist. Die Ubertragung des Impulssatzes auf 
den rotationssymmetrischen Fall wurde von C. B. Millikan? vorgenommen. 

Durch neuere Forschungen von H. Ludwieg und W. Tillmann 3, sowie J. Rotta * tiber die Ge- 
setzmaBigkeiten der turbulenten Strémungen, insbesondere der Wandschubspannung und des 
Energieverlustes in der Reibungsschicht, ist es moglich, die bestehenden halbempirischen Ver- 
fahren besser zu begriinden und zu verbessern. 

Es ist das Ziel der vorliegenden Arbeit, ein Rechenverfahren zu entwickeln, das gleicher- 
mafen fiir die vier Fille laminarer und turbulenter, sowie ebener und rotationssymmetrischer 
Strémung gilt. Dabei sollen die neuesten Ergebnisse Beriicksichtigung finden. 


2. Der Impulssatz und Energiesatz. Die Ausgangsgleichungen zur Berechnung der Reibungs- 
schichten bilden in bekannter Weise die Prandtlschen Grenzschichtgleichungen und die Bernoulli- 
sche Gleichung. Wie wir schon in der Einleitung 
sagten, wollen wir inkompressible Strémungen be- 
trachten. Da sich aber die nachfolgenden Ableitun- 
gen ohne besondere Schwierigkeiten fiir verander- 
liche Dichte 9 angeben lassen, wollen wir den Uber- 
gang zur inkompressiblen Strémung erst im End- 
ergebnis vornehmen. 

Wir nehmen an, dafB im rotationssymmetrischen 
Fall der Kérper- bzw. Rohrradius R gegeniiber der 
Reibungsschichtdicke 6 groB ist. Dann lauten die 
besagten Gleichungen 


Qu, 0u\ ep ox 
e (« du! 5) Meio, hoy. Ce 
rotationssymmerrisch ae dy ? 
Op fay dU 
Agee ae iUe U Woes (3) 


Hierin bedeuten u und v die Geschwindigkeits- 
komponenten innerhalb der Reibungsschicht 


ORS ye) 


in x- und y-Richtung und U die Geschwindigkeit 

auBerhalb der Reibungsschicht (y > 6). Abb. 1. Hierbei ist p der iiber die Reibungsschicht- 

dicke konstante Druck, t die Schubspannung und @ die Dichte des strémenden Mediums. Fiir 

den ebenen Fall ist in (2) sowie auch in allen folgenden Formeln der Halbmesser R zu streichen. 
Die zugehérigen Randbedingungen lauten: 


Abb. 1. Ubersichtsskizze. 


y= OO, Uap 1 | 
a ae 2 By (4) 
Wes Ot We Ul, tO CABS | 
Ist vp == 0, so liegt der Fall der Absaugung (vy) < 0) bzw. des Ausblasens (vp > 0) vor. 
Die drei Gleichungen (1) bis (3) fassen wir jetzt noch zusammen, indem wir zu der mit u” 
multiplizierten Gleichung (1) die mit u™+!/(r + 1)R multiplizierte Gleichung (2) addieren, wobei 
r=0,1,2... beliebig gewahlt werden kann. Es wird dann 


ik @ 5 T 
RROVPR +z owt) = 6+ Yw (eu e +5). (5) 


Die Kontinuitatsgleichung (2) integrieren wir itber den Wandabstand y und erhalten 


ot 
eck ys 
ov — Rx | eURdy + eon. (6) 


0 


' H. C. Garner, The development of turbulent boundary layers, ARC-Rep. No. 2133 (1944). 
> C. B. Millikan, Trans. Am. Soc. Mech. Eng., Appl. Mech. 54 (1932), No. 2, S. 29. 

3 H. Ludwieg u. W. Tillmann, Ing.-Arch. 17 (1949), S. 288. 

4 J. Rotta, Ing.-Arch. 19 (1951), 5. 31 und 20 (1952), S. 195. 
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Nehmen wir jetzt an, dal uns die Geschwindigkeitsverteilung u(y) und die Schubspannungs- 

verteilung t(y) in y-Richtung bekannt sind, dann kénnen wir (5) tiber y von y= 0 bis y= 0 
integrieren. Wir erhalten dann, wenn wir noch (6) einsetzen, 

1 d 

Qa Ui : R dx 

Hierin bedeuten die neu eingefiihrten Abkiinzungen 


ry 2 | gr dU QorVor = 
RG) ot ae a he (7) 


0) 


non we u\r+1 
i rr aad ca 
0 


sen fe] (eo) 


Gleichung (7) gilt ganz allgemein fiir laminare und turbulente, inkompressible und kompressible, 
ebene und rotationssymmetrische Strémung mit und ohne Absaugung bzw. Ausblasen. Wie be- 
reits gesagt, ist im ebenen Fall der Radius R zu streichen. 

Setzt man r= 0, dann ergibt sich die bekannte v. Kérmdnsche Impulsgleichung. K. Wieg- 
hardt} und J. Rotta® haben gezeigt, da man (7) auch fiir r— 1 eine physikalische Deutung 
geben kann, namlich die der Energiegleichung. Wir wollen jetzt annehmen, daf das strémende 
Mittel inkompressibel sei, 9 = konst, auSerdem sei weder Absaugung noch Ausblasen vorhanden, 
Vy = 9. Dann wird 


r= 0: I mpulssat z (v. Kérmén, Pohlhausen). 
igs Bea d*dU — t 
Gre ee ae Pe Oy Ie (°) 
mit 
yore 
O fe | a i = = dy als Impulsverlustdicke, (10) 
0 
e 
Ot == op = 1 r| dy als Verdringungsdicke, (11) 
0 
TL =e als Wandschubspannung. 
r=1: Energiesatz (Wieghardt, Rotta). 
1d spay—9¢tt 12 
B Rae OR) 2 gs me 
mit 
eC é 
haan Pe a L —_ tal dy als Energieverlustdicke, (13) 
ee 5 
d+t 4 pa a als Schubs sarbeit *. 14 
i ae a i 0 UP ay ir) als Schubspannung (14) 


Die Schubspannungsarbeit ist im laminaren Fall gleich der in Warme umgewandelten Energie 
(Dissipation d). Bei turbulenten Str6mungen wird nicht alle Energie in Warme umgewandelt ; 


1 Siehe Fu note 4 auf Seite 211. 
2 Siehe Fufnote 4 auf Seite 212. ‘ 
3 Diese Gleichung wird durch partielle Integration aus (8) gewonnen. 
Ge 
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es verbleibt noch ein Teil als Turbulenzenergie (t), der allerdings meistens vernachlassigbar klein 
ist. Wichtig ist die Tatsache, daB im Fall r = 1 das zweite Glied auf der linken Seite von (7) ver- 
schwindet, da ja bei Annahme konstanter Dichte g, = 0 ist. 

Fithren wir jetzt die Reibungsschichtdickenverhdltnisse 


__ oF Te 15 
H=—= und H=— (15) 

ein, dann kénnen wir fiir (9) und (12) auch schreiben 
] d - | o dU = To é 16 
TRga he ahi a meat ne 
RC RH) —2 Coe 

Wir subtrahieren (16) von (17) und finden 1 

dH pe a et aon 18 
p= (H— 1) | 2 i (18) 


Wie im folgenden Abschnitt gezeigt wird, lassen sich die Schubspannung T,/9 U? und die Schub- 
spannungsarbeit (d+ 1)/9 U? als Funktionen der dimensionslosen GriéBen Ud/y (mit der Impuls- 
dicke gebildete Reynoldssche Zahl) und H ausdriicken. Weiterhin gilt in guter Naherung Ein- 
parametrigkeit fiir die Geschwindigkeitsverteilungen in der Reibungsschicht; d.h. es besteht 
ein fester Zusammenhang zwischen den Reibungsschichtdickenverhaltnissen H und H. Unter 
- Beachtung dieser Tatsachen stellen 
ean die Gleichungen (17) und (18) zwei 
os Gleichungen fiir die beiden Un- 
bekannten @ und H dar. 

Im Gegensatz zu den bisherigen 
Verfahren wollen wir die Impuls- 
verlustdicke nicht aus dem Impuls- 
| Gleicharuck satz (16) sondern aus dem Energie- 
| satz (17) berechnen. Dieses hat fiir 
16 Ablisung | die Ausfithrung der im Abschnitt 4 
0c | angegebenen Integrationen entschei- 
one | i | dende Vorteile. Den Impulssatz 
nehmen wir in Verbindung mit dem 
| Energiesatz zu Hilfe, um den das 
a7 Geschwindigkeitsprofil kennzeich- 


OY 


3 


H | nenden Parameter H zu berechnen, 
y Gleichung (18). 
20 30 40 3. Die Schubspannung, die Schub- 
Abb. 2. Laminare ReibungsschichtgréBen (Hartree- Profile), spannungsarbeit und das Reibungs- 
schichtdickenverhaltnis, Um mit den 
Gleichungen (17) und (18) arbeiten zu kénnen, mu die Abhiangigkeit der Schubspannung und 
der Schubspannungsarbeit von den GréBen U#/y und H bekannt sein, auBerdem mu man den 
Zusammenhang H(H) kennen. 


a) Laminare Strémung. Unter der Annahme der Einparametrigkeit fiir die Geschwindigkeits- 
verteilung kénnen wir schreiben 
Neigh 
eh ( 5H). 


woraus dann folgt 


a ie 
To __ 0 (A) ; _ (ix 
eu Usb cee en oa ere (19) 
=e C= 
v 9/0 


1 Fir den laminaren Fall wird eine entsprechende Formel bereits von A. Walz (FuBnote 5 auf Seite 211) 
angegeben. Auch im turbulenten Fall kann man zeigen, dafs die von E. Gruschwitz, A. E, v. Doenhoff und 
N. Tetervin sowie H. C. Garner angegebenen empirischen Gleichungen zur Berechnung des Formpara- 
meters des Geschwindigkeitsprofiles (bei Gruschwitz 1 = 1 — (uw9/U)?, bei den anderen H) im Aufbau der 
Gleichung (18) sehr 4hnlich sind. 
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6/0 


d Bi(H) 


ue = mit b= - di > (20) 
und 
5/0 3/0 
= Pare ES | a —(e) [45 (21) 
5° 0 


Wir wahlen fiir die Auswertung dieser Formeln die sogenannten Hartree-Profile; das sind die 
exakten Profile fiir den Geschwindigkeitsverlauf U ~ x™. In Abb.2 haben wir die GréBen a, Bs 
H itber H dargestellt 1. Interessant ist, daB f nahezu von H unabhangig ist. Fiir den Fall ohne 
Druckgradient gilt «= 0,220, 8 = 0,173 und H= 2,60. 


b) Turbulente Strémung,. Mit der Bestimmung der Wandschubspannung bei turbu- 
lenten Reibungsschichten mit Druckgradienten haben sich neuerdings H. Ludwieg und W. Till- 
mann sowie J. Rotta? eingehend befaBt. 

Fiir den Bereich der Reynoldsschen Zahlen 1 - 10° << U#/y < 4+ 104 geben Ludwieg- Tillmann 
die folgende Interpolationsformel an: 


To 0,123 —0,678 H 
ou? ==: e apr 5 (22) 
“. 


Rotta findet fiir die Wandschubspannung die Beziehung 


eed ieee B. (23) 


Hierin bedeutet u* = Vmle die Schubspannungsgeschwindigkeit. Es gilt somit 


Bh et . (24) 


@ U? U 


In (23) ist 1/x—= 2,5 und B= B(I,) eine Funktion der GroBe I, = = = : a Die Funktion B 


wurde von Rotta ausgewertet und graphisch dargestellt. 

Wir haben die Schubspannungswerte fiir verschiedene Werte von H berechnet und in Abb. 3 
itiber der Reynoldsschen Zahl U#/y aufgetragen. Die Ubereinstimmung zwischen den Werten 
nach Ludwieg-Tillmann und Rotta ist recht befriedigend. 

Mit der Berechnung der Schubspannungsarbeit (Dissipation und Turbulenzenergie) 
hat sich ebenfalls J. Rotta befaBt. Fiir die Dissipation gilt 


co fle 88) +4 


Die hierin auftretende Funktion G = G(I,) ist von Rotta ausgewertet worden. Unter Beachtung 
von (23) kénnen wir fiir (25) auch schreiben 


i=) p HFe—m) es 


Da eo Ue H), sowie , = oS ! a ist, laBt sich die Dissipation ebenfalls als Funktion 
Vv 


der Reynoldsschen Zahl U#/y und des Formparameters H darstellen, was in Abb. 4 geschehen 
ist. Es zeigt sich, da® die Unterschiede bei verschiedenen Werten H nur gering sind, wahrend 
sie bei der Wandschubspannung sehr erheblich sind. Aus Abb. 4 entnehmen wir, daB sich die 


1 Wir haben die Zahlenwerte von A. Walz (FufBnote 5 auf Seite 211) itbernommen., Und zwar gelten 
folgende Identitaten: Ose fp, =D* und H =H, 
2 Siehe FuBnote 3 und 4 auf Seite 212. 
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Dissipation durch den Ansatz 


sai . 


annahern lift. Ein brauchbarer Wert fiir n ist n = 1/6. Die Abhangigkeit des Wertes /, von H 
ist nur gering, so da®& wir dafiir den konstanten Wert 6 = 0,56- 10° annehmen kénnen. 
Zur Berechnung der Turbulenzenergie gibt J. Rotta eine Naherungsformel an. An Hand dieser 
Formel kann man zeigen, daB die Turbulenzenergie gegeniiber der Dissipation vernachlassigbar 
klein ist, eine Tatsache, auf die bereits Rotta hingewiesen hat. Wir setzen also fiir unsere weiteren 


ae 
U 
4,0 2 
SS 
ES 
60 Sn 
IS SS 
Se — 
~S H=42 
SS 
40 S 
SN a 
08 | 
SS | 16 
06 = <= 
= ~ 
SS > 48 
SS =~ 
Ss + + 
O4 One. = ~ 20 O4 
SS 
= SS 
ieee 
as 
x: a4 Rotta 
a2 7 = gz laa A ee | 
aes Sa 116 
SS UD 
Rotta ov (49) 
———Ludwieg/Tillmann ; 
Ue YY 
iy Oem St ar I a eh 2 
es ee ee yee SP he Suan os 2 Sli, 
Abb. 3. Turbulente Wandschubspannung (nach H. Ludwieg, Abb. 4, Turbulente Dissipation (nach J. Rotta). 
W. Tillmann und J. Rotta). 
t 
Rechnungen me OF (27) 


Das Reibungsschichtdickenverhaltnis H (Energieverlustdicke/Impulsverlustdicke) 
kénnen wir ebenfalls aus den Ergebnissen von J. Rotta bestimmen. Es gilt 


(28) 


Hierin ist I, ein Formparameter, der mit J, in einem festen Zusammenhang steht, der von Rotta 


angegeben worden ist. Wir haben nach (28) die Werte Hin Abhangigkeit von der GréBe H, sowie 
von der Reynoldsschen Zah] U®/y berechnet. Dabei zeigt sich, da der EinfluB der Reynoldsschen 
Zahl verschwindend ist. Das Ergebnis unserer Rechnung haben wir in Abb. 5 dargestellt. Auch 


K. Wieghardt! hat sich mit dem Zusammenhang H(H) beschaftigt. Er findet 


i AH 
fe eer 


mit A=1,269 und B=0,379. (29) 


‘ K. Wieghardt, Zur turbulenten Reibungsschicht bei Druckanstieg. Untersuchungen und Mitteilungen 
der deutschen Luftfahrtforschung, UM 6617 (1944). 
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Die sich hieraus ergebende Kurve haben wir ebenfalls in Abb. 5 eingetragen. Bis auf die gréBeren 


Werte von H ist die Ubereinstimmung mit den Kurven von Rotta gut. 


4, Das Niherungsverfahren, Als Ergebnis des vorigen Abschnittes wollen wir zusammen- 
fassen: Es gelten die Darstellungen 


Ci Ud d+t d U & = 
muiaal| - HY). ou ots a(t) und H=—=h(H), 


welche in Abb. 2 bis 5 angegeben sind. Im laminaren Fall sind die Werte 7/0 U? und d/o U3 um- 
gekehrt proportional der Reynoldsschen Zahl. 


a) Die Berechnung der Im- 7 
pulsverlustdicke. Wie bereits 
oben erwahnt, wollen wir die 
Impulsverlustdicke aus dem 
Energiesatz (17) ermitteln. 

Wir schreiben fiir die Schub- iN | | 
spannungsarbeit nach (20), (26) 
und (27) 


d+t  B(H) 
ie: ie a eas 
ge 


wobei fiir den laminaren Fall 
n= 1 und fiir den turbulenten 
Fall n = 1/6 gilt. 
Diesen Ausdruck setzen wir 
jetzt in (17) ein und erhalten 
mit 47 
1= U3+2 Rit+n Hi+n (2) 9, 
v 


(31) 


o% — 2 (1+ n) B(H) 
H” U3+2 Rit», (32) 


Nehmen wir jetzt an, dafB H 
von x und durch den eindeutigen 
Zusammenhang H(H) auch H 
bekannt seien, dann 1a4Bt sich 
(32) tiber x integrieren, und es 
ergibt sich, wenn wir noch die 


GroBe 


Rota 
| ———Wieghards 


U d\n 
O= icon ean A 
i TG 2 14 16 18 20 22 28 
einfiihren, unter Beachtung Abb. 5. Turbulentes Reibungsschichtdickenverhaltnis 
von (31) (nach J. Rotta und K. Wieghardt). 


“ta (Heo? Ro dig 
O (x) ba ; (34) 


i “H(x)i+* U(x)? + 2n R(x)! n 
Als neue Abkiirzung haben wir eingefiihrt 


E(H)=2(1-+ n) 6(H) H”. (35) 


Diese Funktion haben wir fiir die beiden Falle laminarer und turbulenter Strémung in Abb. 6 


aufgetragen. os 
Naherungsweise wollen wir jetzt fiir E(H) und H konstante Mittelwerte annehmen. Diese 


Annahme trifft besonders gut fiir den turbulenten Fall zu. Wir erhalten dann aus (34) 
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Cet Uet2” RI TH dng! 


o(2) = (2) 9= aos (36) 
Hierin bedeutet 
=an= 2(1 + n) 2 (37) 
einen fiir den Bereich x,< x’ < x brauchbaren Mittelwert. 
Die Integrationskonstante bestimmt sich zu 
C= UpeR (ESY 9, (38) 


mit n= 1/6, wenn vom Punkt x, ab turbulente Reibungsschicht vorliegt. Wir wollen als lami- 


naren bzw. turbulenten Mittelwert fiir A denjenigen Wert wablen, der sich ergibt, wenn wir 
U = U,, konstant annehmen und die 


2 

a a z a2 ebene Platte bei voll-laminarer bzw. 
voll-turbulenter Strémung betrachten. 
| Gleichdruch — Aus (36) wird dann mit x,—0 und 

| ew ae Q= 

we: (Sh Op n 

. turbulent | 0,= ( ; 5 ale (39) 
: | Zwischen der Impulsverlustdicke @, (J) 
40 ob und dem Widerstandsbeiwert cy einer 
. | laminar einseitig benetzten Platte der Lange I, 
: | die mit der Geschwindigkeit U,, an- 
3 | gestrémt wird, besteht der Zusammen- 

t hang 

Dp (I) oi 

| H Shy ae ts cou) 
"19 2,0 30 40 Setzen wir in (39) x= I, dann folgt 


durch Vergleich von (39) und (40) 


Abb. 6. Die Funktion E(H) bei laminarer und turbulenter Strémung. 
we I\" cf 1+n 
Je nach dem vorliegenden Strémungszustand ist c, fiir laminare a im i “(37 a) 
bzw. turbul Stré Abb.7). Di ibwei 
zw. turbulente Strémung zu nehmen ( 7). Die Schreibweise Tahelle le Die Cragin anda 


(37 a) bietet noch den Vorteil, daB man auch die Oberflachen- hat loinares ae 
rauhigkeit leicht beriicksichtigen kann, indem man lediglich Strémung. 


die entsprechenden c,-Werte einer rauhen Platte einsetzt. Fiir lao eR 
glatte Flachen ergeben sich die Werte von Tabelle 1 fiir die | (Blasius) (Falkner) ? 
Konstante A. n 
Wir lésen jetzt (36) nach der Impulsdicke ? auf und er- i ; 6 
halten, wenn wir noch dimensionslose Gréfen einfiihren unter A 0,441 0,760 - 10-2 
Beachtung von (37a) den folgenden Ausdruck: 
in x/1 =i il 
a” [er aynap 
D ue I i 
O(n) afl | 
—_— TR é (41) 
(rz) an 


Die Integrationskonstante bestimmt sich zu 


x, |t = 
cx — [( Ua Ri apttn_ [eps Ae Re ae 
3 Uni) Tata me a, Gs 7) on (42) 
0 
Zusammenfassend wiederholen wir noch einmal: Ist die Strémung vom Anfangspunkt £, =Otan 
voll-laminar bzw. voll-turbulent, dann ist Cf = 0. Es sind jeweils die laminaren bzw. turbulenten 
0,0306 


‘ca (Z= ar : 
Vv 


' V. M. Falkner, Aircraft Engineering, 15, (1943), S. 65. Es gilt: ¢ 
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Widerstandsbeiwerte einzusetzen. Geht der turbulenten Reibungsschicht eine laminare Vor- 
strecke voraus, dann ist der Wert fiir # in (42) der laminaren Rechnung zu entnehmen (zweite 
Formel in (42)). Fir den laminaren Fall gilt n= 1, fiir den turbulenten Fall n— 1/6. Der 
ebene Fall ergibt sich, wenn man in (41) und (42) den Radius R streicht. Wie man aus (41) er- 
sieht, spielt die GréBe n im turbulenten Fall keine besondere Rolle. so daf man fiir iber- 
schlagige Rechnungen auch n = 0 annehmen darf. 


») Die Berechnung des Formparameters, Den das Geschwindigkeitsprofil in der Reibungs- 
schicht kennzeichnenden Formparameter H wollen wir aus Gleichung (18) bestimmen, die wir 
wie folgt schreiben: 
dH 


O (x) — = F(H) I(x) + G(H). (43) 


ce: 10° 


O74 | — 


Yoo i 4 


Prandtl Schlichting, op = We, 
( g =!) 


_turbulent 


s —— 06,0306 
Ubergangsbereich re ee ka at AG 
ee ear 


—= 


+— 


(oe 
v 
a 70° 70’ 1° 709 
Abb. 7. Widerstandsgesetz der glatten langsangestromten ebenen Platte. 
Hier ist O nach (33) gegeben. Fiir die anderen Abkiirzungen gilt 
0 dU 
P(x) = Tas’ (44) 
F(H) =(H—1)H 
me d+ LUONe > 
oe taderee ve 
con (tit 5) 
Wir formen (43) noch um, indem wir die Substitution 
= dH 
“m= {| = L (Hf) (46) 
und die Abkirzung 
Fy 28 CH) 
K (Hf) = — $09 = K(L) (47) 
einfiihren. Dann erhalten wir eine Differentialgleichung fiir den neuen Formparameter L 
@(x) 5 = I(x) — K(L). 7 (48) 


dx 


Durch graphische Integration der Funktion 1/F(H) ttber H haben wir die Grife L bestimmt. 
Den Wert L — 0 haben wir in den Bereich verschwindenden Druckgradienten (Gleichdruck, 
ebene Platte) gelegt. In den Abb. 8 und 9 haben wir die Abhangigkeit H(L), sowie fiir den lami- 
naren Fall auch a(L) und in den Abb. 10 und 11 die Abhangigkeit K(L) fiir die Falle laminarer 
und turbulenter Strémung dargestellt. Wahrend fiir die Funktion K(L) im laminaren Fall kein 
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Einflu® der Reynoldsschen Zahl besteht, ist dieser fiir den turbulenten Fall recht betrachtlich. 
Man kann zeigen, da man den bei E. Gruschwitz, A. Kehl und H. C. Garner auftretenden 
Gleichungen zur Berechnung des Formparameters die Form unserer Gleichung (48) geben kann!. 
Die daraus folgenden Zusammenhinge fiir H(L) und K(L) haben wir in oben bezeichneten Ab- 
bildungen miteingetragen. Die Unterschiede zwischen den einzelnen Verfahren beruhen also auf 
den Abweichungen der Kurven H(L) und K(L). 

Zur Lésung von (48) machen wir 
fiir K(L) einen linearen Ansatz (vgl. 


hierzu die Abb. 10 und 11): 
K(L)= a(L— b). (49) 
Die GréBen a und b ergeben sich 


etwa nach Tabelle 2. Im turbulenten 
Fall ist b noch von der Reynolds- 
schen Zahl und damit von der Lauf- 


lange x abhangig. 


Truckenbroat 
| SC  — — Garner 


-0,02 


— -— Gruschwitz 


2 a lf ei. | re AES Be lal == | re 
-007 Dy a O07 G02 G03 G04 =-02 =O7 0 40 OF 
Abb. 8. Zusammenhang zwischen dem Reibungs- Abb. 9. Zusammenhang zwischen dem Reibungsschichtdicken- 
schichtdickenverhaltnis H, sowie dem Wandschub- verhaltnis H und dem Formparameter L bei turbulenter 
spannungsbeiwert & und dem Formparameter L Strémung. 


bei laminarer Strémung. 


Setzen wir (49) in (48) ein, so erhalten wir eine lineare Differentialgleichung erster Ordnung 
fiir L, die wir geschlossen lésen kénnen. Zu diesem Zweck machen wir die Substitution 


dx fe 


poe 2 (G+ Af Ure RiNaa\2 (50) 


Die letzte Beziehung folgt aus der Tatsache, daB in der Formel (36) fiir 0 der Nenner gerade 
die mit 1/A multiplizierte Ableitung des Zahlers ist. Die Zahlenwerte a/A fiir glatte Flachen 
sind in Tabelle 2 miteingetragen. 

Ohne im einzelnen auf die Zwischenrechnung einzugehen, bei der wir noch eine partielle Inte- 


gration vornehmen, derart, daB die Ableitung dU/d& nicht mehr vorkommt, finden wir schlieB- 
lich 


é 
B U eile 
L=*L1,+In o ' : [eeo—m SP] (51) 


Im Anfangspunkt x, = 0 ist auch & = 0. Es verschwindet dann das erste Glied in (51). Speziell 
fiir den laminaren Fall gilt mit b= 0 der einfache Ausdruck 


i [a Tats (52) 


+ Ein ausfihrlicher Vergleich der einzelnen Verfahren untereinander ist in einem unveréffentlichten 
Bericht des Verfassers durchgefiihrt worden. 
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Uber die Anfangswerte fiir L im Staupunkt x,;== 0 sowie im Umschlagspunkt wird noch be- 
richtet werden. 


Wie man aus (51) ersieht, kann die neue Verdnderliche & mit einem beliebigen Faktor 
versehen werden, ohne da sich dadurch L andert. Wir kénnen dafiir also auch schreiben 


iia x/1 

i a yl ae 3+2n R l+n x! 
ar (3 biti decoah eay ear 
te x,/T 


Tabelle 2. Die GréBen a, b und a/A bei K 


laminarer und turbulenter Strémung. 


= (50a) 


laminar 


eet L<o? | turbulent 
a 2,87 3209 3,04-10-? 
U é 
b 0 0,07 1g —— = 0,23 
6.5 
A > 8.0 4,0 
1 Druckabfall; 2 Druckanstieg. 


Die Berechnung von € gestaltet sich also 
sehr einfach, da der Ausdruck in der 
eckigen Klammer bereits bei der Be- 
rechnung der Impulsverlustdicke (41) vor- 
kommt. 


——— Naherung (49) 


Eine vollstandige Rechnung geht nun =G7. J 
folgendermafen vor sich: Abb. 10. Die Funktion K(L) bei laminarer Strémung. 
Gegeben: 


k- 107 


Allgemeine GréBen 
U(x), R(x), =~. 


: a , y 
Strémungszustand n, cr, =. TE mae 


(laminar oder turbulent) 


Anfangswerte U,, R,, #,, L,- 


Gesucht: 
Impulsverlustdicke (x), Ud\_ 


Formparameter L, H(x). 


Zuerst berechnet man nach (41) 
durch Ausfithrung einer ein- wo'{ 
fachen Quadratur die Impuls- 


Truckenbrodf | 


verlustdicke. Hieraus  bildet ee 
man, sofern man _ turbulente ee Le aie 
Strémungen zu berechnen hat, SPS AOE ANS 
die Reynoldssche Zahl U@d/v, — 194 

mit der man dann nach Ta- 3 ae 

belle2 die GréBe b berechnet. | -140 4 


Als nachstes bestimmt man Abb. 11, Die Funktion K(L) bei turbulenter Strémung. 

nach (50 a) die neue Verander- 

liche £, tiber die man die Funktion b—InU/U, integriert. Nach (51) erhalt man dann den 
Formparameter L und mittels Abb. 8 bzw. 9 das Reibungsschichtdickenverhalinis H. Wegen 
des groBen Exponenten a/A in (50a) wachst die GréBe € sehr rasch mit grifer werdender 
Lauflange x. Das bedeutet aber, das das Glied (&/£) L, in (51) mit zunehmender Lauflange 
immer mehr an Bedeutung verliert, womit gezeigt ist, da® der Anfangswert L, bei gréBeren 


Abstanden x (z.B. in der Nahe des Ablésungspunktes) nur geringen EinfluB hat. Der Vorteil 
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unseres Verfahrens liegt nun darin, daf jeweils nur einfache Quadraturen auszufiithren sind. 
Dariiber hinaus kommen von den Ausgangswerten U und R keinerlei Ableitungen nach x vor. 


c) Die Anfangswerte, Die Anfangswerte sowohl fiir die Impulsverlustdicke # als auch den 
Formparameter L kénnen je nach dem vorliegenden Fall verschieden sein. 
1. Anstrémen eines Kiérpers mit Staupunkt. Die Reibungsschichtrechnung wird 
im Staupunkt, x,— 0, begonnen. Dort ist U=0 und R= 0. Aus (33) folgt dann sofort 
C= 0: (53) 


Sowohl im ebenen als auch rotationssymmetrischen Fall andert sich die Potentialgeschwindig- 
keit linear mit der Lauflange, d. h. 


Ui ex mit ce (ae 5 (54) 
0 


Die Konstante c ist fiir den ebenen und rotationssymmetrischen Fall voneinander verschieden. 
Fiir die Anderung des Halbmessers ergibt sich ebenfalls eine Linearitat in x, namlich 


Re= x. (55) 

Die Integrationskonstante in (34), (36) und (41) verschwindet : 
ae = Oe (56) 
Setzen wir (54) und (55) in (41) ein, dann folgt, wenn wir alle GréBen, die im rotationssymmetri- 


schen Fall zusatzlich in die Formeln eingehen, in geschweifte Klammern {} setzen, nach kurzer 
Zwischenrechnung : 


1—n 
n>: F= : es ae (57) 
2[2(2+n)+ (1+n}]itrerl*” 
Hierin bedeutet c* = i = ; den dimensionslosen Ausdruck fir c. Im laminaren Fall (n = 1) 
ergibt sich ein von Null verschiedener Wert fiir die Impulsverlustdicke in Staupunktsnahe 
c 
(7), "era ee a 


1,328 


l 
Wenn man beachtet, da c* = 7° und cr ——==— ist, kann man fiir (58) auch noch schreiben 


co pus i 
Vice 
£2 I 
V2 = 1,328!° (2) — (59) 
i“ na OF ae aes 
Im turbulenten Fall (n = 1/6) hat die Impulsverlustdicke im Staupunkt selbst (x — 0) den Wert 
Null. 
Fiir das Impulsverlustdickenverhaltnis im Staupunkt bei gleichem Geschwindigkeitsanstieg 
(dU/dx)) im rotationssymmetrischen und ebenen Fall (c= konst) gilt dann 


1 


8 2(2-+ n)\i* 

Do eb c=konst Sie =| j oy 
Setzen wir die Zahlenwerte fiir n ein, dann ergibt sich im laminaren Fall der Wert 0,867 und im 
turbulenten Fall 0,816. Der exakte Wert fiir die laminare Strémung, den man aus dem Hartree- 
Profil errechnen kann, betragt 0,845. 


Den Anfangswert des Formparameters L bestimmen wir nach (48), wenn wir darin nach (53) 
0, = 0 setzen, 


U heed (UA (61) 
Den Wert I[% erhalten wir aus der folgenden Grenzwertbestimmung : 
Peli sole z (62) 


xo U dx cee 22) tk ey 
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Im laminaren Fall ergibt sich mit (59) 


0,22 
DoS OS sera (62 a) 
und aus (49) 
= Ae) To = 0,077 2 
Ly a a 3+ {1} (63) 


(siehe Tabelle 3). Zwar ergibt sich nach (62) im turbulenten Fall ein endlicher von Null ver- 
schiedener Wert fiir I"), der uns aber nicht weiterhilft, den Wert Ly nach (61) zu berechnen, 


da uns die Werte K(L) fiir die Reynoldsschen Zahlen hee 0 nicht bekannt sind. Da die Stré- 
‘ yp 

mung im Staupunkt immer erst laminar sein wird, verweisen wir auf Abschnitt 4,c 3, wo wir 

tiber die Anfangswerte im Umschlagspunkt berichten. 


Fiir die Anderung des Formparameters L mit der Lauflange x im Staupunkt gilt noch, wie 
man aus (48) herleiten kann, 


2. Einstrémenineinen Kanal. Die Reibungsschicht beginnt in einem Punkt x, = x,—0, 
wo die Geschwindigkeit, und ebenso auch der Halbmesser im rotationssyinmetrischen Fall, eine 
endliche GréBe haben. Da wieder 7, = C, = Cf = 0, folgt sofort aus (34), (36) und (41) 


O;, = 0 und My, = 0 5 (65) 
und aus (44), wenn dU/dx + oc, 
i 0. (66) 
Der Formparameter folgt dann aus (48) zu 
K(L,)=0. (67) 
Im laminaren Fall wird daraus 
0s (68) 


d. h. der Wert der langsangestrémten ebenen Platte. Es besteht kein Unterschied zwischen dem 
ebenen und rotationssymmetrischen Fall (siehe Tabelle 3). Fiir die turbulente Strémung gilt 
das bereits bei der Staupunktstrémung Gesagte. 

Die Anderung des Formparameters L mit x erhalt man aus (48) zu 


dg} il fans a\U dx 
1+4(z), geet 


Der Zahlenwert a/A ist aus Tabelle 2 zu entnehmen. 


(32), Gran re € 2) (69) 


Tabelle 3. Laminare Anfangswerte. Die exakten Werte sind, 
sofern sie nicht mit den Naherungswerten iibereinstimmen, in 
Klammern gesetzt. 


3. Umschlag laminar-turbu- 
lent. Nach einem gewissen Ubergangs- 


bereich geht die laminare Reibungs- 
schicht an der Stelle x,— x, in die 


Staupunktstrémung 


Kanalstrémung 


rotations- eben, rotations- 
turbulente Schicht iiber. Aus der Theo- eben symmetrisch symmetrisch 
rie der Turbulenzentstehung, vgl. H. 
Schlichting?, kann man in Ubereinstim- oe 6 0,271 0,235 0 
mung mit Messungen entnehmen, daB \ : (0,292) (0,247) ‘ 
der Umschlag an Stellen liegt, die gegen- 0,0260 0,0195 rs i 
iiber dem Geschwindigkeitsmaximum Ly (0,0292) (0,0208) . 
etwas stromabwarts liegen. Da die Vor- 2,25 2,32 H es 
ginge im Ubergangsgebiet noch nicht Hy (2,22) (2,30) if 3 
erforscht sind, ist die Bestimmung der 0,345 0,320 
fiir die turbulente Rechnung benétigten Mo (0,360) (0,324) “Ys Goee0 


Anfangswerte fiir ® und H bzw. L nur 


annahernd méglich. Das Grundsatzliche unserer Uberlegungen ist in Abb. 12 am Beispiel der 
ebenen Platte gezeigt. Bis zur Stelle x, ist die Reibungsschicht laminar und folgt den Gesetz- 


1 Man kann zeigen, da dieser Wert auch aus (52) folgt. 
2 H. Schlichting, Grenzschicht-Theorie, Karlsruhe, 1951. 


224 Truckenbrodt: Ein Quadraturverfahren zur Berechnung der Reibungsschicht. _ Ingenieur-Archiv 


maBigkeiten nach Blasius, d.h. 0 ~ /%. Der entsprechende H-Wert ist konstant und betragt 
H — 2,60. Von der Stelle x, ab ist die Reibungsschicht voll-turbulent. Hierfiir gilt dann nahe- 
rungsweise ? ~ #,-+ ¢ (x — x,). Der entsprechende H-Wert hangt noch etwas von der im 
Umschlagspunkt herrschenden Reynoldsschen Zahl U.0/y ab. Nach Messungen, z. B. von 
F. Schultz-Grunow 1, sowie aus den ahnlichen Lésungen bei Gleichdruck von J. Rotta 2 findet man 
etwa 1,2 < H < 1,4, und zwar ist der H-Wert umso kleiner, je gréBer die Reynoldssche Zahl 
ist. Der H-Wert fallt also gegeniiber dem laminaren Wert ab. Er muB sich also im Ubergangs- 
bereich stetig von dem laminaren auf den turbulenten Wert andern. In Abb. 13 haben wir die 
Differenz AH iitber UD/y aufgetragen. Der Einfachheit halber lassen wir x, mit %; zusammen- 


40 


Schultz - Grunowlebene Platte) 


ae ahniliche Losun 
a Rotta be/ a) 


05 


yy, 
Dv 


eee Se a 5 2 


Diy, 10 70* 10° 
Abb. 12. Ubersichtsskizze vom Umschlagspunkt Abb. 13. Die Anderung des Formparameters H im 
der ebenen Platte. Ubergangsgebiet von laminarer zu turbulenter Strémung. 


fallen, dann ist der Anfangswert fiir die Impulsverlustdicke der turbulenten Rechnung gleich 
derjenigen Impulsverlustdicke, die sich im Punkt x, ergeben wiirde, wenn die Strémung bis zu 
diesem Punkt voll-laminar ware. Wir wollen daher setzen 


0; (Xp) = — 0 (2) (70) 
Fir den H-Wert schreiben wir 
iH, (4) = = (x5) A (71) 


Wir wollen nun annehmen, daB AH auch fiir die Falle mit Druckgradient aus Abb. 13 ent- 
nommen werden kann. Diese Annahme scheint uns berechtigt, da der Umschlagspunkt in der 
Nahe des Punktes verschwindenden Druckgradienten liegt und hierfiir in brauchbarer Naherung 
die Werte der ebenen Platte angenommen werden kénnen. Haben wir so H,, bestimmt, so er- 
mitteln wir nach Abb. 9 den Wert L, = L(x,), der dann in (51) als L, eingeht. Es sei an dieser 
Stelle noch einmal darauf hingewiesen, da$ der Anfangswert fiir L an Stellen, die weiter strom- 
abwarts vom Umschlagspunkt liegen, nur untergeordnete Bedeutung besitzt, so daB® uns eine 
etwas rohere Abschatzung berechtigt erscheint, wenn es darum geht, den Ablésungspunkt 
zu bestimmen. 


d) Der Ablésungspunkt. Wichtig ist noch die Kenntnis der Lage des Ablésungspunktes. 
Er ergibt sich, wenn die Wandschubspannung den Wert Null annimmt. Wahrend dieser Punkt 
im laminaren Fall entsprechend den vorgegebenen Geschwindigkeitsprofilen durch «= 0 nach 
(19) festliegt, ist man im turbulenten Fall noch nicht in der Lage, eine ganz klare Angabe iiber 
den Ablésungspunkt zu machen. Nach den Wandschub- 
spannungsansatzen von Ludwieg-Tillmann und Rotta 
(vgl. Abb. 3) verkleinert sich die Schubspannung mit 
wachsendem Wert H, ohne jedoch den Wert Null zu er- 
reichen. Als Naherungsregel kann nach v. Doenhoff- 4.038 ° 
Tetervin® etwa gelten, daB die Ablésung frithestens be- 3 —0.018 


Tabelle 4. Werte des Formparameters 
im Ablésungspunkt. 


turbulent 


1,8 bis 2,4 
—0,13 bis —0,18 


1 F. Schultz-Grunow, Luftfahrtforschung, 17 (1940), S. 239. — ? Siehe FuBnote 4 auf Seite 212. — 
3 Siehe FufBnote 11 auf Seite 211. — 4 Diese Angabe entspricht etwa der Gruschwitzschen Annahme 
n = 1 — (us/U)? = 0,8. 
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ginnt, wenn H =~ 1,84 ist, und bestimmt stattgefunden hat, wenn H den Wert 2,4 erreicht hat. 
Tabelle 4 gibt eine Zusammenstellung der auftretenden Werte. 

e) Beispiele. Die Brauchbarkeit des oben beschriebenen Verfahrens wird an je einem Bei- 
spiel laminarer bzw. turbulenter Strémung gezeigt. 


1. Howarth-Strémung (laminar). Als Beispiel fiir die laminare Strémung wahlen wir 
die bekannte Howarth-Strémung !. In diesem Fall ist die Geschwindigkeitsverteilung 


=} 9 


Setzen wir (72) in (41) ein, worin wir die Radiusverteilung R streichen, dann kénnen wir die 
Quadratur geschlossen ausfithren. Es ergibt sich fiir die Impulsverlustdicke 


edaey \2 wlth: al 0,2711 1 = ee = 


T 9V6 ( uA) a oe ORS rs (73) 


(eee, 


—— exakte Werte nach Howarth 


Abb. 14. Die laminaren ReibungsschichtgréBen der Howarth-Strémung. 


In Abb. 14 zeigen wir den Vergleich dieser Naherung mit den exakten Werten von Howarth. Die 
Ubereinstimmung ist im Rahmen der Zeichengenauigkeit vollkommen. 
Zur Bestimmung des Formparameters berechnen wir zunidchst die neue Veranderliche é 


nach (50 a) 
7s init chem He al a a Cake (74) 
Auch der Formparameter selbst 148t sich nach (52) analytisch berechnen. Da wir uns ent- 


sprechend der vorgegebenen Geschwindigkeitsverteilung im Gebiet des Druckanstiegs befinden, 
wahlen wir laut Tabelle 2 den Wert a/A = 8. Man findet nach kurzer Zwischenrechnung 


é Z oo 
Pep eitalrcn sy. 8 1 LA eS zm 
Lae {og e'= sm {ro7 Paes (75) 
0 


0 


1 L. Howarth, Proc. Roy. Soc. London, Ser. A 919, Bd. 164, (1938) S. 547. 
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Abb. 15. Die turbulente Reibungsschicht am Profil NACA 65 (216)—222 
&% = 10,1°, Ups = 2,64. 10°, Messung nach NACA Rep. No. 772. 
‘i 


Unter Beachtung von (74) 
haben wir L(x) ebenfalls in 
Abb. 14 dargestellt. Aus diesen 
Werten haben wir unter Zu- 
hilfenahme der Abb. 8 die 
Werte « fiir die Wandschub- 
spannung bestimmt. Sie sind 
auch in Abb.14 dargestellt und 
dort mit den exakten Werten 
von Howarth verglichen. Die 
Ubereinstimmung ist recht 
befriedigend. Gezeigt ist auch 
die Kurve, die sich nach dem 
Pohlhausen- Verfahren ergibt?. 
Man erhialt, wie schon A. Walz 
feststellt, eine | wesentlich 
bessere Bestimmung des Ab- 
lésungspunktes, «= 0, wenn 
man, wie bei unserem Ver- 
fahren, auBer dem Impulssatz 
noch den Energiesatz benutzt. 


2. Profil NACA 65(216) 
—222 (turbulent). Als Beispiel 
fir die turbulente Rechnung 
wahlen wir das Profil NACA 
65(216)—222 (approx), fiir 
das sowohl Messungen? als 
auch theoretische Rechnungen 
durchgefiihrt worden sind. 
Die Verhaltnisse beziehen sich 
auf die Profiloberseite des unter 
dem Anstellwinkel « = 10,1° 
angestellten Profiles. Die Rey- 
noldssche Zahl betragt U,|/v 
= 2,64-108 Die graphische 
Darstellung der Geschwindig- 
keitsverteilung, mpulsverlust- 
dicke und des Formparameters 
nach der Messung, sowie das 
Ergebnis unserer Rechnung ist 
in Abb. 15 mitgeteilt. Wir 
haben unsere Rechnung indem 
ersten gemessenen Punkt, d.h. 
x/L= 0,075 begonnen. Die 
Ubereinstimmung zwischen 
unserer Naherung und der 
Messung ist als gut zu bezeich- 
nen °, Zum Vergleich haben 
wir noch das Ergebnis nach 


1 Diese Kurve haben wir von 
A, Walz (FuBnote 5 auf Seite 211) 


entnommen, 


2 A. E. v. Doenhoffu. N.Teter- 
vin, FuBnote 11 auf Seite 211. 


3 Die groéBeren Unterschiede 
zwischen Messung und Rechnung 
fiir die Impulsverlustdicke in der 
Umgebung des Ablésungspunktes 
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der v. Doenhoff- Tetervinschen Rechnung, sowie nach den Verfahren von Garner, Gruschwitz und 
Gruschwitz-Kehl miteingetragen. Die letzten beiden Verfahren weichen fiir den FormparameterH 
von den anderen Verfahren und der Messung stark ab. 


5. Zusammenfassung, Zur Berechnung der charakteristischen ReibungsschichtgréBen (Im- 
pulsverlustdicke und Formparameter fiir das Geschwindigkeitsprofil) werden zwei Quadratur- 
formeln angegeben, die sowohl fiir den laminaren als auch den turbulenten Strémungszustand 
gelten. In diesen Formeln ist gleichermafen der ebene und rotationssymmetrische Fall enthalten. 

Die Berechnung der Impulsverlustdicke erfolgt durch eine einfache Integration des Energie- 
satzes. Die Gleichung fiir den Formparameter wird durch Kopplung des Impulssatzes mit dem 
Energiesatz gewonnen. Die Kenntnis der Ableitungen der Geschwindigkeitsverteilung und des 
Kérperradius nach der Lauflange ist nicht erforderlich. 


Anhang. Betrachtung iiker die Formel zur Berechnung der Impulsverlustdicke. 


Wie bereits erwahnt, wurde bisher die Impulsverlustdicke aus dem Impulssatz (16) berechnet. 
Unter gewissen Vereinfachungen ist es méglich, zur Berechnung der Impulsverlustdicke aus 
dem Impulssatz eine unserer Gleichung (36) ahnliche Forme! zu finden. Dabei gibt es zwei 
Méglichkeiten: 

1. Die Wandschubspannung wird nach den Gesetzen der ebenen langsangestrémten Platte 
bestimmt, in denen die Lauflinge x durch die Impulsverlustdicke # und die Anstrémungs- 
geschwindigkeit U, durch die értliche Geschwindigkeit ausgedriickt werden. Sowohl fiir den 
laminaren als auch turbulenten Fall kann man schreiben 


eV? (Zy" 


ae 


Setzt man auBerdem im Impulssatz fiir das Reibungsschichtdickenverhaltnis H einen konstan- 
ten Wert ein, dann kann man, wie E. Truckenbrodt in einem unveréffentlichten Bericht (vgl. 
H. Schlichting, Grenzschicht-Theorie 5.430) gezeigt hat, den Impulssatz geschlossen integrieren, 
und zwar ergibt sich 


GE nas Ut Ro dx! 


ONO ay : 
=| y a Ut Ritn yi 

wobei a = (1-+ n)(H + 2) — n ist. 
Die entsprechenden Zahlenwerte sind aus nachstehender Tabelle zu entnehmen. 


2. Nach K. Pohlhausen fiir den laminaren Fall und nach E. Buri fiir den turbulenten Fall 
sind die Wandschubspannung und das Reibungsschichtverhaltnis Funktionen der GréBe 


@ dU 
Di age PRG 
fiver ese Li = 
D. h. i= (vo und HEF (ry. 
vy 


Setzt man A = OR'+" und x* —f Rtn dx, dann kann man nach einiger Zwischenrechnung 
0 
die folgende Gleichung finden: 


dA - A dU 

aa PD) mit LS Fak 
Hierin stellt ®(1) = (1-+n) f,(’) — [2+n+ (1+ n)f,(L)] LF eine universelle Funktion fiir den 
laminaren bzw. turbulenten Zustand dar, die sich im ersten Fall analytisch berechnen 1aBt 
(z. B. fiir die Pohlhausen-Polynome) und sich im zweiten Fall aus Messungen bestimmen 1aBt. 
In beiden Fallen kann man naherungsweise einen linearen Zusammenhang von ® und ["finden 


(vgl. A. Walz, Lilienthal-Bericht 141 (1941) und E. Burt). 


x/l ~ 0,55 kénnen, worauf A. E. v. Doenhoff und N. Tetervin hinweisen, auf systematische MeSiungenauig- 
keiten zurickzufiihren sein. 


17 
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Setzt man diesen Ausdruck oben ein, dann 1aBt sich die Integration geschlossen ausfithren und 


liefert nach einfacher Zwischenrechnung 


C f UP Rt” dx! 
o= (2%) ieee ‘ 


y 
Die Zahlenwerte c und b sind ebenfalls in nachstehender Tabelle angegeben (vgl. H. Schlichting, 
Grenzschicht-Theorie S. 199 u. 424). Zum Vergleich sind auch die Konstanten nach (36) mit- 


geteilt worden. 


laminar turbulent 
| : i : Falk 
n 1 q (Blasius) % (Falkner) 
A 0,441 — 0,0076 
eee Ne al ad ijlistarives a OMA 0,0160 0.0076 
8 c 0,470 0,0160 = 
Exponent der 3+ 2n 5,0 — 3,33 
Geschwindig- a 8,2 4,0 3.67 
keitsverteilung b 5,0 4,0 | — 


(Eingegangen am 20. Dezember 1951). 


Anschrift des Verfassers: Dr.-Ing. E. Truckenbrodt, Braunschweig, KérnerstraBe 12. 
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Ein Stoffgesetzansatz fiir elastische, anisotrope Medien. 
Von E. Pieruschka. 


1, Einfiihrung. In einer friiheren Arbeit! wurde mit Hilfe des Tensorkalkiils ein Stoffgesetz 
fiir isotrope Korper aufgestellt. Nun soll die Untersuchung auf anisotrope verallgemeinert wer- 
den. Doch wird hier insofern eine Einschrankung gemacht, als nur der elastische Bereich behan- 
delt werden soll. Dieser diirfte aber geniigen, denn Uberbeanspruchung des anisotropen Kérpers 
fiihrt mindestens zur Zerstérung der Kristalle, d. h. der Anisotropie. 

Es sollen hier die gleichen Bezeichnungen wie in der ersten Arbeit gewahlt werden. | Der Ab- 
schnitt ,,Mathematische Hilfsmittel‘* mu aber noch erganzt werden. Das doppelte Multipli- 
kationszeichen (--) zwischen zwei Tensoren soll andeuten, da® die skalare Multiplikation zuerst 
auf die benachbarten und dann noch auf die verbleibenden Einheitsvektoren angewandt wird: 


YB = ayy by, ? (1.1) 


Neben der antimetrischen Deutung des Rotors eines Vektors, der nach der ersten Arbeit mit Rot 
bezeichnet wird, braucht man hier auch die vektorielle Deutung. Dieser Rotor soll mit rot be- 


zeichnet werden. Das Volumenintegral des Rotors eines Vektors 1aBt sich nach der bekannten 
Formel 


fav- ros = 6 [af3] (1.2) 


in ein Oberflachenintegral umwandeln. Die eckige Klammer deutet das vektorielle Produkt 
an. Weiter werden spater noch die Beziehungen gebraucht: 


div (uM) =%W-grad(u) -+u-div%, (1.3) 
div(2l-b)=%W--gradb +56- div, (1.4) 
AS = grad(div$) -div(Rot 8). (1.5) 


2. Das Stoffgesetz. Das Stoffgesetz von Hooke fiir elastische, isotrope Stoffe lautet 
3 
$B = 3 (K9)-e+6G. (2.1) 


Der Ansatz fiir ein Stoffgesetz anisotroper Kérper muf so erfolgen, daB das Hookesche Gesetz 
ein Sonderfall des gesuchten Stoffgesetzes ist. Dies 148t sich so erreichen, dafi man fiir den 
Kugeltensor (K $) einen allgemeinen symmetrischen Tensor setzt. 

Damit diirften aber die Méglichkeiten (2.1) zu verallgemeinern erschépft sein, denn im Sum- 
manden G- &,, der sich auch (G 4): Gy schreiben 1aBt, kann nicht der Kugeltensor G ¥ durch einen 
symmetrischen Stofftensor & ersetzt werden, da & - ©, nicht mehr symmetrisch ware: die Sym- 
metrie der Summanden der Tensoren miiBte aber bei einer Verallgemeinerung von (2.1) wegen 
der Symmetrie von ‘§ erhalten bleiben. ae 

Der einfachste Ansatz fiir ein Stoffgesetz anisotroper, elastischer Medien, aus dem sich bei 
Isotropie automatisch das Stoffgesetz von Hooke ergibt, ist demnach 
= 
ga Pa-e+eG. (2.2) 
Gist der neu eingefiihrte nur vom Stoff abhangige Tensor. Der Faktor von St ist so gewahlt, dal 
bei isotropen Stoffen mit § = k J die Gleichung (2.1) entsteht, wenn man K in (2.1) als Intensitat* 


deutet. Es ist dann namlich K —k- V3 mit k als Mittelkomponente von St. 
Die Aufteilung von (2.2) in Kugeltensor und Deviator ergibt 


p= Bree, (2.2 a) 


1 Ing.-Arch. 19 (1951) S. 271. ; 
= Die te aes @ des Tensors Q{, die durch A? = aababa gemafB (2.4) a. a. O. gegeben ist, lakt 
sich hiernach schreiben -- = A’. 


ila 
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p= 82, +66, (2.2b) 
Formel (2.2) kann auch nach € aufgelést werden. Man erhalt 
‘ae lalels. 1 1 
g=1 (785 _1o)p+en- (2.3) 
Die anisotropen Stoffe, die dem hier entwickelten Stoffgesetz gehorchen, haben zwei wesentliche 
Kigenschaften. 
Ist der Spannungstensor ein reiner Kugeltensor — Deviator ‘% ist Null — so ist der Ver- 


zerrungstensor ein allgemeiner symmetrischer Tensor. Er hangt linear vom Druck p ab. Seine 
Hauptachsen ergeben sich aus der Struktur des Stoffes. Hier wird also ein vor der Beanspru- 
chung kugelférmiges Stoffgebiet, das bei Isotropic auch bei Belastung seine Kugelgestalt behalt, 
zum Ellipsoid, dessen Hauptachsen vom Stoff allein abhangen. Sicher wird wohl diese erste 
Kigenschaft auf alle anisotropen Stoffe im elastischen Bereich zutreffen. 

Ist die Mittelkomponente des Spannungs- und damit auch des Verzerrungstensors Null, dann 
sind nach (2.2a) und (2.2b) der Spannungs- und der Verzerrungstensor — es sind hier Devia- 
toren —koaxial und ahnlich und hangen linear voneinander ab. In diesem Spezialfalle unter- 
scheiden sich nach dem hier gemachten Ansatz isotrope und anisotrope Stoffe nicht. Es muB 
also ein anisotroper Kérper, der dem Stoffgesetz (2.2) gehorcht, bei der Beanspruchung durch 
ein reines Drehmoment, immer dieselbe Verformung haben, unabhangig von der Lage der Dreh- 
momentachsen im Tensorfeld §. Es mu durch Versuche festgestellt werden, welche aniso- 
trope Stoffe die zweite Eigenschaft zeigen. 

Der Spannungs- und der Verzerrungstensor bestehen aus je sechs Komponenten. Der all- 
gemeinste Fall eines Stoffgesetzes wiirde demnach 36 Stoffkonstanten enthalten. Bei Isotropie 
des Stoffes und Elastizitat reduzierten sich diese auf zwei (K und G bzw. © und m). Das aniso- 
trope Stoffgesetz (2.2) hat 7 unabhingige Stoffkonstanten, namlich die sechs Komponenten des 
Stofftensors § und den Schubmodell G. Man kann dies auch so ausdriicken, das neben den drei 
Hauptrichtungen des Tensors §& vier Stoffkonstanten vorhanden sind. Also ist auch hier eine 
sehr starke Spezialisierung des allgemeinst méglichen Stoffgesetzes vorhanden. 

Das Wesentliche des Stoffgesetzes (2.2) kommt am besten zum Ausdruck, wenn man den Spe- 
zialfall betrachtet, bei dem die Hauptachsen des Spannungstensor mit dem des Stofftensors zu- 
sammenfallen. Legt man auch die Koordinatenachsen in die Hauptachsen von § und damit auch 
von S$ und ©, so ergibt sich 

1 


2 
oo ia G Ps 3 (; ae a) : (2.4) 


Bezeichnet man die Koordinatenachsen mit den itiblichen Buchstaben x, y und z, die Haupt- 
dehnungen mit ¢,, ¢, und ¢,, die Hauptspannungen mit o,, ¢, und o,, und die Komponenten des 
Stofftensors § mit k,, ky und k,, und fiihrt man ahnliche GréBen wie den Elastizitatsmodul und 


die Poissonsche Zahl bei isotropen Stoffen ein, so erhalt man . 
es! Oy + oz 
ex Ex (0. Mx f 
1 Oz + Ox 
eg i (% seo (2.5) 


G 
Gein Pe aru 
cele :) 
G 
cs ii oy (2.6) 
3k 3/3 ar 
— 2 
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G 
Me GY yp tty = a G0 ™ Boe Cy 
Fir die hier angenommene Belastung erhilt man die drei Elastizitatsmoduln E,, Ey und E, und 
die drei Poissonschen Konstanten m,, m, und m,, die nicht unabhangig voneinander sind, sondern 
sich durch die Elastizitatemoduln und dem Schubmodul G als vierte Konstante ausdriicken lassen. 
Das Gleichungssystem (2.5) zeigt eine ahnliche Form wie das bekannte entsprechende Gleichungs- 
system bei isotropen Stoffen. 


3. Die Wellengleichung anisotroper Stoffe. Als Beispiel soll hier wieder wie in der ersten be- 
reits erwahnten Arbeit die Wellengleichung anisotroper Stoffe aufgestellt und gelist werden. 
Es werden hier die Voraussetzungen wie bei der ersten Arbeit gemacht und aus dieser Arbeit die 
Beziehungen itbernommen!; 

075 


div $ = Car? (3.1) 
div © = A8 +; grad (div 8). (3.2) 

Forme] (3.2) 1a8t sich durch Einsetzen von (1.5) auch schreiben 
div 6, = grad (div 8) + div (Rot 3). (3.3) 


Entwickelt man div (eS) nach (1.3) und beachtet, daB div § = 0 ist und e = = div 8, so erhalt 
man 


div(e 8) = 2 K- grad (div. 8). (3.4) 


Nun bilde man die Divergenz auf (2.2) und setze in der neu entstandenen Gleichung (3.1), (3.4) 
und (3.3) ein. So erhalt man die Wellengleichung der anisotropen Stoffe, die dem Gesetz (2.2) ge- 
horchen: i 


22 — 
05s = [3 & grad (div 8) ++ 6 grad (div 8) + G div (Rot 8). (3.5) 
Forme] (3.5) laBt sich noch weiter behandeln. Dazu zerlegt man den Vektor § so in die zwei Vek- 
toren 8 und 8, daB® der Rotor vom ersten und die Divergenz vom zweiten Vektor Null werden. 
Es wird so ts) the 
8=843 (3.6) 
mit Rot8=0 und div3=0, (3.7) 


Da 8 und 8 nicht linear voneinander abhiingen, wie sich noch spater zeigen wird, zerfallt die Wel- 
lengleichung (3.5) in die zwei Gleichungen 


023 LE Sis tem An : ae) 

#1 (13. 2 +46 3) grad (divs), (3.8 a) 
as G,. A 

3 © AL 3.8b 
aii ; div (Rot 3) ( ) 


Formel (3.8a) 14Bt sich noch umformen. Dazu bilde man die Divergenz letzterer Gleichung 
und beachtet, da 

= nrg 1 

ist. Unter Verwendung dieser Beziehung und von (1.4), erhalt man schlieBlich, da die Divergenz 


des konstanten Tensors (13 K +26 3] Null ist, 


eee 384269 ++ grad (grad p) . (3.8 c) 
Ot? 0 3 

Damit ist die Wellengleichung (3.5) in die Gleichung fiir die Druck- oder Kompressionswelle 
_(3.8a) bzw. (3.8c) und die Scherungswelle (3.8b) zerlegt. 


1 Vel. die dortigen Gleichungen (9.4) u. (9.3). 
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4, Die Kompressionswelle. Legt man die Koordinatenachsen in die Hauptachsen von KR, so 
da® k,;—= 0 fiir « = B wird, so kann man (3.8c) schreiben 


es, Of # 3 kaa\ Op re 
fe ak: +3 Goa Sa 
Setzt man 
k 
2 af EO, He 
so erhadlt man aus (4.1) 
a ml 2 (4.3) 
Cl a6 Over 


Liegt das Stérungszentrum mit der Frequenz « im Koordinatenursprung und kénnen sich die 
Wellen frei nach allen Richtungen ausbreiten, so lautet die Lésung der Wellengleichung (4.3) 
— wiiren die y, reine Raumkoordinaten, so kénnte man von Kugelwellen sprechen — 


p diss oe a(t—R/v)) (4.4) 
mit 
x2 
i (4.5) 
3 +13 kaa/G 
und 


vy == ye ; (4.6) 


Um die Form der Wellenfronten zu erkennen, setze man in (4.4) und (4.5) R = Konstant. Man 
erkennt dann, daf die Wellenfronten auf den Ellipsoidflachen liegen, deren Gleichungen (4.5) 
angibt. Man kann also hier von Ellipsoidwellen sprechen — analog Kugelwellen. Die Welle 
breitet sich strahlenférmig um den Koordinatenursprung aus. Die Fortpflanzungsgeschwindig- 
keiten in den Hauptachsen von § sind 


ps =) Te (4 + 13 kaa) ae (4.7) 


Die Geschwindigkeit in einer anderen Richtung gibt die Lange des vom Ursprung nach der frag- 
lichen Richtung gezogenen Vektors t. Der Endpunkt des Vektors r liegt auf dem Ellipsoid 


2 
Xo 


= 
olg+ls kaa/6) 


=1. (4.8) 


Die Richtung von 3 ergibt sich aus der Formel (1.2). Statt 3 setze man hier 3. Da der Rotor von 
3 Null ist, erhalt man 


h [df 3] =0. (4.9) 
Betrachtet man eine Kugelflache im y,-System, so hat df die Richtung von tr. Aus Symmetrie- 


griinden muf fiir jeden Punkt der Kugeloberflache das vektorielle Produkt aus df bzw. t und 3 
den gleichen absoluten Wert haben. 


Deshalb kann (4.9) nur erfiillt sein, wenn 8 in Richtung von t liegt, d. h. senkrecht zu der Tan- 
gentialebene der Kugel steht. Dies bedeutet, da 3 im x,-System senkrecht zu der Tangential- 
ebene der Ellipsoide steht, auf denen die Wellen@outens liegen. 


5. Die Scherungswelle. Da die Divergenz von 3 Null ist, erhalt man aus (1.5) 
div (Rot 8) = 48. (5.1) 
Die Scherungswellengleichung nimmt damit die Form an 


23 G = 
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Eine Lésung von (5.2) ist die Kugelwelle 


Sga€a io — 
Sesto i a(t + t,—r/v) (5.3) 
r 


= 


mit 
= 
1 = Xe und t= = : (5.4) 
e 
Die Richtung von 8 ergibt sich aus dem Gaufschen Satz, erweitert auf einen Tensor zweiter Stufe 


(siehe die erste Arbeit, auf die hier Bezug genommen wird): 


v 


[ avs Gary (5.5) 


(Als Volumen nimmt man hier zweckmaBig eine Kugel). Da nun div 3 Null ist und der Wert 
von $-t auf der Kugeloberflache aus Symmetriegriinden immer gleich sein muB, folgt 
3 steht senkrecht auf r. 


(Eingegangen am 9. Januar 1952.) 
Anschrift des Verfassers: Dr, E, Pieruschka, Kassel-Wilhelmshohe, Bergstr. 15. 
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Zur Kombination thermischer und elektromagnetischer Felder 
im Falle der ebenen Platte. 


Von M. Kriiger. 


1. Einleitung. In der Warmetechnik finden in steigendem Mafe Verfahren Anwendung, welche 
die Erwairmung von Leitern und Nichtleitern durch hochfrequente Wechselstréme benutzen. 
Der Hauptvorteil dieser Verfahren, abgesehen von einer Reihe weiterer Vorteile rein technischer 
Art, liegt in der Warmeproduktion im Materialinnern. In allgemeiner physikalischer Formulie- 
rung lautet dabei die Aufgabenstellung: Es sind thermische und elektromagnetische Felder so 
zu kombinieren, da®B zu einer vorgegebenen Zeit im Material ein vorgegebener Temperatur- 
verlauf erreicht wird, meist noch unter der einschrankenden Bedingung, daB nirgends vor der 
Endzeit die Endtemperatur iiberschritten wird. Unter thermischen Feldern sollen dabei solche 
verstanden werden, wie sie sich bei Vorhandensein der sogenannten auBeren und inneren Warme- 
leitung ausbilden. Zu dem iiblichen inhomogenen Warmeleitungsproblem mit Anfangs- und 
Randbedingungen tritt also hier noch eine ,,Endbedingung* 
mit gegebenenfalls weiteren Einschrankungen. Im Interesse 
eines méglichst einfachen Einblicks in den Problemkreis mégen 
hier die elektromagnetischen Felder zunachst auf quasistationare 
Wechselfelder in dielektrischen Stoffen mit nicht zu groBem 
Verlustfaktor eingeschrankt werden uud das Material die Form 
einer ebenen, unendlich ausgedehnten Platte endlicher Dicke d 
Abb. 1. Ebene Platte der Dicke d. haben (Abb. 1). 


2. Allgemeine Lésung des Problems bei beliebiger Warmeproduktion. Es handelt sich darum, 
unter den eben gemachten Einschrankungen eine allgemeine Lésung des Warmeleitungspro- 
blems méglichst in geschlossener Form zu finden, soda daraus eine Bedingungsgleichung fiir 
den vorgeschriebenen Temperaturverlauf zu einer vorgeschriebenen Zeit gewonnen werden kann. 
Ks ist hier noch nicht notwendig, die Funktion f(z), welche, mit der Dichte und der spezifischen 
Warme multipliziert, die Warmeproduktion je ccm und sec im Innern der Platte darstellt, zu 
spezialisieren. Bezeichnet man mit ¢t die Zeit und mit (z, t) die Temperatur in ° C, so lautet also 
die zugrunde zu legende Warmeleitungsgleichung 

ao ab 

oF at SP + pe), (1) 
worin a, unabhangig von Ort, Druck und Temperatur angenommen, die Temperaturleitzahl 
bedeutet. Zur Zeit t = 0 herrsche die Temperaturverteilung 


Hz, 0) = yz) . (2) 
und an den Randern werde entsprechend dem Newtonschen Abkihlungsgesetz gefordert: 
pre on 
[P,(t) — H(0, t)] a, = — x ao (0, t), (3) 
[95() Odea) np = 2, ian) 4 
2 d Oz 9 ? ( ) 


worin ((t), Jy(t) zeitliche Temperaturverteilungen an den duBeren Randern der Platte bedeuten 
und x sowie M%, Aq die konstant angenommenen Koeffizienten der inneren bzw. duBeren Warme- 
leitung darstellen. Mit 
so te __& 
one eer Oy (5) 
lauten die Randbedingungen endgiiltig 


g(t) — (0, 1) = 798" (0,1), (6) 


Bult) — Hd, 1) = + 452 (dt) (7) 
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% = 90 bedeutet einen bei z = 0 vorgeschriebenen oder gesuchten (siehe spater) Temperatur- 
verlauf, x) = 00 liefert die adiabatische Randbedingung. Fiir x, gilt Entsprechendes. Sind dann 
f(z) und gz) fiir 0<z<d stetige, #,(t), Du(t) stiickweise glatte Funktionen fiir t>0, wobei die 
letzten beiden tiberdies noch der Bedingung geniigen sollen, daf sich ein a, finden 1aBt derart 
daB fiir alle «>a, die Integrale 


f Bo(t) e—% t dt, I Du(t) e+ dt (8) 


absolut konvergieren, so findet man mit Hilfe der Laplace-Transformation die eindeutig bestimmte 
Lésung!, die weiter unten verifiziert werden soll: 


D(z, t) 
¢ d 
1 eV et dy E, (z) 2 S 1 ‘ evet ay E, (z) a 
— a NES, a) HS | aes = —xEéE Jé 
2 mi | ee Ok Lr san n? (oe? ¥# 6) 1) [se dé 
L,+L_, 0 nial F 
a? aay = CMa ° e  (4—2) 
ee — 2 art ped ta) S50 —x? a 9 
ni od ey 1) 1 tmx Bo(T) e dt+5 cere | Oa(t) e dr (9) 
Deh, 0 Fs 
2 f 
1 4 
22 | dé, [760 a. 
0 0 
Hier ist 
spon, = ee NSS 
Ez) =e Wee aee e (10) 


gesetzt, und [ist das Integral iiber einen Integrationsweg der komplexen x-Ebene, welcher durch 


Pe Abbildune 
y=@r?, y=atip, a>a>0 (11) 


aus der Parallelen im Abstande « zur imagindren Achse der y-Ebene hervorgeht. Diese Parallele 
ist in Richtung von negativen zu positiven Imaginarteilen durchlaufen zu denken. Integral ih 
bedeutet das Integral iiber einen Weg der x-Ebene, welcher aus L, durch Spiegelung am Null- 
punkt der x-Ebene hervorgeht, also die entgegengesetzte Durchlaufrichtung wie L, hat. Ferner 
ist 
_ (l=) (1a) 
eae (A En (2) 

Um zu sehen, ob die Liésung (9) der Warmeleitungsgleichung (1) geniigt, hat man nur zu be- 
achten, daB t und zin den Exponenten von e-Potenzen in den ersten drei Summanden von (9) 
derart auftreten, daB diese Summanden nach den (1) entsprechenden Differentiationen, in (1) 
eingesetzt, verschwinden, wahrend der zweimal nach z ditferenzierte vierte Summand von (9) 
sich gegen f(z) in (1) weghebt. 

Zum Nachweis der Erfillung der Anfangsbedingung (2) setze man im ersten Summanden 
von (9) t = 0 und betrachte 


d 
< em ie(u) dx ard te ee Ss “Eg 
I ed acme -— e ) { me eras. (13) 


e274 oy) 1 1—xo% 


4 Hw 


Setzt man hier x = x, +i; ,,,; reell, so entsteht daraus: 


d 
eee a ih ; 
end IE aoa ea | p(é) e%§ d& (ew Dieses a Seger nd) ze , 


en Foe) 1— x) % 
Liy is é \ i 
1 Die Existenz und Eindeutigkeit der Lésung sowie ihre Herleitung mit Hilfe der Laplace-Transfor- 


mation ergibt sich nach Hilbert-Courant, Methoden der math. Phys., Bd. II, 5. 205ff., Berlin 1937, wenn 
man den Randbedingungen (6), (7), entsprechende etwas allgemeinere zulassige Voraussetzungen macht. 
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wobei x° der kleinste Wert von x, auf L,,, ist, welcher nach (11) sicher gréRer als 0 ist. Wie man 
leicht sieht, ist das Integral aber jedenfalls fiir alle z aus 0<2<d absolut beschrankt, und da 
x? wegen (11) beliebig groB sein kann, verschwindet (13). Dieselbe SchluBweise trifft auch auf 
Hal ) dx fiir den Summanden mit e—** zu. Fiigt man zu dem verbleibenden 


x 


d 
1 ‘ et** dy 
ei d 
204 [aoe [ve : 
ie, 0 
das verschwindende Integral 
d 
1 e-2 xd e(*) ev t 
eee cus e—*5 d 
2900 eer p(s) é 
E,, 0 


hinzu, so erhalt man fiir den ersten Summanden in (9) 


d 
1 big: 
gay | ere [e o(é) dé, 
- 0 


und das ist nach Verlagerung des Weges L,, in die imaginare Achse die Fouriersche Integraldar- 
stellung fiir y(z) im Intervall 0 bis d. Es bleibt also noch zu zeigen, da fiir t = 0 die drei letzten 
Summanden in (9) verschwinden. Zunichst sieht man sofort, daB fiir 0<z< d der dritte Summand 
in (9) fir t= 0 aus ganz ahnlichen Griinden wie eben verschwindet, wenn man die Voraus- 


d 
setzung (8) beriicksichtigt. weimalige partielle Integration von sf enaet FCG) dé im zweiten 
0 


Summanden von (9) liefert 
d d : d eT d é &y 
[ett f(E) de = or 4 | FE) dey weet f den hf (6) de, femne |e, Nhe) ce ne 
0 0 0 0 0 0 0 


Setzt man dies in den zweiten Summanden von (9) fiir t = 0 ein, so liefern fiir 0< z < d die ersten 
beiden Summanden der rechten Seite von (14) wieder verschwindende Bestandteile, wahrend der 


dritte Summand nach derselben SchluBweise wie oben bei dem q-Integral die Fouriersche 


Integraldarstellung von = | d&, [ f@) dé, ergibt, die sich gegen den vierten Summanden in (9) 
0 0 
weghebt. Somit ist bewiesen, da (9) die Anfangsbedingung (2) erfiillt. 
Um die Giltigkeit der Randbedingungen (6), (7) nachzuweisen, bedenke man, daB gilt 
dE, 
| Ex(2) — my cal =0, (15) 


a 


Bei z = 0 ergibt sich mithin aus (9) 


oo 


ad i ie ety dy ‘ Shope he a? [ e? ty dee? *4 (x) 
HO, t) — xq a (0, t) = nf em | Do(T) e dr + =| 24 oy) 1 


4 


co co 


ae a? 42 t d: —xd = 
x [ ae ev att dr ni xdaxe qd % Xo) [ ear dr 
0 


ee 2 ex) )—1)(1+ x4 x) 
0 


a2 eb dice? *4 6 (x) er + 2 *) P ea 
+ wi (2 *4 e(x) —1) (1—xq x) [0 a 
Tip 0 


Hier heben sich wegen (12) die beiden letzten Summanden auf der rechten Seite gegeneinander 
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weg, und es bleibt 
a0 2 P 
H(0, t) — x S° (0,1) = | ent tng doe | Oj(t)e— 8 dr. 
L, 0 


iad 


Dies aber liefert nach (11): 


a0 1 
HO, 8) — 95 0.1) =a; | ody | Hye) ede = LY (17) 
Ly 0 
nach bekannten Satzen iiber die Laplace-Transformation!. Also wird die Randbedingung (6) 


durch (9) erfiillt. 
Mit Hilfe von (16) ergibt sich ferner aus (9), daB dort bei der Bildung von #(d, t) + x4 . (d, t) 
z 
der erste Summand verschwindet, da sich e—?*¢ ¢(x) —1 im Zahler und Nenner des Integranden 
weghebt und der Integrationsweg L,, fiir t += 0 in die imaginare Achse verlagert werden kann, also 


L_,, die Integration iiber die imaginare Achse in entgegengesetzter Richtung bedeutet. Vom zwei- 


ten Summanden in (9) bleibt nur das Residuum bei x = 0 iibrig. Dieses ist nach (14) und (16) 
d 


d & 
i 1 
Ee f(G) dy +4 | dé, | f(&) d&, hebt sich also bei der Bildung von 6(d,1)+ xu 9° (d.t) 
0 0 0 

gegen den letzten Summanden von (9) fort, soda schlieBlich nur noch bleibt : 


co 


ao a2 etn and (Gxt a 6 e905) | 
O(d, t (0) d —xt att | 
(d, t) + xa=- (d, t) | ar oo ST | eure t 
6 


pyar bp pn 


fo) 


a? : et be du (1 + % a) ; J 

— Oy(t)ie—* Lady. 

| em ot) 1) + % 4). a(rye a 
0 


Lyt+l_, 


Hier verschwindet nun der erste Summand, wie man wieder durch Verwandlung des | in [ 


“ —H 


zeigt, und der zweite Summand ist aus denselben Griinden, wie sie oben beim Nachweis der Er- 
fillung der Randbedingungen bei z= 0 dargelegt wurden, identisch mit #,(t), womit dann ge- 
zeigt ist, daB (9) auch die Randbedingung (7) erfiillt. 

Nachdem so (9) als geschlossene Lésung des allgemeinen inhomogenen Warmeleitungspro- 
blems erwiesen ist, handelt es sich jetzt darum, die eingangs erwahnte Endbedingung — @(z, t) 
soll zu einer vorgeschriebenen Zeit t = t, einen vorgeschriebenen Temperaturverlauf 


(2, t.) = g(2) (18) 
annehmen — mit Hilfe von (9) zu formulieren. Damit (18) lésbar ist, wird von g(z) gefordert, 
daB es zweimal stetig differenzierbar ist und die folgenden Bedingungen erfiillt, die aus den 
Randbedingungen (6), (7) hervorgehen: 


Bult.) — (0) = —r0 (55) (19) 


Oults) — ald) = 4 (SE) (20) 


Zur Aufstellung der Endbedingung soll (9) umgeformt werden in der Absicht, eine Gleichung zu 
erhalten, welche g(z), f(z), Oo(t), Ou(t) in ihrer gegenseitigen Abhangigkeit leicht auswertbar enthalt. 

In (9) tritt in den ersten drei Summanden der gleiche Nenner e—?*4e(x) —1 auf. An seinen 
Nullstellen x,, welche also der Gleichung 


ore (1 —xo x) (1—*%q %) = My == CO's 9 
ees en, wes ey 


geniigen, liegen jedenfalls Singularitaten der Integranden. x = 0 ist Nullstelle, die ibrigen Null- 


1 Courant-Hilbert, a. a. O,. S. 202ff. 
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F ‘ : - Sad 
stellen miissen verschwindende Realteile haben, denn sonst kénnte der Betrag von e *n E(Hn) 
nicht gleich 1 sein, d. h. man kann setzen 

H. == 10, Mit Wy reel UR N., =a Ono (22) 


damitc auch 745, ,, eine Nullstelle ist. n durchlauft alle ganzen Zahlen von 1 bis+- «, da, 
=z 
wie man nach leichter Umformung von e oe &(%n)= 1 erkennt, o, der transzendenten Gleichung 


genugt 


Sin (Qo or a) (23) 


tz On d = 


welche abzaihlbar unendlich viele beliebig groBe Lésungen besitzt. Da fiir die folgenden allgemei- 
nen Betrachtungen die Kenntnis obiger Eigenschaften der ~, geniigt, braucht hier auf die Még- 
lichkeit der graphischen Ermittlung der «, aus (23) im konkreten Falle nur hingewiesen zu 
werden}. 

Wegen der Lage der x, auf der imaginiren Achse lassen sich der erste und zweite Summand in 
(9) fiir t4 0 in Summen itber die Residuen der Integranden an den Nullstellen des Nenners ver- 
wandeln, gekennzeichnet durch das Symbol 5’ f. Bezeichnet man also die vier Summanden in 
(9) der Reihe nach mit I, II, ITI, IV, so wird 


242 2 ) 
Pa be * du E,(z) pees: e® 4 dye Ex(z) 
I+ ll=; end e(%) ® fre dt rece Dag 2(¢—2*4 (4) )—1) | 
(24) 
x j fer | 
0 J 
Wegen 
d —2 xd —— 2(xo = %a) ¢! — %o X%q Hy Ee : 
FAC (x) =e Wye oe — 2, d — (ee 5a xa) (l er x2) = N(#n) (25) 
mit der Abkiirzung E,, (2) = E,(z) erhalt man dann aus (24) 
x—=+0o Ont 6) d ) 
e in Z ele 
Teor = wos | o(é) ef dé 
n=-+co a? “i t 4 d (26) 
a ve” En(z) eh a" de, Bu (2) a4 
Tat coy Re 2 N (xn » ff dé aie eae l [ s@¢ Fdé. 
0 =O) 0 J 


Die Entwicklung des Integranden des dritten Summanden nach Potenzen von x liefert unter Weg- 
lassung uninteressanter Glieder 


ee 2 : 2 . ? : 
[f04|-3 3s arm waa), bla le, (Oe, 
10 0 


0 


athe 
a? x N(0) 


Hierin bedeuten nach (10) und (21) 


rae 2 (z + xp) 29) 
ae = 4 (2+); (28) 
(as), = N20) ; (29) 


also liefert der dritte Summand in (26) 


Z+ HX d é, 
a [s, ino 5) dé peices J fE2) i, : (30) 


* Vel. dazu die Lésung von(23) fiir x)=x, bei M. Planck, Theorie der Warme, S.135ff., Leipzig 1932, 
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Die ersten beiden Summanden in (26) lassen sich leicht in Summen iiber n — 1 bis n — + co 
umformen, wenn man bedenkt, daB 


Ea(s) SE Bi) (31) 


1— x9 *n 


ist und die E,,(z)/V N(xn) (n = 1,2, ...) ein vollstandiges, normiertes Orthogonalsystem bilden, 
in dem mit 


d d 
@ ‘ _n(& E_n(é) 
B= {oS dg, Ef = | f(g) © ag 32 
4 ) Nn) J NG) oo 
die n-ten Komponenten von (2), f(z) bezeichnet werden mégen!. Dann wird zusammen mit (30): 
; on: 
SEO al we Ens) 
[+11 = e " — E+ - —:-— Bf 

ma /N(%n) ee Co) (33) 


a+ x g fs 
+ 2G) (oe (46) a8 + | as 716) a8). 


Da sich nun oben beim Beweis der Giiltigkeit der Randbedingungen gezeigt hat, daB If + 1V 
bei Anwendung der Operatoren 1—x, 0/dz bzw. 1+ x, 0/dz an den Rindern z= 0 und z=d 
verschwindet, dies aber auch wegen (15) und (16) fiir E,,(z) zutrifft, so ist wegen des Entwicklungs- 
satzes II+ IV nach E,(2)/N(%n) entwickelbar, d. h. es sind in (33) noch der letzte Summand 
und IV in diesen Eigenfunktionen auszudriicken. Eine langere, aber einfache Rechnung, die hier 
iibergangen sei, lehrt dann, daf 


oO 


ai Saxe t En(z) = 1 1 En(z) f ax? t 
Anny Absa? TPR ! Bee Efe ae (34) 


Bleibt noch der Summand III in (9). Man sieht sofort, daB es die Integrale f sind, welche einer 
0 
Verlagerung des Integrationsweges in die imaginare Achse im Wege stehen. Teilt man aber auf: 
oo t co 
| = a, + | , so kann man zeigen, dafs der vom zweiten Integral herriihrende Bestandteil 
0 0 t 


fiir alle z aus einem abgeschlossenen Intervall in 0< z< d verschwindet und also nur der vom ersten 
Summanden herriithrende Bestandteil iibrig bleibt, welcher die gewiinschte Verlagerung zulaBrt. 
Es wird also behauptet, daf 


a fs oe sey, de Same t oe *(4—-2) ‘ ges 
——= mae oF a |! yy, ma Card cal 35 
ni | e245) —1| 1 txox Dee alias eS rpg POI A33) 
LL, * t 
gilt fiir z ~0,d. Zum Beweise, der zunachst fiir den Wegbestandteil L, gefiihrt werde, fiihre 
man entsprechend (11) y = a?x? ein mit x = x, + ix; und also 
te =——Joa+yoer+ pf; wi=x%—a, (36) 
\2a | ie a 


d. h. Lae al fiir |B| 00. Dann ergibt sich fiir den Absolutbetrag des ij in (35) 
2a L 


+ co 
a | a Ae Ie e' dp eek ras 7 i 
my | 2m J le 2% 4 o(¢)—1| \ (LF Hq Hr)? + (A+ momi)® = (+ ea er)? + (CL + mame)? 
Li os 
Nun hat |e—* x4 e(%)—1| fiir « >a) >0 und f zwischen den Integralgrenzen sicher einen klein- 
sten positiven Wert B. Die Wurzeln in den Nennern sind sicher gréBer oder gleich 1, wobei das 


ort fa (fe dr of [84(x) 8" *) 
t t 


PE, : . 
1 Die En(z) sind namlich Eigenfunktionen des Problems ne ~2H, mit den Randbedingungen 


Ey, — % = =0 bei z= 0 und E, ee —( bei z —d. Dies aber bedeutet nach (15) und (16), daB 
; t % 
x = xn die Eigenwerte sind. Daraus folgt auch die Vollstandigkeit der E,(z). 
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Gleichheitszeichen bei x) = x = 0 sogar fiir den ganzen Integrationsweg gilt. Es ist also 


a 
EL 
L,, 


Hier kann man nun bei Betrachtung des ersten Summanden fiir z + 0 nach (36) setzen 


++ 00 fo) co 
= = | dp oll | 9o(z)| onl!) dr + ae | 4(z)| elt) i] ; (37) 
t 


+00 % 


— 2 Va+ Vorth z 
dpe“? =2]e¢**dB=2] e '*° dp . 
0 0 


ae 


Mit Hilfe der Substitution t = Va 4- i/o? + B? gelangt man dann fir hinreichend groBe « zu der 
Abschatzung 


+90 2 


V2a 


yo 


px: SS € eee 
dB ere Le at (o2 +2 2a-0+2a) do. 


== 


Man sieht also, da fiir z+-0 das linke Integral verschwindet, da auch die rechte Seite fiir 
o-oo zu Null wird. Das Integral | (9,(z)|e—*—9 dr in (37) bleibt jedenfalls wegen der Voraus- 
t 


setzungen (8) unter einer von « unabhangigen Schranke, womit dann das Verschwinden des ersten 
Summanden in (37) nachgewiesen ist. Dasselbe zeigt man auch fiir z==d von dem zweiten Sum- 
manden mit genau denselben Mitteln wie eben. In III aus (9) 1aBt sich dann aber der Integra- 
tionsweg in einen solchen um die Nullstellen des Nenners auf der imaginaren Achse verlagern, 
sodaB sich zusammen mit (34) fir t> 0, 0< z<d aus (9) ergibt: 
Sekt ale) pe 2 S11 Fale) py (geht 4) 
O(z,t) = e BES Got eee pa a 
?) pa VN(en) hye Da nV N(en) ( 
$ t (38) 


—xp, d 


co a*® net 4 ; ss 
2 é Yn pe —a unt -_ e —x*y ar 
pet See Be) oe [6.0 ode — fe ws) 
n= 0 3 


Die Funktion g(z) der Endbedingung (18) mu8 nun (38) fiir t = t,, 0< z< d geniigen. Wenn die 
Entwicklung des g(z) nach den E,(z) gleichmafig in z konvergiert, was angenommen werden soll 
bzw. im Einzelfall nachgewiesen werden muf, so folgt aus der Vollstandigkeit des Systems der 
E,(z), daB diese Entwicklung jedenfalls in 0< z< d die Funktion g(z) auch wirklich eindeutig dar- 
stellt und (38) also in den Komponenten der E,,(z) geschrieben werden kann. Das liefert schlieB- 


lich die gewiinschte auswertbare Endbedingung (Ef sind die Komponenten von g(2)) 


ax t 


ED 7 o2 
Ea = e n 8 BY? n i *n'g__ 


eC 
t t 39 
2 a? tn or “n ‘g 1 : —x3 at Pee ; —x2 ar 
f /NGn) 1 Bare, i O4(T) e dt —— Tebygies i 04(T) e dt A 
0 0 


Sind ¢, f, Jo, Ja gegeben und g gesucht, so liefert (39) natiirlich nichts Neues. Es ist dann ®(z, t,) 
= g(z) nach (9), und die Bedingungen (19), (20) sind erfiillt. Schreibt man aber g(z) vor, so stellt 
(39) eine notwendige Bedingung dar, der dann 9g, f, 9), 04 gehorchen miissen. Im allgemeinen wird 
(39) nicht ausreichen, um diese Funktionen zu bestimmen. Sind aber drei von ihnen gegeben, so 
wird sich die vierte aus (39) bestimmen lassen. Wenn unter den drei gegebenen #, und #, sind, 
so ist die vierte eindeutig bestimmbar. Ihre Existenz ist auch dann allerdings noch an den Nach- 
weis der Konvergenz der Reihenentwicklung nach den Eigenfunktionen gekniipft. Ist aber Dy 


oder #4 gesucht, so zeigt (39), daf® im Einzelfall noch eine weitere Behandlung dieser Gleichung 
t 


notig wird. Diese hat sich dann darauf zu griinden, dah fo@ . 


2 2 
a*r 
"dr unter der eingangs gemach- 


0 
ten Annahme der Endlichkeit von (tr) mit «, iiber alle Grenzen wichst, es sei denn, daB Ht) 
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identisch verschwindet. Es wird so z. B. weiter unten gezeigt werden, daB im Falle der Warme- 
produktion f(z) durch spezielle quasistationare elektromagnetische Felder und unter der Foxr- 
derung einer konstanten Endtemperatur g(z) — 9° die Funktionen P, Bo, Oa sogar eindeutig durch 
(39) bestimmt sind. Auf weitere Kombinationsméglichkeiten, da® z. B. zwei der im Spiele stehen- 
den fiinf Funktionen gesucht werden, wahrend drei gegeben sind, usw., mag hier nur hingewiesen 
werden. Wegen der bisher noch vdllig allgemein gelassenen Warmeproduktion f(z) umschlieBt 
(39) natitrlich auch Falle, in denen diese nicht durch elektromagnetische Felder geleistet werden 
kann oder mindestens die neue Aufgabe entsteht, eine aus (39) gewonnene Warmeproduktion 
durch solche Felder zu approximieren. Es soll hierfiir noch ein einfaches, anschauliches Beispiel 
gegeben werden. Legt man, ausgehend von der Anfangstemperatur y(z) = 0 die beiden Begren- 
zungsflachen der Platte auf die konstante Temperatur 9°, so mu sich nach unendlich langer Zeit 
die Platte iiberall auf die Temperatur 9° erwairmen. Stellt man sich aber die Aufgabe, dies schon 
zu der beliebig vorgegebenen Zeit t = t, zu erreichen, so erhalt man aus (39) eindeutig nach etwas 
Rechnung (x, = xq = 0) die Warmeproduktion 


2 (2 n+1)? 2° a? ty i! 
4 a® x 9 : 2 
fe) = 3 Sent tein n+l. 4 afi (40) 
n=0 
und damit aus (9) 
w sin(2n+1)77/ (2n+1) ata t, j- (2 n-+1)? at a? 
A eS (4g) 
Bet) = 452 2n+1 Sp ia by Bee to : ay 
n=0 


Man sieht aus diesen Ausdriicken sofort, daB f(z) = 0 fiir t, = co und #(z, t,) = 0° wie verlangt 
wurde. Uberdies kann man auch zeigen, da stets 0(z, t)< Oz, t,). Inwieweit man das gefundene 
f(z) durch elektromagnetische Felder oder vielleicht als atomare Warmeproduktion realisieren 
kann, soll hier nicht untersucht werden. Jedenfalls ist (40) mit einem quasistationaren Wechsel- 
feld in einem Dielektrikum geringen Verlustfaktors unvertraglich, wie man aus dem Folgenden 
entnehmen wird. 
Es mu8 noch auf den in (21)ff. ausgeschlossenen singularen Fall eingegangen werden, daB 
Xo, %q einzeln oder beide unendlich groB sind. 
Trifft dies zu, also x) = x4 = ©, so liegt bei z = 0 und z = dnach (6) und (7) die adiabatische 
Randbedingung vor, und das Orthogonalsystem der E,(z)/J/N(x,) muB fiir x,=0 bei n=0 durch 
sy ceed 
yN(0) yd 
zu einem vollstandigen erganzt werden. Im ersten Summanden der rechten Seite von (34) ist 
dann die Summation von n — 0 bis n = © zuerstrecken. Fiir das Residuum von II in (9) an der 
Stelle x» = 0 erhalt man statt (30) den Ausdruck 


4. Resthy = [re dé oe ate 4 | 2 fas [as | IG) ay. (43) 
( 0 0 0 


a | 
0 


Die Summe von (43) und IV aus (9) geniigt nun, wie man sofort sieht, der Randbedingung 


= —0 bei z — 0 und z = d, ist also nach dem durch (42) erweiterten System der E,(2)/V N(xn) 
e 
entwickelbar: 


(42) 


d fo) if 

= _1 S12) Fs as 
COR Aleka mee . 
0 


mit N(x,) = 2d nach (25). . : 
Man erhilt also schlieBlich fiir x, = ~q — © statt (38), da IIT in (9) wegfallt, 


d d 
ule!) =z | 9) datz [Sl de+ 


0 
= tt Ene) pe, 1S PO ie 
~ 1B . e 5 
pee ange ot 2, oi Ne) 


n=1 


(45) 
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Wie man sieht, unterscheidet sich dieser Ausdruck fiir den beiderseitig adiabatischen Fall vom 
allgemeinen Fall in (38) nur durch das Hinzukommen eines in t linearen Summanden, der bei 
t—> co allein fiir die Temperaturverteilung in der Platte maBgebend ist. Kin mit t wachsendes Glied 
mu hier bei ununterbrochener Warmeproduktion ja auch anschaulich physikalisch erwartet 
werden. Nach (19) und (20) wird fir die Entdemperaturverteilung verlangt: 


Ce. huethes 20 Weandeale =ni0 (46) 


os 


und zu der Endbedingung (39) mit entfallenden -Gliedern tritt die weitere Bedingung 
d d d 
f a(z) dz = f (2) dz +1, [ f(z) dz. (47) 
0 0 0 


Ist nur %) = © aber x, endlich, so lehrt eine Betrachtung, die Schritt fiir Schritt der eben durch- 
gefithrten analog ist, daB (38) und (39) bestehen bleiben, wobei aber die (7) enthaltenden Glie- 
der wegfallen ebenso wie in (9). Dasselbe gilt auch im Falle x; = 00, aber x, endlich, nur daB hier 
die Glieder mit 8 (7) in (38). (39), (9) zu streichen sind. 


3. Die Warmeproduktion in einem Dielektrikum geringen Verlustfaktors durch quasistatio- 
naire Wechselfelder. Hier soll die Frage beantwortet werden, wie denn ein quasistationdres 
Wechselfeld aussehen mu, damit es in der dielektrisch angenommenen Platte die Warmeproduk- 
tion f(z) liefert, bzw. welche f(z) mit einem solchen Felde vertraglich sind. Die in dem Volumen- 
element dV der Platte in der Zeit dt produzierte Warme ist nach der Maxwellschen Theorie 
ganz allgemein oG?dVdt, wobei o die Leitfahigkeit und € die elektrische Feldstarke bedeuten. 
Mit diesem Ansatz bei einem leitfahigen Dielektrikum durchzukommen ist aber nur so 
moglich, dafS man € sinusférmig mit der Zeit t schwanken ]4Bt und dann im Falle isotroper 
Stoffe, um die es sich hier nur handeln soll, o ebenso wie die Dielektrizitatskonstante ¢ noch 
von der Frequenz abhangig annimmt. Die Dispersion von o und € ist experimentell er- 
mittelt zu denken. Durch zeitliche Mittelung erhalt man dann fiir die pro Sekunde und 
ccm erzeugte Warme o ©, wobei hier © nur noch von den Raumkoordinaten abhangt und 


sein Betrag gleich dem - -fachen des Scheitelwertes der Feldstarke ist. Nimmt man ferner den 


Verlustfaktor des dielektrischen Plattenmaterials nicht zu hoch an, so daB © = — grad @ gilt, 
man also wie mit einem statischen Felde in einem idealen Isolator operieren kann, so gelten die 


beiden Gleichungen 
ed a@ 0@\2 0®\2 
aa + 3 — °° (55) ei (=) pele ~~ 


wobei c eine positive Konstante ist. Es wird sich zeigen, daB (48) nur fiir bestimmte @(x. z) und 
f (2) méglich ist, welche es jetzt zu bestimmen gilt. 

Unter der Annahme eines differenzierbaren f(z), folgt durch Differentiation der zweiten Glei- 
chung (48) nach z und x 


DD, = OFD; 05 
aA (49) 
20, ©,,+2 0,0, =f. 
Zusammen mit der ersten Gleichung (48) folgt weiter aus (49) wegen ®,,— ®,, unter der Voraus- 
setzung @,=- 0 das Gleichungssystem: 
2 
2 ORD A Pm ed weg | 
D; (50) 
1: OF 1 @,, —0': | 
Da die Determinante proportional f(z) ist, dies aber + 0 sein soll, ist (50) lésbar zu 
z QD, & @, 
@., ul OG aa LS . (51) 


vA is 
Hier entnimmt man aus der ersten Gleichung 


®, = ¢,(x) /f. (52) 


2f 
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Zusammen mit (48) und @,, — @,, folgt aber dann 


Si Giles : 2 
af Ot ey(x)<ec wegen 3? = (c—c?)f+0, (53) 
Ve Cy 


und hieraus entnimmt man, da die linke und rechte Seite nur einer gemeinsamen Konstanten 


k/2 gleich sein kénnen: 
9S) Sai le Ck? (54) 


Die Potentialfunktion @ 1aBt sich herleiten, wenn man nach (53) integriert : 


é rs k 
(2 ~ 48 == «+k, = — are eos 2). (55) 
: 1 


Nach (52) und (53) erhalt man damit die Gleichungen 
®, =Yocos(k x +h). Ff, | 
Fie ool 
@,, = —f Je sin (ie x + a 2 | 


(56) 


also 


an ee ee A oe k 
D=k, f/f dé ee ye cos(4 x + hy) i, +k. (57) 


Man bestatigt umgekehrt sofort, da (57) die Gleichungen (48) tatsiichlich lést. Setzt man f (2) 
aus (54) in (57) ein, so erhalt man schlieBlich 


k 
2 oe k 
O=— hice cos(& x + hy) + hy. (58) 


Ist k — 0, so entnimmt man aus (54), da eine konstante Warmeproduktion vorliegt, fiir welche 
nach (56) ®, und ®, konstant sind, ® also ein homogenes Feld liefert: 


® = Yc cos k,- Yk, z — Ye sin ky ky x + hy. . (59) 
Die dielektrische Erwarmung durch eine innere Warmeproduktion, welche allein von z ab- 
hangt, ist also keinesfalls auf homogene Felder beschrankt mit f(z) = konst. Vielmehr ist dies 
nur — mit sin k, = 0 — der praktisch naheliegendste Fall, weil die Elektrodenflachen dann"mit 
den Begrenzungsflachen der Platte zusammenfallen. Die Anwendungsméglichkeiten von Elek- 
trodenflachen ® = konst. nach (58) und (59) sollen hier, wo es sich nur um das Grundsatzliche 
handelt, nicht diskutiert werden. Dagegen sollen noch die allgemeine Lésung (9) der Warmelei- 
tungsgleichung und die Endbedingung (39) fiir das gefundene f(z) ausgewertet werden. 


4, Der Temperaturverlauf und die Endbedingung bei f(z) =k, e**. Obwohl eben f(z) = k, e** 
nur fiir ein Dielektrikum geringen Verlustfaktors hergeleitet wurde, tritt eine solche Warme- 
produktion doch auch in anderen wichtigen Anwendungen auf. Man denke etwa an eine Stahl- 
platte, welche einseitig durch induzierte Wirbelstréme erwarmt und dann durch Abschrecken 
gehartet werden soll oder eine absorbierte ebene Welle. Zunachst findet man aus (9) die Antwort 
auf die Frage, wie die Temperatur sich mit der Zeit und dem Ort dndert, wenn die Anfangs- 
bedingung und die Randbedingungen gegeben sind. Man hat nach (33) Ef zu bilden und f(z) in 
IV aus (9) einzusetzen. Es ist unter Zuhilfenahme von (21) 


d 
1 ne We ee We ee 2 unky 1—xk 1+ qk en) ‘ 60 
ie ss VNGen) i e E_,,(2) dz (k face pe) VNGn) F — Ho tn nee Xa ne ( ) 


und also nach - und (9) 


atx t 

oo FEn(z) pe 2 ky, 1 En(z) (. via 1) i — xk 1-— Hq k (k—x,) d 
t) : o EB, @ 

B(z, t) =p ies "hs n=l N(%n) (kh? — x7) %n\1— xo tn 1-+%q %n 


Se 


~ (61) 


ME ey et | 3 F aa eo *(d—2) i 9, ( hate dy 
Tt ev et dr = UPN) Cmca 4 
eit e 2 *4 ox) ast] L+n9x / o(T) e 1+ xy % J 1 
x 
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Fitr %) = xa = © (adiabatische Randbedingungen) entfallt der dritte Summand, aber gema (45) 
hat man noch zu addieren 


d 
k 
a | pls) de + 74 (e4— 1) 
0 
sodaB in diesem Falle gilt 


arnt 
Sy da Tn(z) eae 2 ak 1 En (2) fs 1 1) olk—*n) “) 
Bail? Vay ne Noa” ae 1 N Gen) (HP — x2) 22 (1 


nm 


(62) 


ee 


d 
1 k 
+ [ ol) det py (es) 
0 


Es muB hier darauf verzichtet werden, die allgemeinen Gleichungen (61), (62) auf ee An- 
wendungen zu spezialisieren. Statt dessen soll etwas ausfihrlicher ein Beispiel fiir den oben er- 
wahnten Fall durchgerechnet werden, da8 @,(t), du(t) gesuchte Funktionen sind, wahrend g(z) 
sowie f(z) = k, c** vorgegeben sind und g(z) offen bleibt. Trivialerweise kann dabei der Fall 
adiabatischer Randbedingungen auf er Betracht bleiben. 

Die ihn betreffenden allgemeinen Fragestellungen beziiglich g(z) und @(z) erledigen sich mit 
Ef nach (60) gema8 der Endbedingung (39) und der Zusatzbedingung (47) sehr einfach. Ist 
(2) gegeben, so ist mit (62) schon alles beantwortet, ist aber g(z) gegeben mit erfiillter Bedingung 
(46), so findet man aus (39) E? und aus (47) E?, also 

ns Et ; 

Fee eg baa ee Bt 2 (63) 
Das aber ist wegen der Eigenschaften (22), (23) der x,, die auch fiir %) = %q = oo erhalten 
bleiben, nicht fiir beliebige E8 méglich, die auf konvergente Entwicklungen fihren, sondern nur 
fiir solche, welche entsprechend dieser Gleichung auch auf E? fiihren, die eine konvergente Ent- 
wicklung nach den Eigenfunktionen definieren. Z. B. erweist sich die Endbedingung g(z) = konst. 
als unméglich, welche zwar (46) befriedigt, aber wegen Ef nach (60)keine konvergenzerzeugen- __ 
den E? nach (63) ergibt, es sei denn Ef=0. | 

Es soll nun, nachdem solcherweise der Fall adiabatischer Randbedingungen erledigt ist, die 
folgende Aufgabe gelést werden: 

Gegeben: g(z) = 0° = konst., 
Gesucht: mégliche 9,(t), Pa(t), @(z) bei endlichen x), xq. (64) 


Bei gegebener Warmeproduktion f(z) sollen also mégliche Rand- und Anfangstemperaturen 
bestimmt werden derart, da zur Zeit t = t, eine konstante Temperaturverteilung herrscht. Man 
kann natiirlich nicht erwarten, daB gerade bei diesem speziellen Beispiel sich solche gesuchte 
Temperaturfunktionen ergeben, die eine leichte praktische Realisierung erméglichen. Das Bei- 
spiel ist aber deswegen recht instruktiv, weil es die Lésungsmethode demonstriert und iiberdies 
ein Fall ist, in dem 9)(t), Ja(t), p(z) eindeutig durch die Endbedingung und die Warmeproduktion 
bei festen x), xq bestimmt sind. 

Zur Auswertung von (39) braucht man neben Ef nach (60) noch Eg. Es ist nach (21) 


d 
Hs po Nee = 2 90 ond lt Ho %n 1 
| NG) J ppetele VNGén) %in (1 — xo xn) Tg en i: 


Dann lautet (39) 


2 
© er GM 


2M@e "8 Hind 1 — 6 ohn 
—— —_—l(e ae ey | 
/N(%n) #n(1 —xo | 1 + %4 %n 


—gq?2 2 
See 2 ky (e mere 


Sa eee Ya 


a® | N(xn) xn (k2 —2xn) (1 — 2p 2n) (1+ %q %n) (1— xk) ‘ (65) 


od, 


t 
2 a? tn Be eer e (1 — x %n) ts ryan 
| = a) n d 0 a Hv 
: VN(«n) (1 ats) %n) (| oft) 1 e 1 ap %q %n | Balt) 7 ; 
0 0 


XX. Band 1952. Kriiger: Zur Kombination thermischer und elektromagnetischer Felder usw. 245 


Um hier die mit wachsendem x, unbegrenzt wachsende linke Seite zum Verschwinden zu bringen, 
setze man 


Bolt) = P+ OM), — Balt) = O° + OM(z) . (66) 
Dann wird aus (65) 


0— E42 kyle“ = 1) (ok =I) mf __ (Lx %n) (1+ xq k) F—%n) 4 
a? / Nin) xn (k2 —xn) (1 —xp %n) (1+ 2% vn) (1 — x ky | 


/ 


2 O81 aoe eee ni) 
ad a Hq tn (67) 
V Nin) Xn (1 —= 6 Xn) 
2 a2 x f 2 x, 4d] ‘5 2 
ie: f Oe) ert nee Ue ee (1) AP ene , 
Nn) ol 0 (T) ase Oy (zt) e dt 
0 0 
Bedenkt man nun, da8 im zweiten Summanden 
Li ape t 1 
22 (k? — x2) bk (i — «2 r a 
gesetzt werden kann, so fiihrt die Substitution 
eee te 8): | a 


O'})(r) ait fl SF xq k) ekd | k,(1+ 2g k) ekd si ek? a*(t—t,) ae @2)(z) | 


a? k? a? he 


zu einer weiteren Vereinfachung von (67), bei der der dortige zweite mit wachsendem x, un- 
beschrankt wachsende Summand beseitigt ist. Man erhalt 


2 po (1 __ L-xo%n uit 


Res gase 1+ xg%n 2 un ky (% k — 1) Cuan 
ss VNGén) %n(L — x9 %n) | N(2en) (1 —%q %n) 2 2(k® — 23) 
| 2 %n ky en) “ (1+ xik)e ka: 
' VNGén) (1 + %q %n) a2 (hk? — x2) 
oy 2; k,(%o k— 1) Yin ‘es a ecg tn) qd 4- %aq k) on) | 
a? |/ N(%n) k2(k®? — xn) (1 — xo xn) (1 + %4%n) (1 — xo k) 
oder 
1 —— ein nt } 
D0 a eek ees Oe ae 
= a 1+ Wain | 
‘ Nn) 2n(1 — %9 %n) ' (69) 


SIG kel) yl laen 8) 1 LH Hn) (1 tHa k) he | 
a? VN(xn) k2(1 — x9 n) (k® — xn) (1 + x4 %n) (1 —%k) 


Hier und in der vorhergehenden Gleichung ist bereits das sich aus (67), (68) ergebende Glied 


i 
g 
2 a2 xn r 1— 2% %n —%nd —a8 x27 
—— O?)(r) — e azN)|e dt (70) 
N(%n) (1 — %o%n) {| ) Dire ra oe: i 
0 


weggelassen. Zur Begriindung bedenke man, dafi nach (21) gilt 


Er 
spe es —r_d 1—HXo%n Vi HG Hy, 
GO Noe 2G = 


114 %n 1+ %4%n 1 — 23 x2” 
der Integrand in (70) also reell ist und (70) selbst mit wachsendem x, iiber alle Grenzen wachst, 
wihrend die Glieder in (69) alle beschrankt bleiben. Soll also das Problem lésbar sein, so miissen 
Q(z) und O)(r) identisch verschwinden. 

18* 
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Man sieht in (69), daB wegen der Eigenschaften der x, nach (22), (23) E? auf eine konvergente 
Entwicklung fiihrt, also ein ¢(z) definiert, das weiter unten nach der Berechnung von (2, t) ge- 
mi (9) sofort erhalten werden kann. 

Nach (66) und (68) ergibt sich fiir die Randtemperaturen 


k—1 
Seer 


[1 ew k(i—t,)) : | 
(71) 
Hi) Po ieee %4 k) ) ean — gat H(t t)) | 


welche also fiir ¢ = t, identisch mit 0° werden. 


Geht man mit (71), (69) und f(z) = k,e** in (9) ein, so folgt nach langerer Rechnung, die aber 
nichts Neues mehr bietet, 
ue ad t—ty 

Kz, t) = 99 —— eh [1 — )]. (72) 
Man sieht sofort, daB (72) fiir t = t, die Endbedingung ((z, t,) = 9° erfiillt und fiir t<t, iiber- 
dies nirgends die Endtemperatur iiberschritten wird. Die Randbedingungen (6), (7) sind ebenfalls 
unter Beriicksichtigung von (71) unmittelbar abzulesen, wahrend sich fiir t — 0 die Anfangs- 
temperatur zu 


ws ek #(1 — ev tg) (73) 


a* k? 


H(z, 0) = (2) = 9 


ergibt. Aus ihr erkennt man iibrigens, da8B nur fiir k — 0, d. h. eine raéumlich konstante Warme- 
produktion, eine rdumlich konstante Anfangstemperatur méglich ist. Nach (71) sind in diesem 
Falle die Randtemperaturen 0,(t) = 04(t) = 0°-+ k,(t—t,). Es liegt nahe, diese Uberlegungen 
auf allgemeinere Felder und mehr Dimensionen auszudehnen. Doch sollten hier nur an den 
einfachsten Fallen die charakteristischen Fragestellungen gezeigt werden. 


5. Zusammenfassung. Das aus der Anwendung hochfrequenter Wechselstréme auf die Er- 
wirmung von Materialien wachsende allgemeine Problem lautet: Durch Kombination elektro- 
magnetischer und thermischer Felder ist zu einer vorgeschriebenen Zeit ein vorgeschriebener 
Temperaturverlauf herzustellen. Unter der Annahme des Newtonschen Abkiihlungsgesetzes zu 
beiden Seiten einer unendlich ausgedehnten ebenen Platte der Dicke d wird zunichst bei all- 
gemein gegebener, ebenfalls nur von einer Variablen abhangiger Warmeproduktion f(z), eine all- 
gemeine Lisung des inhomogenen Warmeleitungsproblems in solcher Form gegeben, da daraus 
die Endtemperaturbedingung abgeleitet werden kann. Als erstes Beispiel wird aus dieser die 
Warmeproduktion berechnet, die nétig ware, um schon zu vorgegebener endlicher Zeit die Platte 
auf eine vorgeschriebene konstante Temperatur 7 zu bringen, wenn die Anfangstemperatur Null 
ist und an den Randern die konstante Temperatur 7° herrscht. Die Frage nach dem Aussehen 
quasistationarer elektromagnetischer Felder in Dielektriken geringen Verlustfaktors, die mit 
f (2) vereinbar sind, fiihrt auf f (z)= k, e**(k,, k Konstante), eine auch fiir andere wichtige Anwen- 
dungen giiltige Form der Warmeproduktion. Fiir dieses f(z) wird der allgemeine Temperatur- 
verlauf angegeben und die Endbedingung fiir das Problem ausgewertet, Anfangs- und Randtem- 
peraturen so zu bestimmen, da®B sich zu einer vorgeschriebenen endlichen Zeit eine konstante 
Temperatur in der ganzen Platte einstellt. Es zeigt sich, daB bei diesem Beispiel die gesuchten 
Temperaturverlaufe eindeutig bestimmbar sind und nur bei homogenem Feld (k = 0) eine kon- 
stante Anfangstemperatur mdglich ist. 


(Eingegangen am 11. Januar 1952.) 
Anschrift des Verfassers: Dr. Manfred Kriiger, Dachau, Obb., Grébenzellerstr. 6. 
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Uber Biegewellen in Staben. 
Von H. Schirmer. 


1. Aufgabenstellung. Die vorliegende Arbeit! befaBt sich mit Biegewellen in Staben, ins- 
besondere solchen, die durch stoBartige Beanspruchungen ausgelést werden. Im ersten Teil 
(Ziffer 2 bis 5) wird der Mechanismus der Biegewellen durch Zerlegung in zwei Elementarwellen 
anschaulich gemacht. Der zweite Teil kniipft an eine Arbeit von W. Fliigge® an, der den Stab 
als eindimensionalen Wellenleiter behandelt, der der von S. Timoshenko? angegebenen Gleichung 


is E 2 ke 
EF Wy xxx — O ke (1 + a6) Wy xtt +o Weeet wrt! W,, = 9 (1) 


gehorcht. Darin bedeuten E den Elastizitatsmodul, G den Gleitmodul, 9 die Werkstoffdichte, 
k den Tragkeitsradius des konstanten Balkenquerschnitts, A eine dimensionslose, bei der Be- 
rechnung der Querkraft aus den Schubspannungen auftretende Konstante, w die transversale 
Auslenkung des Balkens, x die Langskoordinate, t die Zeit und die Indizes in bekannter Symbolik 
partielle Ableitungen. Wir schreiben (1) sogleich in die fiir das Folgende zweckmaBigere Form 


1 1 1 wo? 
Wy%%% — |, ee Wx xtt v ae We tit et W,, = 0 (2) 


um, die die Ausbreitungsgeschwindigkeiten 
[E [GA 
cj = | —. C¢, = ae 
malig } 0 


der Longitudinal- und Transversalwellen im Stabe enthalt. Die Bedeutung der Frequenz w, = c,/k 
wird noch hervortreten. 

Die Gleichungen (1) und (2) beriicksichtigen sowohl Querkraftverformung als auch Dreh- 
tragheit der Querschnitte. W. Fliigge hat gezeigt, da® dies unbedingt erforderlich ist, wenn man 
zu endlichen Ausbreitungsgeschwindigkeiten kommen will. Die charakteristischen Geschwindig- 
keiten sind -} c,;, -+ c,; die erste ist die Laufgeschwindigkeit eines Biegemomentensprunges, die 
zweite die einer Querkraftunstetigkeit. Beide Spriinge breiten sich mit senkrechter Stirn und 
konstanter Sprunghéhe aus. Diese Ergebnisse gewinnt W. Fliigge durch eine physikalische 
Betrachtung des Wellenkopfes; man kann sie aber natiirlich auch mit Hilfe der Charakteristiken- 
theorie erhalten (vgl. F. Pfeiffer‘). 

Da die Gleichungen (1) und (2) neben dem Hauptteil cin Glied zweiter Ordnung enthalten, 
ist die Ausbreitung der Biegewellen mit Dispersion verbunden; ein aufgebrachter Stof 
verformt sich also langs des Ausbreitungsweges. Diese Verformung wird im zweiten Teil 
dieser Arbeit ermittelt. Die entsprechende Aufgabe hat A. Kromm® fiir die Kreislochscheibe 
behandelt, deren Lochrand einen Radialspannungssto& crfahrt. Hier bleibt die Sprunghéhe 
wahrend der Ausbreitung jedoch nicht erhalten, sondern nimmt entsprechend der radialen 
Ausbreitung mit wachsendem Laufweg ab. 

In den genannten Arbeiten ergeben sich die senkrechten Wellenstirnen stets als Folge der 
Vereinfachungen, die in der eindimensionalen Behandlung der Stabe und zweidimensionalen 
Behandlung der Scheibe enthalten und nur solange berechtigt sind, wie die Wellenlange des 
Vorganges groB gegen die Querschnittsabmessungen des Stabes bzw. die Dicke der Scheibe ist. 
Darauf hat L. Cremer® hingewiesen. Jeder geniigend steile Sto enthalt aber Teilwellen beliebig 
kurzer Wellenlange. Die senkrechte Wellenstirn ist also nur eine Annaherung; in Wahrheit ist 


1 Gekiirzte Fassung einer am Lehrstuhl fiir Mechanik entstandenen und von der Fakultat fiir Maschinen- 
wesen der Technischen Hochschule Hannover genehmigten Dissertation 1949 (D 89). Referent: Prof. Dr.- 
Ing. O. Flachsbart; Korreferent: Prof. Dr.-Ing. habil. A. Pfliiger. 

2 W. Fliigge, Z. angew. Math. Mech. 22 (1942), S, 312. 

3S, Timoshenko, Phil. Mag. (6) 41 (1921), S. 744. 

4 F. Pfeiffer, Z. angew. Math. Mech. 25/27 (1947) S. 83. 

5 A, Kromm, Z. angew. Math. Mech. 28 (1948) 5. 104. 

6 L. Cremer, Z. angew. Math. Mech. 23 (1943) S. 291. 
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der Wellenkopf verschliffen. Fiir den Fall der Telegraphengleichung hat H. Pleijel! die Ver- 
schleifung der Stirn berechnet, die dadurch entsteht, da die Stromverteilung itber den Draht- 
querschnitt frequenzabhangig ist. 

2. Zerlegung der Biegewelle in zwei Elementarwellen. Da wir sowohl Querkraftdeformation 
als auch Drehtragheit der Stabelemente beriicksichtigen, hat jeder Querschnitt zwei Freiheits- 
grade, transversale Auslenkung w und Drehwinkel yp 
(Abb. 1). 

Fir Biegemoment M und Querkraft Q gelten bekannt- 
lich? die Gleichungen 


M=—EFF y,, Q=GAF(w,—y), (8a, b) 
die das vollstandige Elastizitatsgesetz des Stabes dar- 
stellen. Hinzu treten die dynamischen Zusammenhange 

Ona 20 Lewy, M,=Q—o0FFYy, (4a, b) 


fiir Transversal- und Drehbewegung. Aus den vier Glei- 
chungen (3) und (4) 1aBt sich die Gleichung (1) oder auch das ihr gleichwertige System 


Abb. 1. Balkenelement, 


1 (oy) 
Wax x a rt = Opp) =. W,, 
1 C] 


(5) 


Wa a Ces Yo | 
& 


herleiten, das sich anschaulich verstehen la®t. Die linke Seite der ersten Gleichung (5), gleich 
Nuil gesetzt, beschreibt die reine Drehwelle des Stabes, die sich ausbilden wiirde, wenn man die 
transversale Auslenkung durch eine geeignete Abstiitzung verhinderte. Das rechts stehende 
Glied stellt die Kopplung zwischen Dreh- und Querwelle dar. Entsprechendes gilt fiir die untere 
Gleichung, deren gleich Null gesetzte linke Seite die reine Transversalwelle beschreibt, die ent- 
steht, wenn man die Drehung aller Stabelemente unterdriickt. Dreh- und Querwelle beeinflussen 
sich gegenseitig iiber die Kopplungsglieder der rechten Seiten der Gleichungen (5). 


blatttederntcy ov 


0 C5-ty b LZ hyy WAL 
Gs Q 
Abb. 2. StoBerregung einer Querwelle: Abb, 3. Schwingerkette als Modell der Querwelle: 
a) bei x = 0 aufgebrachter StoB Q,(t), a) Aufbau der Kette, b) Einschwingvorgang der 
b) Querwelle nach Ablauf der Zeit f,. Kette bei StoBerregung des Anfangsgliedes. 


Daraus erkennt man, daB in einem Stabe, in dem Biegewellen einer bestimmten Frequenz 
angeregt werden, zwei verschiedene Teilwellen laufen, den beiden Freiheitsgraden der Quer- 
schnitte entsprechend. Diese Teilwellen sind jedoch nicht Quer- und Drebwelle, sondern wegen 
der Kopplung zwischen ihnen Wellenformen, in denen in jedem Stabelement die vier Gré®en 
M,Q,y,w gleichzeitig schwingen. Dennoch ist es niitzlich, zundichst Quer- und Drehwelle 
getrennt zu betrachten. Um keine Riicksicht auf Reflexionen nehmen zu miissen, nehmen wir 
stets einen einseitig unendlich langen Stab an, der im Anfangsquerschnitt x = 0 erregt wird. 


3. Eigenschaften der Querwelle. 
Die reine Querwelle 


if 
Wax — 3 Wis = 0 
s 


1 H. Pleijel, Teknisk Tidskrift (E) 48, (1918) S. 129, vgl. K. W. Wagner, Operat h ipa 
1940, S. 238.) — ® W. Fliigge, a. a. O. : cochh oyun 2 aoa ae 
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gehért wie die Seilwelle zum einfachsten Wellentyp. Phasen- und Gruppengeschwindigkeit 
fallen zusammen und sind frequenzunabhingig. Ein dem Anfangsquerschnitt x = 0 antes. 
driickter QuerkraftstoB Q(t) (Abb. 2a), lauft daher unverformt als Kopfwelle mit senkrechter 
Stirn den Stab entlang (Abb. 2b). 

Die Querwelle 1aBt sich durch die in Abb. 3a dargestellte oo lange Schwingerkette veran- 
schaulichen, die bekanntlich einen Tiefpa® mit der Grenzfrequenz Cee Vcu[m darstellt?. Da 
hohe Fre quenzen in der Kette nicht fortgeleitet werden, schwingen die Glieder der Kette im Gegen- 
satz zu den Stabelementen stetig ein? (vgl. Abb. 3b), entsprechend der Aussage der Unscharfe- 
relation, nach der fir jedes Ubertragungssystem das Produkt aus Bandbreite und Anstiegs- 
dauer At von der GréSenordnung 1 ist®. 

Von der Gleichung des yten Schwingers der Kette 


e c 
WwW, ae (Wy 41 are 2, w, + W,_1) = 0 


gelangen wir zu der Gleichung der Querwelle zuriick, wenn wir die Federkonstante c, durch die 
Elastizitat GA/Ax, die Masse m durch die Tragheit 9 F Ax des Balkenstiickes Ax ersetzen und 
den Grenziibergang 4x — 0 durchfiihren. Dabei riickt die Grenzfrequenz w, nach oo: aus dem 
TiefpaB wird ein AllpaB. 


4, Eigenschaften der Drehwelle. 
Die Gleichung 


w? 


1 
Pit - a (6) 
1 


Wx “3 


c? 
der reinen Drehwelle stellt den Sonderfall «a =0 der allgemeinen Telegraphengleichung 
C2 Uy, — UWj,;—XU,—fPu=O0 dar, deren Eigenschaften bekannt sind*. Neben den zweiten 


i] 


aly 
i 


ezeta) 


b Durchlalbereiah— 
7) 


<= @z—> 0 0, QO, 
Abb. 4. Zur Drehbewegung der Abb. 5. Drehschwingerkette als Modell der Drehwelle: 
Stabquerschnitte. a) Aufbau der Kette, b) DurchlaBbereich der Kette. 


Ableitungen tritt die Funktion y selbst auf; daran ist zu erkennen, daB jedes einzelne Stab- 
. element fiir sich schwingungsfahig ist, und zwar, wie Gleichung (6) zeigt, mit der Resonanz- 
frequenz W,) = ¢,/k. Diese Fahigkeit kommt offenbar durch das Zusammenwirken von Quer- . 
kraftdeformation und Drehtragheit zustande. 

Das in Abb. 4 dargestellte Stabelement, in welchem die Stabachse um den Winkel w, = yo 
gegen die Ruhelage geneigt sei, erfahrt namlich durch die Querkraft das Drehmoment Q dx = 
GAF (wy) — yp) dx [Gleichung (3b)], das dem Moment —g Fk? pdx der Massentragheit das 
Gleichgewicht halt. Das Momentengleichgewicht liefert die Schwingungsgleichung 

P+ OY = Yo 
mit der Resonanzfrequenz @). Die Querschnitte schwingen dabei um die Ruhelage yo, d. h. um 
die Stellung senkrecht zur Balkenachse. Benachbarte Querschnitte sind iiber das Biegemoment 
miteinander gekoppelt. Die Drehwelle wird also durch das Biegemoment vorwarts getrieben. 

Zur Veranschaulichung diene die in Abb. 5a dargestellte unendlich lange Kette, in der Dreh- 


schwinger der ungekoppelten Eigenfrequenz 9 = /D/O iiber Federpaare vom Direktions- 
moment D, miteinander verbunden sind. Diese Kette ist ein einfacher Bandpa, dessen Durch- 


1 R. Grammel, Ing.-Arch. 14 (1943/44) 5, 213. 
2 J. R. Carson, Proc. Inst. El. Eng. (1919) 5. 407. i 
3 K, Kiipfmiiller, Elektr. Nachr. Techn. 1 (1924) 8. 141, 5 (1928), 5. 18; K. Franz, Arch. elektr. Uber- 


tragung 5 (1951) S. 10. é : é 
4 7. B. G. Doetsch, Theorie und Anwendung der Laplace-Transformation, Berlin 1937, Kapitel 21. 
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laBbereich durch die Frequenzen 
@,=YVDIO, ow =a@f1+4D,/Do 
begrenzt wird! (Abb. 5b). Die Bewegung des y-ten Schwingers wird durch die Gleichung 


ns D D 

D +e w— G41 —2 4+ pA) = 0 
beschrieben. Ersetzen wir, um zum kontinuierlichen Stab zuriickzukehren, O und D,, durch die 
entsprechenden GréBen 9 Fk? Ax und EF k?/Ax eines Balkenstiickes der Lange Ax sowie 


/D,/O durch c,/k, so erhalten wir nach dem Grenziibergang 4x — 0 wieder die Gleichung (6) 

der reinen Drehwelle. Beim Grenziibergang riickt @) 
oy ey wieder nach oo; aus dem BandpaB wird ein Hochpab. 
Mit dem Ansatz einer laufenden Sinuswelle 


p(x, t) es Ww et o(t— xc) 


crhalten wir die Phasengeschwindigkeit 


C(O) 


yey 


sowie die Gruppengeschwindigkeit 


Abb. 7. Drehwelle bei der Resonanzfrequenz «,. 


0 crt, (Abb. 6a). Hierbei ist s «w/e die Wellenzahl. Beide 
ABELe: Bigeacuartude: Desh wclles mk remnene: Geschwindigkeiten sind nur fiir @ > Wg reell, entsprechend 
gang der Phasen- und Gruppengeschwindigkeit, der HochpaBeigenschaft des Leiters. An der unteren Grenze 
4) MC), a Pepa de eee: ce eae des DurchlaBbereiches ist ¢ unendlich groB, was unmittel- 

bar anschaulich ist: alle Stabelemente schwingen in ihrer 
ungekoppelten Resonanzfrequenz mit gleicher Amplitude im Gleichtakt (Abb.7). Zu diesem 
Zustand gehért eine unendlich groBe Phasengeschwindigkeit und zugleich eine verschwindende 
Gruppengeschwindigkeit, da das Biegemoment iiberall verschwindet, ein Energietransport also 
nicht erfolgen kann. Ubrigens mu nach einem von H. G. Baerwald? angegebenen Satz die 


Gruppengeschwindigkeit an allen Haufungsstellen der Eigenwerte, das sind hier die Grenzen — 


des DurchlaBbereiches, verschwinden. Im DurchlaBbereich nimmt c mit wachsender Frequenz ab 
und v zu, da die Bewegungen benachbarter Querschnitte zeitlich um so mehr gegeneinander 
phasenverschoben sind, je héher die Frequenz ist. Im Sperrbereich schlieBlich kann keine fort- 
laufende Welle bestehen. La®t man némlich den Querschnitt x = 0 nach der Gleichung 


19 (0:1) ieee eh t (w < Wp) 


schwingen, so ist die Bewegung benachbarter Querschnitte durch 
4p (&, 1) == 2p (0, te? * mit B= c1V (ao)? —1 


gegeben; die Querschnitte schwingen phasengleich und mit értlich abnehmender Amplitude. 
Da itberall v < cist, ist die Dispersion normal. 

Der Frequenzgang der Phasen- und Gruppengeschwindigkeit kommt, wie oben dargestellt, 
durch kontinuierlich im Ausbreitungsgebiet der Welle verteilte Resonatoren zustande, die auf die 
Welle zuriickwirken. Das gleiche Dispersionsverhalten ist in der Elektrodynamik bei Wellen in 
Medien mit freien Elektronen bekannt, z. B. bei Rundfunkwellen in der Heavisideschicht3. 


1 R. Grammel, a.a.O. 
2 H. G. Baerwald, Elektr. Nachr. Technik 8 (1931), S. 224. 
3 W.O. Schumann, Elektrische Wellen, Miinchen 1948, S. 89 und 91. 
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Auch die beim Durchgang von Licht durch Materie gemessenen Dispersionskurven lassen sich auf 
diese Weise erklaren!; die Abweichungen dieser Kurven von’ den in Abb. 6a dargestellten 
erklaren sich aus der Absorption des Lichtes sowie daraus, da® das Licht in der Materie nicht 
eine, sondern eine Vielzahl von Resonanzstellen vorfindet. ; 

Erteilt man dem Stabanfang x = 0 den in Abb. 6b dargestellten Biegemomentensto8, so ist 
das Biegemoment im Querschnitt x gegeben durch 
OL fur. ath ati es 

a ; e 
A seen Flo Ve? = (/a)*) 

“1 [P= @/ay 


(= 

M(x,t)) — 

ke ‘le t fitr t> x/c,. “ 
U 


a/c] 


Dieses Ergebnis, in welchem J, die Besselfunktion erster Ordnung ist, ]aBt sich aus der bekannten 


Lésung des Einschaltproblems der allgemeinen Telegraphenleitung herleiten?. In dimensions- 
losen Koordinaten 


e x 
6 = Wy, 1] =Wyt 
Cy 


nimmt (7) die Form 


[= 0 fiir 7) << 's ° | 
J i - /-o__, £2 | 

MEE, ” = Mate URES) fir y DE { iy 
oa | 


an. Die laufende Welle (7) ist in Abb. 6c dargestellt. Die Wellenstirn bleibt bei der Ausbreitung 
erhalten; da die hohen Frequenzen dispersionsfrei iibertragen werden, ist dies verstandlich. Die 
Dispersion wirkt sich nur auf den Wellenriicken aus; dieser beginnt im Wellenkopf mit der 
Anfangstangente 

ee eed ee 

Oxo eae Perey, 
die im Laufe der Ausbreitung 
also immer steiler wird. 

Die StoBwelle (7) kann in 

der Fourierdarstellung 


M(x, t) 


+ co 
hgeo S(x, w) 
20 o 


et dy 
— oo 


geschrieben werden; die Spek- 
tralfunktion ist von der Form 


(4,0). == a | M(x, £)] 


. “ne, -f 
S(xo)=e AV % S(awj=e 2? 


Abb, 8. Spektralfunktionen: 
— Wy eed V1—(o/o,)? a) der Drehwelle (Laufzeit Ty, = x/cy, Parameter w/a). 
4! 


— 6 b) der Querwelle (Laufzeit T) = x/¢c,). 


+) 


worin das Zeichen & die Laplace-Transformation kennzeichnet, und in Abb. 8a dargestellt*. 
Sperr- und DurchlaBbereich treten auch in dieser Darstellung deutlich hervor. Zum Vergleich 
zeigt Abb. 8b die Spektralfunktion der reinen Querwelle. 

Wenn wir, wie Abb. 9 andeutet, den Balken auf seiner ganzen Lange elastisch stiitzen, wobei 
die Elastizitat kontinuierlich verteilt zu denken ist, tritt auch fiir die Querwelle eine Resonanz- 
frequenz Q, auf, die der Frequenz m, der Drehwelle entspricht. Dann sind beide Elementar- 
wellen vom gleichen Typus, und Gleichung (2) nimmt die Form 


1 il 1 (OR oO. => 2? we 22 
Wy «xx |r =~ —5 |) Wx xtt De Aap Re 2 Wxx Se 2 2 w= 0 
cy cs ay Mee Cs Tacs Gees 


an, die bis auf das Glied mit w,, hinsichtlich Quer- und Drehwelle symmetrisch gebaut ist. 


1 R. W. Pohl, Optik, 2. Aufl., Berlin 1941, § 106. 
2 7. B. G. Doetsch, a. a. O. S. 374. — * G. Doetsch, a. a. O. S. 368. 
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5, Eigenschaften der Biegewelle. Wir gehen nun zur Biegewelle iitber, wobei wir uns zunachst 
auf eine reine Sinuswelle 


w(x, t) — W eilow(t—4/e)] (8) 
beschrinken. Gleichung (2) liefert mit dem Ansatz (8) fiir die Phasengeschwindigkeit ¢ den 


Ausdruck 
1 el lg! yi myo: (9) 
Ge 2 \ cj Dy Ne ks cp e2\@ 
Die Abhangigkeit der Phasengeschwindigkeit von der Frequenz ist in Abb. 10 dargestellt. Es 
bestehen zwei Zweige c, und c,; die beiden zu einer Frequenz gehérenden Wellenformen unter- 


scheiden sich also hinsichtlich der Laufgeschwindigkeit und damit der Wellenlange. Form 2 
hat nur fiir @ > w, Wellencharakter, da cy fitr @ < wo imaginar ist. 


Zu der Welle (8) gehéren nach den 4° 
Gleichungen (3) und (4) die GréBen 
x0) 
w, = — I— Ww, 
c - 
reine Drehwelle 
yp=—ivew ( 2 == 1 w, Glegewelle ohne Querkrativertormung uLrehtraighet 
s NS ee 
22/1 I 
M=EFRo aa se 
€ 4) ——— 
S/ a 
at 
° WMC ele 
Q=i1CAF—w. a reine guerwelle 
c 
s Cs a 
i @ 
Seregeeeeeee ete @ 
i Se 
if a 3 
Abb. 9. Elastisch gestiitzter Balken. Abb. 10. Frequenzgang der Phasengeschwindigkeit der Biegewelle. 


Diese Beziehungen zeigen die Phasenwinkel zwischen den Verschiebungs- und KraftgréBen auf. 
In Abb. 11 sind die beiden zu einer Frequenz w > w, gehérenden Bewegungsformen im Augen- 
blick 1, wm (t — x/c) =a, d. h. beim Nulldurchgang der transversalen Auslenkung w in Richtung 
wachsender w, fiir ein Stabelement dargestellt. 

Bei Form I (c, < ¢,) haben die Drehwinkel w, und, da y < w, ist, auch die Biegemomente 
Qdx und M, dx den gleichen Drehsinn. Bei gegebener Auslenkung sind diese beiden Momente 
daher einzeln relativ klein, so daB die Welle 
von relativ geringen Kraften vorwarts ge- 
trieben wird, was die geringe Laufgeschwin- 
digkeit erklart. Die Wellenform liegt zwi- 
schen der reinen Querwelle (y = 0) und 
der Biegewelle nach der elementaren Glei- 
chung 

EP Wane amo rp (p=) 5 
und zwar fiir hohe Frequenzen naher bei 
Abb. 11. Die beiden Formen J (links) und 2 (rechts) der Biegewelle fiir der ersten, fiir niedrige Frequenzen naher 
® > @); der gezeichnete Augenblick entspricht dem Durchgang des bei der zweiten (Abb. 10). Fiir sehr hohe 
Staben durch dic Ruheloge die singetragenen GroBea haben thrt> Fyequengen wird die Wellenform mit der 

reinen Querwelle identisch und kann dem- 
entsprechend hier nur eine Querkraft, jedoch kein Biegemoment iibertragen. 

Fir Form 2ist c, > ¢,, und w, und y sind nun um 180° phasenversetzt. Dies fihrt zu groBen 
Werten der Querkraft. Noch gréfer als das Moment Q dx ist jedoch das diesem entgegen wir- 
kende Moment M, dx. Es sind starke Spannungen vorhanden, die die groBe Laufgeschwindig- 
keit c, und damit eine groBe Wellenlange bewirken. Dementsprechend ist die Neigung w, der 
Balkenachse relativ gering. Die Form unterscheidet sich nur geringfiigig von der reinen Dreh- 
welle, fiir die w, = Oist (Abb. 10). Fir sehr hohe Frequenzen wird dieser Unterschied itberhaupt 
zu Null; die Wellenform vermag dann nur ein Biegemoment, nicht aber eine Querkraft zu 
itbertragen. 


j Q Ruhelage dep 
Balkenachse(r0=0) 
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mie von W. Pliigge gefundenen senkrechten Wellenstirnen, die bei der Ausbreitung eines 
Querkraft- und eines Biegemomentenstofes auftreten, werden im ersten Fall ausschlieBlich vou 
der Form 1, im zweiten nur von der Form 2 vorwirts getragen. 

Wir fiithren nun die Abkiirzungen n = ¢/c,, b = (n? —1)/(n? +1 
Frequenzmas 


) und das dimensionslose 


n2—1 @ 


z= 


ZN Wo 


ein. Fiir Stabe mit Rechteckprofil und der Poissonschen Zahl m— 10/3 ist n=2 und 


3 
ay (@/@,). Die Wellenzahlen werden 


ips ie ae (10) 
Die Hilfsfunktion 


_ iim Ely — 
h;,. (2) Sareea zg? bzyl + 7? (11) 
ist in Abb. 12 dargestellt. Fir die Gruppengeschwindigkeit v = dw/ds finden wir 
hy,» (2) 
eee 2 We epee 2 
val?) = Pe eG) te) 
2 
mit der weiterem Hilfsfunktion g,,(z) = z-{ b* ey 
Vi+2 


Den Frequenzgang der Gruppengeschwindigkeit zeigt Abb. 13. Der Ast v, hat fiir alle Frequenzen 
die Eigenschaft v, > c;, was anomale Dispersion bedeutet. Fiir die elementare Biegewelle ohne 


\ 4 h(z) / 5 


a 


aa 
Za 


a Giegewelle ohne Querkrativertormung 
und Drehiraghert 


Abb. 12. Die Aste h, und h, von h (2). Die Werte ober- Abb. 13. Frequenzgang der Gruppengeschwindigkeit 
halb der z-Achse sind positiv reell, die unterhalb der- der Biegewelle. 
selben negativ imaginar zu verstehen. Die beiden rest- 
lichen, zuh, undh, beziiglich der z-Achse spiegelbildlichen 
Aste —h, und —h, sind nicht gezeichnet. 


Querkraftdeformation und Drehtragheit ist diese Eigenschaft lange bekannt; schon Lord Ray- 
leigh hat gezeigt, da fiir diese Welle v = 2c ist 1. Dieses Ergebnis ist auch anschaulich: bei 
gegebener transversaler Auslenkung sind die Kriimmung der Balkenachse und damit das trei- 
bende Biegemoment um so griéfer, je kleiner die Wellenlange ist. Daher ist dce/di < 0 und, was 
gleichbedeutend ist, v > c. 

Ubrigens ist anomale Dispersion auch bei anderen Wellen, deren vorwiarts treibende Kraft 
einer Kriimmung proportional ist, wie z. B. den Kapillarwellen einer Flissigkeitsoberflache, 
bekannt?. Anschaulich ist auch, daB der Ast v, unterhalb der Gruppengeschwindigkeit der 
elementaren Biegewelle verlauft, da die Beriicksichtigung der Querkraftdeformation eine Ver- 
minderung der Balkensteifigkeit bedeutet. Der Ast v, zeigt dagegen wie die reine Drehwelle nor- 
male Dispersion. 


6. Ausbreitung eines Biegemomentenstofes. Wir behandeln nun die Ausbreitung des in Abb.6b 


dargestellten Biegemomentenstofes 
+ co 
e 
Moo) = aes = do , (13) 


— co 


1 Lord Rayleigh, Theorie des Schalles, Anhang zu § 273. 
2 A. Sommerfeld, Lehrbuch der theoretischen Physik, Band 2, § 25. 
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der im Anfangsquerschnitt « = 0 des einseitig unendlich langen und bei x = 0 gelenkig gestiitzten 
Stabes aufgebracht wird, also den gleichen Fall, den W. Prager fiir die Gleichung E k? wy, xx 
++ 0 w,, = 0 behandelt hat!. Der Stab sei im StoBaugenblick in Ruhe: 


w(a00) ==10.(x,.0)i= 0 stir alles. 

Fiir den Querschnitt « 0 bestehen die Stiitzbedingung w(0,t) = 0 und die Momenten- 
bedingung 

EF keg, (0.0) BE hw, (051) =) (14) 
Der Ansatz 

w(x, t) — A ei(wt—s x) 
fiihrt wieder auf die Wellenzahlen (10) und auf die Lésung 
w(x, t) = (Ae * + Bets + Cette + Dest *) er, 


In dieser wollen wir unter s, (@) eine fiir alle @ > 0 positiv reelle Gréfe und unter s, (w) eine 
fiir © > W, positiv reelle, fir 0 << w < wm, negativ imaginare GréBe verstehen (vgl. die beiden 
Aste von h(z) in Abb. 12). Dann stellt A e'(@*+s:*) eine in Richtung abnehmender x laufende 
Welle dar, die wegen des Fehlens einer Reflexionsstelle verschwinden mu. Das Gleiche gilt 
beziiglich C e(@*+%*) fir w > w,. Fir 0 < w < wy kann dieses Glied Ce? * e'%' (Bf = —is, 
positiv reell) geschrieben werden; es bedeutet also eine Schwingung, deren Amplitude mit 
wachsendem x ansteigt, Diese kann ebenfalls nicht bestehen. Somit ist A= C=0. Aus 
w (0, t) = 0 folgt auBerdem B= — D. Die verbleibende Teillésung 


W(X, 8) == Blea hae ee 25 8) or 
laBt sich der Randbedingung (14) dadurch anpassen, daB unendlich viele Teillésungen ttberlagert 


werden: 
W(X; Dt i B(@) (e—? 81 * — e—* 52%) eh Ot da. 
Hier werden s, und s, auch fiir negative w bendtigt. 
Die Funktion B(q@) 14Bt sich aus den Gleichungen (13) und (14) bestimmen. Die voll- 
standige Lésung lautet dann 


+ co 


w(x, t) = — ont EFE e it day. 


(sj — 83) 


Fir das Biegemoment 


entsteht der Ausdruck 


UV) eiot day (15) 
mit der Spektralfunktion S(«,@), die aus zwei den beiden Wellenformen entsprechenden An- 
teilen S, und S, besteht: 

S (x, 0) = Sy(x, 0) + Sq (x, «) 


Qe ih ew 7 in n? —1 @ 


iG ays ements 21 Wo 


Cane 8, % 


1 | lees : 
2 Pe yeaah WR enh 2 di fasts | n? — 1\2/a@\2 
Pe eee ae pele Sega 

Oe deel eee eRe 

Mit den oben eingefiihrten Abkiirzungen wird einfacher 
—ih,(z) «/k —ih,(z) x/k 
S (a2) 1S) (aa Sh eee (1 tla pig’ Saif (ee 
il ) o( ) 9 yi+ 2 | D If V1+ 2 ~ (16a) 

Die Spektralfunktion ist in Abb. 14 fiir x/k — 5 dargestellt. S, beginnt auf der reellen Achse im 


(16) 


1 W. Prager, Ing.-Arch. 4 (1933) S. 16. 
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Punkte (1/2, 0) und durchlauft mit wachsender Frequenz eine Spirale, die mit @—> oo in den 
Nullpunkt einmiindet, was physikalisch bedeutet, daB Wellenform I bei hohen Frequenzen kein 
Biegemoment tibertragen kann. S, beginnt ebenfalls im Punkte (1/2, 0), bleibt im Sperrbereich 
@® < Wo jedoch auf der reellen Achse. Bei w = @, geht S, in eine Spirale iiber, die sich mit wach- 
sender Frequenz asymptotisch dem Einheitskreis nahert und den Nullpunkt wie ¢ *° 71 
(Ty, = x/¢) umlauft. Dieses Verhalten ist bekanntlich kennzeichnend fiir Ubertragungssysteme, 
deren Ubergangsfunktion nach Ablauf der Laufzeit T;, mit der Sprunghdhe Kins einsetzt). 
Wegen S,(x, 0) + S,(«, 0) =1 ist ferner nach bekannter Grenzwertrelation? lim M(x, i) = M, 
fiir alle x, wie es sein muB. tro 

Da an der Bildung des Wellenkopfes nur Wellenform 2 beteiligt ist, haben wir dort weit- 
gehende Ubereinstimmung zwischen Biegewelle und reiner Drehwelle zu erwarten. In der Tat 
istim Wellenkopf fiir beide Arten Q = w = vw, = py =0, 
M = M,, ars (/EFR, Wt = C1 Wxe 

Kin Unterschied besteht dagegen bei Q,; es ist 
namlich Q, = M, w3/cj fiir die Drehwelle, Q, = M, w2/ 
(c? —c%) fiir die Biegewelle, wie sich aus den vier 
Gleichungen (3) und (4) finden lat. Die genannten 
GréBen sind im Wellenkopf von « ganz unabhiangig. 
Dagegen wichst die GréBe M, bei der Drehwelle mit 
x, wie oben im Anschlu8 an Gleichung (7a) festgestellt. | 
Um auch fiir die Biegewelle den Wert von M, im \ 
Wellenkopf zu finden, fiihren wir in Gleichung (15) ein 
neues Zeitmah t, = t — x/c; ein, so da im Wellen- 
kopf 4,0 ist, und entwickeln die damit entstehende 
Unterfunktion 


+4Smag. 


ad Krels 


—=—2Z=0 2=076\+Reell 
oO o 
(O-®) 


S(x, @) oe w x[eq 


Abb. 14. Spektralfunktionen S, und S, fiir x/k = 5 
oO 


und n = 2, 

nach Potenzen von 1/w, wobei wir alle Glieder von héherer als erster Ordnung in 1/q fortlassen, 
da wir nur eine fiir die Umgebung des Wellenkopfes giiltige Naherungslésung suchen. Fiir diese 
ist es aber ohne Nachteil, da® die Unterfunktionen sich im Unendlichen singular verhalt und daher 
ihre Potenzreihenentwicklung nach 1/@ nicht konvergieren kann*. Wir erhalten so die Beziehung 


+ co 
7 | 
peed SJR ER, 
M le er Os 
M(x, ty) —* = Sati eles 
A] fa 2 @ 


die Korrespondenz! 4,74 (dort s =1i@ setzen) liefert schlieBlich 


M(x, t)) = M, Jy (Baa) , (17) 


wobei J, die Besselfunktion der Ordnung Null bedeutet. Fiir den Wellenkopf erhalten wir somit 
lim M(x, t,) = M, und 
t},70 v. 0. 


oj M, 
y 


x, (18) 


Preven Saye 


ty, 0 v. 0. Ox cj a cs 
Die drtliche Abteilung des Biegemomentes im Wellenkopf ist bei der Biegewelle gegeniiber der 
Drehwelle also um denselben Faktor gréBer wie die ortliche Ableitung der Querkraft. Abbau 
des Biegemomentes und Aufbau der Querkraft hinter dem Wellenkopf verlaufen somit bei der 
Biegewelle rascher als bei der Drehwelle. Dies ist auch anschaulich, da ein Stabelement bei der 
Biegewelle unmittelbar nach Durchzug des Wellenkopfes transversal beschleunigt wird. Bei 
der Drehwelle ist die Querbewegung dagegen unterdrickt. 


1 W. Oppelt, Grundgesetze der Regelung, Wolfenbiittel 1947, Tafel 15, S. 49. 
2 G, Doetsch, Tabellen zur Laplace-Transformation, Berlin 1947, 5. 7. 

3 G, Doetsch, Tabellen, S. 8. 

4G. Doetsch, Tabellen, S. 111. 
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7. Charakteristikenverfahren. Die Eigenschaften des Kopfes der sich ausbreitenden Welle 
sind damit bekannt. Um auch den Wellenriicken zu kennen, muf das Integral (15) ausgewertet 
werden. Dies sté®t nun auf so bedeutende Schwierigkeiten, da% wir uns nach einem Naherungs- 
verfahren umsehen miissen. Hier bietet sich die von F. Pfeiffer! angegebene Methode dar, ein 
langs der Charakteristiken fortschreitendes Differenzenverfahren. Da dort der Lésungsgang 
auch fiir den vorliegenden Fall skizziert ist, konnen wir uns kurz fassen. 

Wir breiten in der (x, t)-Ebene das in Abb. 15 dargestellte, nach den Charakteristiken der 
Differentialgleichung (1) orientierte Gitternetz aus und wahlen die kleine Maschenweite k. Von 
den bekannten Punkten der Charakteristik x = c,t bzw. der Ordinatenachse x = 0 ausgehend, 


So; 


Sana 
k ék 3k 
Abb. 17. Biegemoment und Drehwinkel fiir verschiedene Augenblicke (zweite 
Abb. 15. Gitternetz Naherung). Die strichpunktierten Linien sind die Tangenten im Wellenkopf 
fiir das Charakteristikenverfahren. nach Gleichung (18). 


1 
Abb, 18. Vergleich verschiedener Wellenformen fir den Augenblick t= aa 
1 
Iu. 2: Biegewelle nach Charakteristikenverfahren erster und zweiter 
Naherung; 3: Drehwelle, Lésung nach Gleichung (7); 4: Biegewelle, Lésung 
nach Prager; 5: Wellenkopf-Tangente zu 1 und 2 nach Gleichung (18) 


Abb. 16. Das durch- 
gerechnete Gebiet der 
(x, t)-Ebene. 


2 


6: Wellenkopf-Tangente zu 3 (m= wb ay - . 

die dem Wellenkopf bzw. dem Stabanfang entsprechen, kénnen wir die GréBen w, w,, w,, 
WY. Wx: Ys, M, Qin weiteren Punkten ermitteln, indem wir langs der Charakteristiken fortschreiten. 
Die in Abb. 15 ausgezogenen Linien entsprechen den charakteristischen Geschwindigkeiten —{ c;, 
die gestrichelten + c,. In dieser Weise wurde das in Abb. 16 angedeutete Gebiet durchgerechnet 
und zur Fehlerkontrolle die Rechnung fiir das halbe Gebiet (unterhalb der strichpunktierten 
Linie) mit halber Maschenweite 3k wiederholt. Eine solche Kontrollrechnung ist um so mehr 
erforderlich, als die Fehler sich unginstig fortpflanzen und mit zunehmender Entfernung vom 
Anfangsquerschnitt « = 0 eine wachsende Unsicherheit des Ergebnisses bewirken. 

Der Zeitaufwand, den das Verfahren verlangt, ist nicht unerheblich. Unter Benutzung einer 
halbautomatischen Rechenmaschine benétigen zwei Rechner, die zweckmabig zusammenarbeiten, 
fiir jeden Punkt etwa 5 Minuten zuziiglich etwa der halben Zeit fiir Kontrollen. Gerechnet 
wurden etwa 300 Punkte in erster und 500 Punkte in zweiter Naherung. Das Ergebnis ist in 
Abb. 17 und 18 dargestellt. Abb. 17 zeigt die Gré®en M und yw iiber der Lingskoordinate x 
fiir verschiedene Zeitpunkte. Im Laufe der Ausbreitung nehmen die Verformung des Riickens 
der M-Welle und die Tangentenneigung im Wellenkopf zu; letztere verhalt sich, soweit die 


1F, Pfeiffer, a. a. O. 
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beschrankte Genauigkeit des Naherungsverfahrens diesen Schlu8 zulaBt, wie Gleichung (18) 
angibt. Die y-Kurven sind bis auf einen Faktor das Integral der M-Kurven. Auch diesen Zu- 
sammenhang gibt Abb. 17 befriedigend wieder. Die Wendepunkte und Minima der y-Kurven 
liegen an den Stellen, die die Momentenlinien ihnen zuweisen. Dagegen liegen die Maxima etwas 
zu weit links, entsprechend der starken Kriimmung in der Nahe des Wellenkopfes. Die erste 
Naherungsrechnung gab diesen Zusammenhang nur sehr unvollkommen wieder. 

In Abb. 18 sind fiir den Zeitpunkt t = 10k/c, die Ergebnisse der ersten und zweiten Naherung 
miteinander verglichen (Kurve I u. 2). Die erste Naherung gibt die Verformung des Wellen- 
riickens also zu ausgepragt wieder. Die genaue Kurve verlauft danach oberhalb der zweiten 
Naherungskurve; ihr Abstand von dieser kann zu etwa 1/4 des Abstandes der ersten und zweiten 
Naherungskurve angegeben werden. Au®erdem sind in Abb. 18 zum Vergleich die Biegemomen- 
tenkurve fiir die Drehwelle (Kurve 3) nach Gleichung (7) sowie fiir die Biegewelle ohne Quer- 
kraftverformung und Drehtragheit nach der von W. Prager! angegebenen Lésung aufgezeichnet 
(Kurve 4). Bei dieser letzten Welle ist die Laufgeschwindigkeit des Wellenkopfes unendlich groB, 


8. Zusammenfassung. Die vorliegende Arbeit behandelt die Biegewellen der Stabe unter 
Beriicksichtigung der Querkraftverformung und Drehtraigheit. Der Mechanismus dieser Wellen 
wird durch Zerlegung in zwei Elementarwellen anschaulich gemacht. Fiir sehr hohe Frequenzen 
wird das Biegemoment ausschlieBlich durch die eine dieser Wellenformen iibertragen, weshalb 
deren Eigenschaften fiir den Kopf der Welle bestimmend sind, in der sich ein BiegemomentenstoB 
ausbreitet. Ein solcher Ausbreitungsvorgang wird untersucht. 


(Eingegangen am 22, Januar 1952) 


Anschrift des Verfassers: Dr.-Ing. Hans Schirmer, Mannheim, Heckerstrafe 11. 


1 W. Prager, a. a. O. 
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Stabilitatsprobleme bei tordierten Stiben und Wellen. 
Von A. Trésch. 


1. Einleitung. Die vorliegende Arbeit befaBt sich mit der Knickung von Staben (Teil T) 
und kritischen Drehzahlen von Wellen (Teil II), die durch ein Moment auf Torsion beansprucht 
sind. Beim Knickproblem werden wir auSerdem eine axiale Druckkraft mitberiicksichtigen. 
Wir setzen voraus, daf% der Vektor des aiuReren Momentes bei der Deformation seine axiale 
Richtung beibehalt (vgl. Abb. 2 und [1] *), und gerade diese Voraussetzung ist fiir uns wesentlich, 
indem ein solches Moment im allgemeinen kein Potential besitzt ([2], S. 279). Durch die An- 
nahme einer nichtkonservativen Belastung wird das Vorgehen grundlegend beeinfluBbt. Von 
den drei Stabilitatskriterien ([3]) diirfen das statische (als Frage nach den nichttrivialen Gleich- 
gewichtslagen) und das energetische Kriterium (als Frage nach der nicht mehr positiv definiten 
potentiellen Energie) nicht mehr ohne weiteres verwendet werden; wir miissen auf die urspriing- 
liche Bedeutung des Stabilititsbegriffes zuriickgreifen und die Schwingungen der Stabe und 
Wellen untersuchen. Ein auBeres Moment fester Richtung ist bei der Knickung wohl kaum 
realisierbar und tritt auch bei den kritischen Drehzahlen nur in seltenen Fallen auf; die Annahme 
liegt jedoch den meisten bisherigen Arbeiten zugrunde und soll hier beibehalten werden (vgl. [4)]). 

Die Fragestellung lautet demnach: Bei welchen axialen Belastungen oder bei welchen Dreh- 
zahlen entfernen sich die Elemente des Stabes oder der Welle unzulassig weit von der Ruhelage, 
so da der Stab oder der Rotor gefahrdet ist. Die Zusammengehérigkeit des Knickproblems und 
der Frage nach den kritischen Drehzahlen ist offensichtlich, wenn wir bei konservativen Be- 
lastungen die tiefste Knicklast als die Belastung ansprechen, bei der erstmals eine nichttriviale 
Gleichgewichtslage méglich ist, wahrend eine kritische Drehzahl dann auftritt, wenn auch bei 
verschwindender Exzentrizitat der Scheibe sich eine stationdre Auslenkung aufrecht erhalten 
kann ([5], 5S. 782). Wir werden, wie erwahnt, nicht diesen statischen Standpunkt einnehmen, 
sondern die Schwingungsrechnungen durchfihren; hier geht die Zusammengehorigkeit hervor 
aus einem Vergleich von (2.2) etwa mit (8.2). 

Unsere Aufgabe stellt vor allem ein Elastizitatsproblem dar. Bei der Knickung liegt ein Stab 
vor, der durch das auBere Moment, die 4uBere Druckkraft und auSerdem durch die Tragheits- 
krafte seiner eigenen Elemente belastet ist. Bei den kritischen Drehzahlen dagegen — wir ver- 
nachlassigen die Wellenmasse — besteht die Belastung in 4uBerem Moment und Druckkraft und 
auBerdem in den Tragheitskraften der rotierenden Scheibe. Wenn wir die Scheibenebene als 
senkrecht zur Ruhelage der Welle voraussetzen, so sehen wir damit von der Kreiselwirkung ab 
(Teil II, A); wahlen wir sie senkrecht zur ausgelenkten Welle, so wird die Kreiselwirkung mit- 
beriicksichtigt (Teil II, B). Die Knickaufgabe stellt ein Eigenwertproblem dar, wahrend die 
kritischen Drehzahlen auf ein Anfangswertproblem fithren. 

Wir beschranken uns iibrigens durchwegs auf Stabe und Wellen mit einer einzigen Biege- 
steifigkeit, so da8 wir mit Vorteil in einem raumfesten Bezugsystem rechnen. 

Im Teil I untersuchen wir die Knickung von geraden Staben in vier verschiedenen Lagerungen 
(Abb. 1), fiir die wir aus der Schwingungsdifferentialgleichung die Eigenwertgleichungen auf- 
stellen (Abschnitte 2 und 3). In den Fallen II und III verhalten sich die Stabe schon bei be- 
liebig kleinem Moment instabil ([2], S. 288); dies 4uBert sich darin, da® die Stabelemente sich 
auf Spiralen aus der Ruhelage entfernen. Wir nehmen dabei das Moment als klein an, was sich 
dadurch rechtfertigen 148t, daB sich die wesentliche Aussage auf kleine Momente bezieht. Hine 
Druckkraft beriicksichtigen wir vorlaufig (Abschnitt 4) nicht. Im Fall I besitzt das Moment 
ein Potential, und wir erwarten das Knickmoment, das sich aus der statischen Rechnung ergibt. 
Uberhaupt lassen sich die Ergebnisse des statischen Kriteriums in den Rahmen unserer 


* Die eckig eingeklammerten Zahlen, wenn nétig mit hinzugefiigter Seitenzahl, beziehen sich auf das 
folgende Schrifttum: 

[1] A. G. Greenhill, Proc. Inst. Mechan. Engr. (1883) S. 182; 

[2] H. Ziegler, Z. angew. Math. Phys. 2 (1951) S. 265; 

[3] H. Ziegler, Z. angew. Math. Mech. 31 (1951) S. 265; 

[4] H. Ziegler, Z. angew. Math. Phys. 3 (1952) S. 96; 

[5] C. B. Biezeno u. R.Grammel, Technische Dynamik. Berlin 1939. 
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Schwingungsrechnungen einfiigen. In dem Sonderfall, der der statischen Behandlung entspricht, 
brauchen wir uns iibrigens nicht mehr auf kleine Momente zu beschranken. Wir erhalten fiir den 
Fall IV als Ergebnis das Greenhillsche Knickmoment, wobei hier — die Belastung ist nicht 
konservativ — gezeigt werden mu8, daB es sich um das kleinste kritische Moment handelt. 

Da eine Druckkraft ein Potential besitzt, wird im folgenden Abschnitt 5 nur an einem Bei- 
spiel das bekannte Ergebnis der statischen Methode bestiitigt. 

Fiir die Belastung mit Druck und Torsion gleichzeitig (Abschnitt 6) lieBen sich die Rech- 
nungen mit kleinem Torsionsmoment in dhnlicher Weise durchfiihren wie fiir Torsion allein; wir 
werden jedoch auf die Wiedergabe verzichten und nur die Ergebnisse erwahnen. Kinzig fiir den 
Fall I berechnen wir die kritische Belastung, wobei reine Torsions- und reine Druckbelastung 
im Resultat mitenthalten sind. In der Art einer Greenhillschen Knickformel ({2], S. 276) 14Bt 
sich das Ergebnis nur als Naherungsformel darstellen. 

Es kénnen verschiedene Griinde dafiir geltend gemacht werden, da das auffallige Verhalten 
der beiden instabilen Falle bis jetzt nicht beobachtet worden ist ({4], S. 98). Im Abschnitt 7 
wenden wir uns einer Erscheinung zu, die dabei eine Rolle spielen kann, indem wir eine der 
Geschwindigkeit proportionale Dampfungskraft einfithren, Bei fester Dimpfungszahl laBt sich 
dann fiir das Moment eine Grenze angeben, unterhalb welcher sich die Stabe IJ und III stabil 
verhalten. Wie sich dieses Grenzmoment bei einer zusitzlichen Druckbelastung verandert, ist 
lediglich mit einem Zahlenbeispiel dargelegt. 

Im Teil II, in dem kritische Drehzahlen untersucht werden, sehen wir vorerst von der Kreisel- 
wirkung ab (Teil II, A). Wir fiihren in bekannter Weise ([5], S. 796) das Problem der kritischen 
Drehzahlen zuriick auf die Frage, welche Auslenkung eine Welle unter einer zu ihr senkrechten 
Kraft erfahrt (Abschnitt 8). Es ergibt sich, daB nur in zwei symmetrischen Fallen (Abb. 5a 
und d) die einzige bekannte kritische Drehzahl auftritt ([5], S. 800), wahrend in den iibrigen 
untersuchten Fallen jede Drehzahl kritisch wird; dies rithrt daher, da die Auslenkung nicht in 
die Kraftrichtung fallt und die Drehzahl tberhaupt nicht in die Lésung eingeht. Es ist deshalb 
auch nicht verwunderlich, wenn eine Dampfung — wir wahlen sie wiederum der Geschwindigkeit 
proportional — die Welle stabil machen kann, aber auch dies unabhangig von der Drehzahl 
(Abschnitt 10). 

Man kann sich fragen, ob das instabile Verhalten nicht durch die Kreiselwirkung doch noch 
wesentlich beeinfluBt wird. Diese Aufgabe, der die folgenden Abschnitte (Teil II, B) gewidmet 
sind, fiihrt auf zwei simultane inhomogene Differentialgleichungen; sie entstehen aus dem 
Schwerpunktsatz und dem Drallsatz (in der Gestalt der Eulerschen Gleichungen) fiir die Scheibe 
und dem Zusammenhang zwischen Auslenkung und Belastung der Welle (Abschnitte 11 und 12). 
Fiir die fliegende und die beidseitig gelenkig gelagerte Welle (Abb. 4) werden die allgemeinen 
Ergebnisse weiter verfolgt und dabei festgestellt, daB die Kreiselwirkung zwar die Welle nicht 
stabilisieren kann, aber immerhin die Amplitude der Schwingungen fiir Scheiben mit gréBerem 
Tragheitsradius weniger rasch anwachst (Abb. 8, Gl. (16.3)). Wenn wir das Torsionsmoment 
vernachlassigen, so kommen wir zu den bekannten kritischen Drehzahlen mit Kreiselwirkung 
zuriick ([5], S. 803), wobei sich allerdings nur der Gleichlauf als kritisch erweist, jedoch nicht der 
Gegenlauf. 

Weiter untersuchen wir die Wirkung einer — wiederum der Geschwindigkeit proportionalen — 
Dampfung; dies fithrt hier jedoch nicht zum gewiimschten Erfolg, indem auch eine grofe 
Dampfungszahl die Welle nicht stabilisieren kann (Abschnitte 14, 15, 17). 

Die als zweites Beispiel untersuchte gelenkig gelagerte Welle zeigt grundsatzlich das gleiche 
Verhalten wie die fliegende Welle: die ohne Torsionsmoment auftretende kritische Drehzahl ist 
nicht mehr eigentlich kritisch, sondern nur noch gefihrlich; dagegen wird tiberhaupt jede Dreh- 

zahl kritisch. Die Kreiselwirkung verlangsamt das (exponentielle) Anwachsen der Schwingungen. 


I. Die Knickung von Stiben. 


2, Die Schwingungsdifferentialgleichungen. Mit dem kinetischen Stabilitatskriterium, d. h. 
durch die Betrachtung der Schwingungen um die Gleichgewichtslage, sollen hier die Stabe bei 
Belastung durch ein Torsionsmoment W und eine Druckkraft D in folgender Lagerung unter- 
sucht werden: 

Fall I , der beidseitig eingespannte Stab, 

Fall IT , der einseitig eingespannte. Stab, 

Fall III, der auf einer Seite gelenkig gelagerte, auf der andern Seite eingespannte Stab, 
Fall IV, der beidseitig gelenkig gelagerte Stab. 


We 
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Wahrend der Fall I ein konservatives System darstellt, handelt es sich bei den drei tbrigen 
Fallen um nichtkonservative Systeme. 

Die Stabe werden homogen vorausgesetzt mit der Masse wu je Langeneinheit; ferner besitzen 
sie auf der ganzen Linge und unabhangig von der Achsenrichtung im Querschnitt dieselbe Biege- 
steifigkeit « = EJ und dieselbe Torsionssteifigkeit GJ,. 

Die Dyname in einem beliebigen Schnitt z besitzt is folgenden Komponenten ((2], S. 281): 


ee dx oe 
ee insti, a 
= d3x d*v dx ane dy | d?x dy : 
so nar Fe) Wor ae a ie Ie ue dz? ey. ; 
und die Sins ances lauten 
0? x Ox ory 08x ory ory ‘ 
Sat WS TDS thea Ox Le ee gtk 
sie lassen sich mit r= x-+1y vereinfachen zu 
anY —iWr'’4+ Dr’+ur=0. (2.2) 
Die Randbedingungen sind nach (2.1) 
p= 0, r=0, 


ar’—iWr =0, ar’ —iWr’+ Dr=0, 


LLL W 
D 
Z 
Z YZ L y 
I It 
Abb. 1. Die Lagerung in den vier Abb. 2. Die Belastung 
Knickfallen. des Stabes. 


je nachdem am betreffenden Ende die Ausbiegung der elastischen Linie, ihre Neigung, die Mo- 
mente M,, M, oder die Querkrafte K,, K, verschwinden. 
Mit dem Ansatz 


r(z, t) = R(z) T(t) 
und dem Separationsparameter ax‘ erhalten wir die beiden Differentialgleichungen 


ao RIV —iWR’’ + DR’ —axntR=0 


und 
uTt+axntT=0, 
die mit den Abkiirzungen Y= w , = yi aaundya ee kan 
o 
RIV —iwR” +bR” —xAR=0 (2.3) 
und 
T+#T=0 (2.4) 
iibergehen, wahrend wir als Randbedingungen 
R=, f eaeeee Ue 


R” —iwRk’=0, R'’ —iwR"’+bR =0 
erhalten. Aus (2.4) folgt als Lésung fiir den Zeitanteil 
T(t) = Aexp (tAt) + Bexp (—iAt). 


Wir ersehen daraus, da die Schwingung beschrankt bleibt, wenn J reell, aber verschieden von 
Null ist. Ausgedriickt in den x‘ heif®t das: die Schwingung ist genau dann beschrankt, wenn x4 


XX.Band 1952. Trésch: Stabilitatsprobleme bei tordierten Stiben und Wellen, 261 


auf der positiven reellen Achse liegt oder wenn, in Formeln, 
Imx4=0 und Rext>0 (2.5) 


gilt. Die beiden Beziehungen stellen zusammen also die Stabilitatsbedingungen dar. 


Der Ansatz 
R(z) == A exp(o 2) 
liefert uns als charakteristische Gleichung der Differentialgleichung (2.3) 
o*—iwo s+ bo®—x!=—0. (2.6) 


Die vier Wurzeln sind abhangig von x4 und der Belastung. Mit den Randbedingungen an den 
Stabenden erhalt man eine Kigenwertgleichung, die, wenn man sich die o aus der Gleichung (2.6) 
berechnet denkt, eine Beziehung nur noch zwischen der Belastung und x! darstellt. Fiir die 


Untersuchung der Stabilitat bendtigen wir also diese Eigenwertgleichungen; sie sollen im nachsten 
Abschnitt aufgestellt werden. 


3. Die Eigenwertgleichungen. Bei der Aufstellung der Eigenwertgleichungen der vier Knick- 
falle schlieBen wir von vornherein Doppelwurzeln aus — die Untersuchung der Mehrfachwurzeln 
hat auBer bei den unten besonders besprochenen Fallen auf kein erwahnenswertes Ergebnis ge- 
fiihrt —, so daB sich der Ortsanteil der Schwingung darstellen laBt als 


R(z) == Aj, exp(o; 2) . 


Fiir die Formulierung der Randbedingungen benétigen wir die Ausdriicke 
R'(z) = > Ax 6, exp(ox 2). R'(z) —iw R'(z) =  A,(ozf — i w ox) exp(ox 2) 
und 
R’'’(z) —iw R'(z) + b R(z) = DS) A, (of — i w of + b ox) exp (0% 2) 
wobei es unter Umstanden zweckmafiger ist, die beiden letzten Gleichungen mit Hilfe von (2.6) 


zu schreiben als 
R''(2) — iw R'(2) =o S Ay exp(o, 2) —b S Ax exp(or 2), 
k 
R" (2) — iw R"(2) + BR’ (2) = 24 S Ay = exp(or 2). 
Im Fall I lauten die Randbedingungen 


A, 0. >, A, 6, =”; > Ay exp (Go; !)i=— 0, > A, o, exp (o, 1) = 0. 


Damit nichttriviale Lésungen A, vorhanden sind, mu8 die Determinante des Gleichungssystems, 
die die Eigenwertgleichung darstellt, verschwinden: 


1 1 1 1 

O7 i) 03 O*4 

exp(o,!) exp(a,l)  exp(o3l) exp (o,!) 

0, exp (0,1) 0, exp (dy!) 63 exp (03!) 04 exp (04!) 

Sie ergibt ausgerechnet 

+ (6, — 0;) (64— 63) {exp [(o,+ 0) 1] + exp[(o3+ 04) I} 

== (0; Gj) (04 Gp) {exp [(o, + 95) 1] + exp[(o,+ 0%) UI} 

+ (64 — 0) (03 — 92) {exp [(o,+ 04) 1] + exp [(o2+ 93) I} =0 


oder 


Gin (o,—9,) 1/2 Sin(o,—o,) 1/2 ee Gin (o,—9,) 1/2 Sin(o,— 4) 1/2 
(5 cs 0;) (4 ae >) (05 = 0;) (o, co 63) 
Fiir den Knickfall IJ lauten die Randbedingungen 
eA == 0, > A,(of — i w o;) exp (ox. 1) = 0, 
1 
D> A, o, = 0, He 2: Ave exp (ti!) = 0 


19% 


262 Trésch: Stabilitatsprobleme bei tordierten Staben und Wellen. Ingenieur- Archiv 


und die Eigenwertgleichung 


O7 
% : 0% 
o? (6, —iw)exp(o,l)... 
exp (0, /) | 
Thre Ausrechnung — z. B. mit dem Laplaceschen Entwicklungssatz — lohnt sich in der all- 


gemeinen Form nicht. 


Im Fall IIT ist 
1 
2) Ap 0) 2 Ay oy he exp (G4) 10, 2 Ai, = exp (0x!) ==) 
k 


und 


07 exp (0; !) 
Im Fall IV lauten die Randbedingungen 
1 
SA=0, LAZ=0, SApexp(nl)=0, 2 Arzexp(al) = 9. 
k 


Als Eigenwertgleichung erhalt man 
1 
oF 
exp (0; /) 
oj exp (0; /) 
und daraus 
+ (05 — 0%) (of — 63) {exp [(03+ 04) 1] + exp[(o,+ 05) UT} | 
— (93 — 93) (of —o3) {exp [(0,+ 04) 2] + exp[(o,+ a5) U} 
+ (of — of) (63 — 03) {exp[(o.+ 95) 4] + exp[(o,+ 04) 1]} = 0 


oder 
Sin (o3 —9,) 1/2 Sin(oz—oy) 1/2 Sin(o, — 9,) l/2 Gin(o,— og) 1/2 
(03 — 03) (oti02) 5 1 ere (Go) (oj — 03) 


4, Knickung durch Torsion allein, Wenn auf den Stab einzig ein Torsionsmoment wirkt, so 
vereinfacht sich (2.6) zu 


ot —iwor—xi= 0. 
Dabei besitzt w gegeniiber x einen kleinen Wert ([2], S. 283), so daB wir fiir die Wurzeln der 


charakteristischen Gleichung eine Entwicklung nach w ansetzen kénnen, und zwar, wie man 
sich leicht iiberzeugt, als 


= iw 
o=ot Pte, 


wo o der w-freien Gleichung o4— x4 = 0 geniigt. 
Fir den Ortsanteil der Schwingung ergibt sich dann, immer unter Vernachlassigung héherer 
Glieder in w, 


4 : 4 
R(z)/e= > A), exp (0% 2) = exp ea) > Aj, exp (Gx 2) , 
be heal 


und je nach der Lagerung muB R (z) die entsprechenden vier Randbedingungen erfiillen. 
Der Knickfall I darf als konservatives System mit dem statischen Stabilitatskriterium be- 
handelt werden; er kommt erst in Abschnitt 6 zur Sprache. 
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Fiir den Knickfall IT erhalt man die Kigenwertgleichung ([2], S. 285) 
1 +- cos #1 Gp) x1 — 5, (cos “1 Sin xl -- sin x1 oj xl) = 0. (4.1) 
Thr kleinster Eigenwert betragt, ebenfalls nach w entwickelt, 
#1 = 1,875 —0,144iwl. (4.2) 
Wir schlieBen daraus, daB jedes noch so kleine Torsionsmoment mit fester Richtung eine kritische 


Last darstellt. 
Beim Knickfall III erhalten wir als Eigenwertgleichung in derselben Niherung 


sin x1 Wo) x1 — cos x1 Gin xl — e( — cos x1 Gp) x1) = 0, (4.3) 
die als tiefsten Eigenwert 
xl = 3,927 — 0,067iwl 
besitzt. Es erweist sich also auch der einseitig eingespannte und auf der andern Seite gelenkig 
gelagerte Stab als labil unter einem beliebig kleinen Torsionsmoment. 


Beim Knickfall IV lautet die Eigenwertgleichung 
8 x? sin xl Gin xl + w? |1 — cos x1 Coj x1 + * x1 (cos x1 Gin xl —sin x1 Coj x1)| = 0 


mit dem tiefsten Eigenwert 
%1 = (1 — 0,006 wl?) 


Der erste nichtverschwindende Teim ist hier quadratisch in wl. Der Stab schwingt also stabil, 
wird aber durch das Moment weicher und schwingt mit einer tieferen Frequenz. 

Um uns ein Bild von der Art der Instabilitat der Falle IT und III machen zu kénnen, sei die 
Schwingung noch etwas naher untersucht. Die Liésung der Differentialgleichung lautet fiir eine 
Eigenschwingung 


r(z, t) = [A exp(tAt) -+- B exp(— iANS Aj, exp (0% 2) . 


Demnach wird der Stab wahrend der Bewegung ahnlich verformt. Zerlegen wir / in Real- und 
Imaginarteil, 7 = 1, i/,, so 14Bt sich die Bewegung eines festen Punktes, z. B. des Endpunktes 
in Fall IJ, der Mitte in Fall III, beschreiben durch 
r (I, t) = A, exp[(tA,—/,) t] + B, exp[— (tA,—A,) ¢]. 
A, und B, sind durch Anfangslage und Anfangsgeschwindigkeit bestimmt. Die Bahn des Punktes, 
die fiir verschwindendes A, eine Ellipse darstellen wiirde, wird fiir kleine A, zu einer elliptischen 
Spirale. Die Ellipse wird dabei immer runder und gehi, sobald exp (—|/,|t) geniigend ab- 
geklungen ist, mit 
r (I, t) ~exp[(sgn Ay) Ay — tA) #] 
in eine Kreisspirale iiber. Wir fragen nun, auf welchen Wert die Amplitude angewachsen ist, 
wenn der Punkt in der Zeit t einen vollen Kreis durchlaufen hat. Es ist dann A, 7 = 2% und also 
pet = 200 A,//,. 
Damit erhalten wir im Knickfall IJ 


7014 
yt = 2 2( 1,875 jwl~po. 


Wo Q, den totalen Winkel der Verdrehung bedeute, die der Stab unter dem Torsionsmoment W 
erfahrt. Das Amplitudenverhaltnis nach einem vollen Umlauf wird damit 
An+ 1 
An 
d. h. die Amplitude wichst bei zulassigen Torsionsmomenten recht langsam an. 


Im Knickfall [JJ wird 


~expYyo~1+ o> 


An+1 14h: 


An 
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das Anwachsen erfolgt noch langsamer. Vergleicht man jedoch das Anwachsen der Amplitude 
in gleichen Zeiten, so zeigt sich, daB es fiir Fall IT und III ziemlich genau gleich rasch vor sich 
eht. 

Knicklasten, die sich aus dem statischen Kriterium fiir die Falle J und IV berechnen lassen, 
erweisen sich bei unserem Problem ebenfalls als kritisch, wenn wir die Stabschwingungen unter- 
suchen. Denn die statische Behandlung stellt nichts anderes dar als die Vernachlassigung der 
Tragheitskrafte, was wir durch w = 0 und damit auch x = 0 beschreiben kénnen; dieser Wert 
von x fiihrt auf eine linear mit der Zeit anwachsende ,,Schwingung‘‘. 

Im Fall IV erhalt man mit der Methode des letzten Abschnitts fiir x = 0 leicht die Beziehung 
wl =2an (n ganze Zahl), was nichts anderes als die Greenhilische Knickformel fiir D = 0 dar- 
stellt ({2], S. 276). Es erhebt sich nun jedoch die Frage, ob dieses Knickmoment fiir n = 1 auch 
das kleinste Knickmoment ist. Durch eine nahere Untersuchung von x4 in Abhangigkeit von w 1 
1aBt sich tatsachlich zeigen, daB fiir O< wl< 22 die GriBe x* niemals auBerhalb der positiven 
reellen Achse liegen kann. 

Das statische Kriterium liefert bei diesem nichtkonservativen System den richtigen Wert fiir 
das Knickmoment. Da dieses Ergebnis keine unerwartete Aussage darstellt und der Beweis der 
Behauptung ziemlich weitlaufig ist, soll er hier nicht wiedergegeben werden. 


5. Knickung durch eine Druckkraft. Eine Druckkraft mit fester Richtung im Raum besitzt 
ein Potential, so da das statische Kriterium hier zulassig ist. Die Ergebnisse der Schwingungs- 
rechnung sollen fiir den Fall IV nur als Beispiel angefiithrt werden. 

Bei einer Druckbelastung schreibt sich (2.6) als 


o#+boe?—xi=0, 


und ihre vier Wurzeln sind 


i b b2 
Cla ae oa \/-= +e +e. 


Fir den Fall IV erhalten wir aus (3.1) leicht 
Gin o,l- Sin o,1 = 0 


mit den Lésungen 
G,/=ina oder o,l=1nn (7 ganze Zahl). 
Sie fithren beide auf dieselbe Bedingung 
Gem? . wee 


i b 7 m 


und der Stab wird nach (2.5) labil, wenn 
Ton 
12 
ist. Als Grenze erhalten wir somit die Eulersche Knicklast 
D A 


iF Teas 


b> 


6. Knickung durch Druck und Torsion. Wenn wir zur Torsion die Druckbelastung mit- 
beriicksichtigen, so sind wir wiederum gezwungen, uns auf kleine Torsionsmomente zu be- 
schranken (|w/x| <1). Die Eigenwertgleichungen lassen sich dann immer noch in der Art von 
(4.1) oder (4.3) schreiben; die Schliisse, die man daraus zieht, seien hier aber nur angedeutet. 
Die Falle IJ und IIT, die schon beim Fehlen einer Druckkraft labil sind, bleiben labil, wenn wir 
die Stabe zusatzlich mit einer Druckkraft belasten und sogar auch dann, wenn diese Kraft eine 
beliebig groBe Zugkraft darstellt. Hier miissen wir also, um die Wirkung der Druckkraft zu 
erfassen, offenbar die Frage im Zusammenhang mit der Dampfung aufwerfen, was im nachsten 
Abschnitt geschehen soll. 

Beim Fall IV erhebt sich einzig die Frage, ob man ihn wie einen Knickfall mit konservativer 
Belastung behandeln darf. Da dies fiir den nichtkonservativen Teil der Belastung (das Mo- 
ment W) der Fall ist, diirfen wir wohl erwarten, daB die Druckkraft, die ein Potential besitzt 
daran nichts andert. Ein Beweis miiBte in dahnlicher Weise gefiihrt werden, wie wir es fiir ane 
Torsionsmoment allein angedeutet haben. 
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Obschon wir fiir x40 die Ergebnisse des statischen Stabilitatskriteriums erwarten miissen, 
wollen wir hier die Ergebnisse der Schwingungsrechnung fiir die Falle I und IV wiedergeben, 
und zwar deshalb, weil wir im Fall I zu einer Knickformel gelangen, die von derjenigen Greenhills 
abweicht ([1]; [4], $.97). Nicolai gibt das richtige Ergebnis an, berechnet jedoch fiir Druck- 
krafte keinen Punkt mit zusammengesetzter Belastung (Torsion und Druck gleichzeitig) der 
kritischen Grenzkurve (vgl. [5], S. 547). Diese kleine Liicke wollen wir noch schliefen, 

Die Wurzeln der Gleichung (2.6) sind jetzt 


O12 => (w® + uw? +45), OC, == 0,0. 
Fitr die Formulierung der Randbedingungen benétigen wir au®er 
R(z) = A, exp (0,2) + A, exp (a, 2) + As + Ay 
und seinen Ableitungen den Ausdruck 
RY — iw R = — [Aj 0,0, exp (0,2) + A, 0, 0, exp (022) + Ay(o, + Oy). 
Im Fall IV erhalten wir als Sakulargleichung einfach 


exp (0; !) = exp (03!) 
mit den Liésungen 


o,/=o0,1+2inn (7 ganze Zahl), 
was unmittelbar auf die Greenhillsche Knickformel 
2 D 2 
a gore ia 
fiithrt. 
Fiir den Knickfall I schreiben wir als Eigenwertgleichung 
exp (— 01/2) exp (— 01/2) 1 0 
o,lexp(—o,1/2) ol exp (— 01/2) 0 Us es 0 
exp (0; 1/2) exp (01/2) 1 1 
0, ! exp (9, 1/2) 021 exp (dz 1/2) 0 1 


aus der wir nach einigen Umformungen die Gleichung 
Gin(o,—o,) 1/2 Gino, W/2 Sino, 1/2 


(6,— 6,) V2 o, 1/2 a, 1/2 


oder 
sin (wo? + 46) 1/2 — sin(w + w? +45) 1/4 sin(w —|w? + 45) 1/4 
(/w? + 4b) 1/2 (w+ Vw? +464 (w —w® + 45) 1/4 
erhalten. Im Grenzfall fiir verschwindendes Moment ergibt sich daraus 
sin bl ( Vo nl 
Jo-1 Vo ,- 2! 


mit der Eulerschen Knicklast 
D 47? 
7 ee 
als tiefster Knicklast. Fir verschwindende Druckkraft wird die Higenwertgleichung zu 
tg (wl/2) = wi/2 
mit der tiefsten Wurzel 
Wi 
~— == 4,4934 . 
2 & 
Statt einen Grenziibergang fiir kleine 6 zu vollziehen, gehen wir hier besser direkt von der nun- 
mehr giiltigen Lésung 
R(z) = A, exp (i wz) -+ A,+ 43% + Ay? 


aus, 
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Der die beiden Punkte w = 0 und b = 0 verbindende Kurvenzweig wird am einfachsten aus 
einer Wertetabelle sin x/« erhalten. Daraus la6t sich weiter in einem Koordinatensystem mit 
den Achsen DI?/x und W21?/4 «2 die kritische Grenzkurve zeichnen ([5], S. 547). Dagegen 1aBt 
sich der Zusammenhang nicht in einfacher Form als Knickformel hinschreiben. Allenfalls konnte 
man die Knickformel 


= 47? 


3 |WIN DE 
a a T 


beniitzen, die fiir D >0 auf wenige Prozent genau ist. 


7. Die Dimpfung. Der Hauptgrund, weshalb das instabile Verhalten der Knickfalle IJ und II] 
unter Torsion nicht beobachtet worden ist, diirfte darin liegen, daB die Annahme eines Momentes 
mit fester Richtung dem wirklichen Kraftangriff nur in seltenen Fallen entspricht und man auf 
die das Moment erzeugenden Krafte zuriickgehen mu, um den in Wirklichkeit meist vorliegenden 
Verhaltnissen gerecht zu werden ([4]). Da wir hier jedoch die Voraussetzung eines Momentes 
mit fester Richtung beibehalten wollen, sei als stabilisierende Erscheinung die Wirkung einer 
Dampfung untersucht. Fiir die Dampfungskrafte wahlen wir den einfachsten Ansatz, namlich 
eine der Geschwindigkeit proportionale Dampfungskraft in jedem Stabelelement. 


War ohne Dampfungskraft die Belastung des Stabes je Langeneinheit durch die Tragheits- 
kraft 


07% ory 
Pa aaa Py=—H# Or 
gegeben (Abschnitt 2), so wird jetzt mit der Dampfung 
a 02x Ox Ory oy . 
a a ERR oD ee Em ne POT aT Te A er 


y hat dabei die Bedeutung einer Dampfungszahl. Die Gleichung (2.2) lautet damit 

antY —i{Wr'’ + Dr’ +urtyuyr=0 
und separiert sich wiederum in 

“a RIVY—i WR” + DR’ —axtR=0 
und 
lt 4h +uy iT +ax'T = 0 
oder 
T Ly T 17270. 
Die Lésung der zeitabhangigen Gleichung, die uns im Augenblick vor allem interessiert, lautet 
T(t) == A, exp [(— y + Vy? — 47?) t/2| + A, exp (ie y— Vy?— 422) 1/2]; 


dabei ist y reell und, weil wirklich eine Dampfung vorhanden sein soll, y > 0. In diesem Fall 
ist die Schwingung stabil, wenn die beiden Ungleichungen 


Re(— y + Vy? — 47?) <0 
erfiillt sind. 


Re(— y + Vy? — 422) = 0 
stellt offenbar die Grenze des stabilen Gebietes dar. Sie 1aBt sich auch schreiben als 
(Im. /?)?= y? Re 7? , 
und die Schwingung ist genau dann stabil, wenn 
(Im /?)? = y? Re 2? 
ist. Wir sehen daraus z. B., daf fiir Re? < 0 die Schwingung immer aufgeschaukelt wird. 


In den beiden unter kleinem Moment labilen Knickfallen besitzt das 4 einen gegeniiber dem 
Realteil A, kleinen Imaginarteil i/,. Fitr A =/,+ id,, wobei 22< 2 ist, vereinfacht sich die Sta- 
bilitatsbedingung zu 


ABZ. 
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Bei schwacher Aufschaukelung ist die Schwingung genau dann stabil, wenn 
2|ImAl=y (y > 0). 
Im Knickfall IT gilt, falls keine Druckkraft vorhanden ist, nach (4.2) 


xl = 1,875 —0,144iwl, 
und damit wird 


1y/a (969 __ 
eG tee rine ee 
Nach W aufgelést lautet dann die Stabilitatshedingung 


Wear <Wo= 1,91 Va u y/2 . 


Diese Stabilitatsbedingung mit demselben Faktor (innerhalb der angegebenen Genauigkeit) 
erhalten wir auch fiir den Knickfall IJI. Bei gegebener Dampfungszahl] schwingt also der Stab 
bei Momenten, die kleiner als W, sind, stabil, bei gréBeren labil. 

Obschon wir iiber die Gréfe der Dampfungskraft keine Angaben machen kénnen, schlieBen 
wir aus der Betrachtung der Schwingungen in Abschnitt 4, da8 offenbar auch bei gréBtem zu- 
lassigem Moment eine recht geringe Dampfung schon stabilisiert. 

Beriicksichtigt man nun noch den Einflu8 der Druckkraft, so 1aBt sich bei fester Dampfungs- 
zah] das Grenzmoment W, fiir jede Druck- oder Zugkraft berechnen. So erhéht sich z. B. im 
Fall II das W, auf annihernd den doppelten Wert, wenn wir eine Zugkraft vom Betrage der 
Knicklast (4 DI? = az?) wirken lassen. Driicken wir andrerseits den Stab mit einem Viertel der 
kritischen Last (16 Dl? = 47”), so sinkt W, um einen Sechstel. 

Mit diesen zwei Zahlenbeispielen schlieBen wir die Untersuchungen der Knickung von geraden 
Staében unter Druck und Torsion ab. 


Il. Kritische Drehzahlen. 


A. Kritische Drehzahlen ohne Beriicksichtigung der Kreiselwirkung. 


8. Der Zusammenhang der kritischen Drehzahlen mit der statischen Durchbiegung der Welle. 
Wir untersuchen zuerst die kritischen Drehzahlen einer kreiszylindrischen Welle, die eine diinne 
’ Scheibe tragt. Die Scheibenmasse m sei so gro}, daB die Masse der Welle vernachlassigt werden 
darf. Es ist bekannt, da8 fiir kleine Auslenkungen r, des Schwerpunktes 
S und fiir eine im Vergleich zum Tragheitsradius kleine Exzentrizitat der’ 
stationare Drehzustand (g=w —konstant) und stationarer Betrieb (d.h. 
resultierendes antreibendes Moment M =0) nahezu zusammenfallen ([5], 
S. 777; [2]. S. 287). Diese Aussage wird durch den Drallsatz geliefert, und 
wir brauchen hier, wo wir die Bewegung des Schwerpunktes ohne Riick- 
sicht auf die Kreiselwirkung untersuchen wollen, nur noch den Schwer- 
punktsatz anzuschreiben. Die einzige Kraft, die auf die Welle wirkt, 
rihrt her von der Auslenkung rz der Welle aus ihrer Ruhelage O und ist 
der Auslenkung proportional. Wir fassen nun die GréBen in der Scheiben- app. 3, Krifte und Bezeich- 
ebene als komplex auf, so da®B wir zwischen der Kraft und der Aus- bunges/ ax! der Scherbe- 


lenkung die Beziehung 


K=—mar, 


anschreiben kénnen. Die Masse m haben wir willkiirlich aus dem Proportionalitatsfaktor heraus- 
gezogen. Der Schwerpunktsatz lautet dann 


r,t ar,=aeexp(iwt), (8.1) 
wobei wir beniitzt haben, daB wegen y = wt die Beziehung 
Tg =T, —eexp(i a t) 


gilt. Als Lésung der Differentialgleichung erhalten wir 
t, == Ay exp (i Jat) + A, exp (— iya aes —— exp(tmt). 


Es treten nun zwei Faille auf, die sich grundsatzlich unterscheiden: 
1. Wenn die GroéBe a reell ist, so daB® also die Auslenkung und die Kraft in dieselbe Richtung 
fallen — dies wird der normale Fall sein —, so stellen die beiden ersten Glieder der Lésung 
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harmonische Schwingungen dar, wahrend die partikulare Lésung bei der kritischen Drehzahl 
W,p—= Va ein r, ergibt, das linear mit der Zeit iiber alle Grenzen wachst. 

2. Ist dagegen a komplex, was z. B. vorkommen kann, wenn die Welle auf Torsion beansprucht 
ist, so bewirkt eines der beiden ersten Glieder exponentielles Anwachsen von r, fiir jede Dreh- 
zahl und zwar derart, daB die Geschwindigkeit jederzeit senkrecht auf der Winkelhalbierenden 
zwischen Auslenkung und Kraft steht. Es mu somit jede Drehzahl als kritisch angesprochen 
werden; dagegen bleibt jetzt die erzwungene Schwingung beschrankt. Ihre Amplitude kann, 
wie man sich leicht iiberzeugt, immer noch betrachtliche Werte erreichen; fiir eine fliegende 
Welle von 4m Lange ist z. B. Anax 2, 200 e (vgl. [2], 5S. 288). 

Wenn wir nun das Verhalten von Wellen, die auf Torsion beansprucht sind, bei verschiedener 
Lagerung untersuchen, so handelt es sich nach dem, was wir oben ausgefiihrt haben, lediglich 
darum, den Zusammenhang zwischen einer Belastung und der Auslenkung des WellendurchstoB- 
punktes D zu ermitteln. Damit wir bei Beriicksichtigung der Kreiselwirkung diese Fragen nicht 
nochmals aufnehmen miissen, seien sie gleich hier mit der nétigen Allgemeinheit untersucht. 

Wir belasten die Welle am DurchstoSpunkt D durch eine Kraft K und ein Moment M, deren 
Vektoren in der zur Welle senkrechten (x, y)-Ebene liegen. AuBerdem sei die Welle durch ein 
Moment W auf Torsion beansprucht. Fiir die Auslenkung erhalten wir dann zwei Differential- 
gleichungen vierter Ordnung, die bei bekannten Randbedingungen sofort zweimal integriert 
werden kénnen. Die Berechnung der Auslenkung kann wegen des Einflusses von W bei gewissen 
Randbedingungen recht miihsam sein, bietet aber nichts Interessantes, so daB wir uns damit 
begniigen, die Gleichungen und die Lésungen fiir zwei Falle als Beispiele anzufihren. 

Im Falle einer einseitig eingespannten Welle, die die Belastung am freien Ende tragt, lauten 


die Gleichungen ([2], 5. 281) 
On” Wy) = Kol 3)- MO 


ay” —Wx' —K,(l—2) +M,=0. 


Dieses System la8t sich mit der komplexen Schreibweise in die einzige Gleichung 


orl Sb Wor oi Kee a0 (8.2) 


zusammenfassen. Wir erhalten mit den Randbedingungen daraus fiir kleine Torsionsmomente W. | 


die Auslenkung des Endpunktes r (1) und die Projektion des Tangenten-Einheitsvektors auf die 
(x, y)-Ebene, die wir mit v bezeichnen wollen, beide in Abhangigkeit von K und M. Fiir die 
Kraft und das Moment ergibt sich dann 


K=2[12r—(6 iw ol), | 
iM= fe [(6 +i) r—(4—"7) wl, | (8.3) 


sofern wir nur die linearen Glieder in w mitnehmen. Die Gré8e w steht dabei als Abkiirzung 
fiir W/x und stellt ungeféhr den Drehwinkel dar, um den die Welle durch W je Langeneinheit 
verdreht wird. Fiir M — 0 kommen wir auf das schon bekannte Ergebnis 

__ 3a 3 
K=F(1 —yiwl)r 


3 
zuriick ([2], S. 287). 
Als zweites Beispiel betrachten wir die beidseitig gelenkig gelagerte Welle, die in ihrer Mittel- 
ebene durch K und M belastet wird und bei der nur die eine Halfte unter Torsion steht. Die 


Differentialgleichungen fiir die Auslenkung, komplex geschrieben, lauten dann fiir die beiden 
Stabhalften verschieden, namlich 


ar’ —iWr +4 (Ks iM)—0 


ar Ces iM)=0. 


Fiir den Zusammenhang zwischen Belastung, Auslenkung und Neigung der Tangente in der 


und 
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Stabmitte erhalten wir schlieBlich 


48 « ae el lla. 
Ko = ale givl)r+sgiwlol ; | 
8.4) 
8al9. Oe 
iMag itor —3 (1 —g5 vel) 2] | 


oder fir M=0 


48 
K = (Ltgy ier. 


9. Beispiele. Wir wollen nun die Ergebnisse zusammenstellen, die man fiir verschiedene 
Lagerungsarten der Welle erhalt. 

Nur in zwei der in Abb. 5 dargestellten Falle tritt in der bekannten Weise eine einzige kritische 
Drehzahl auf, namlich in den Fallen a und d, und zwar ist hier auch fiir groBe Werte von W, wie 
die Rechnung zeigt, die Auslenkung parallel zur 
angelegten Kraft. Obschon auch hier die Welle ~*~ — a 
in eine Raumkurve verbogen wird, ist dies Er- % 
gebnis aus Symmetriegriinden zu erwarten. Die 

ae | Ce 
Y 


a 


I 
' 


I 
. 


an 


i 


— = 
it 
DAD Vv Vv 
eee 
W 8 
Abb. 4. Die Belastung der gelenkig gelagerten Welle. Abb. 5. Verschiedene Lagerungen der Welle. 


Erniedrigung der kritischen Drehzahbl durch das Torsionsmoment ist von R. Grammel ([5], 5.800) 
angegeben worden. In allen iibrigen Fallen wird jede Drehzahl kritisch, indem die Auslenkung 
nicht in die Kraftrichtung fallt. Besondere Erwahnung verdient vielleicht noch der Fall e, bei 
dem das exponentielle Anwachsen nicht mit wl erfolgt., wie in allen andern Fallen, sondern nur 
mit (w1)%, wobei die Amplitude nach einer Schwingung auch bei einer 10 m langen Welle héchstens 
um ihren 10-8-ten Teil ansteigen kann. Auch der Beiwert von (w/)® weicht tibrigens etwas von 
demjenigen ab, den R. Grammel fiir diesen Fall angibt. 


10. Die Beriicksichtigung einer Dampfungskraft. Der Einflu8 einer Dampfungskraft labt 
sich leicht tibersehen, wenn wir sie proportional der Geschwindigkeit des Scheibenschwerpunktes 
annehmen und setzen 

R= my f,*, 
wobei y eine Dampfungszah] bedeutet. 

Die Bewegungsdifferentialgleichung (8.1) lautet dann 


7,-yr,tar, = aeexp(iwt). 


Wir fragen uns nun, welche Dampfung gerade geniigt, um das exponentielle Anwachsen der 
Schwingungen zu verhindern. Dazu miissen wir die Lésung der homogenen Differentialgleichung 


oder ihre charakteristische Gleichung 
etyota=—0 
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untersuchen. Man zeigt leicht, daB der Realteil von o negativ wird, daB also die Losung abklingt, 
wenn die Ungleichung 
(Im a)? <y? Rea 
erfiillt ist. Bei einer festen Dampfung kommt es also darauf an, ob der Wellenquerschnitt, der 
durch eine Einheitskraft belastet wird — die Scheibenmasse m werde gleich Kins gesetzt —, das 
stabile Gebiet im Innern der Parabel (Abb. 6) verlaBt oder nicht. Dies hangt aber einzig von der 
Lagerung und wl ab, dagegen nicht von der Winkelgeschwindigkeit w. 
Bei verhaltnismaBig kleinem Imaginarteil von a, was bei unsern Fallen immer zutrifft, ist 
die Ungleichung gleichbedeutend mit 
2|Imya|<y. 


B. Die kritischen Drehzahlen bei Beriicksichtigung der Kreiselwirkung. 

11. Koordinatentransformation und kinematische Beziehungen. Durch die Kreiselwirkung 
wird unsere Aufgabe wesentlich verwickelter. Auch hier diirfen wir nicht einfach nach den 
Drehzustinden fragen, in denen die Welle bei verschwindender Exzentrizitat der Scheibe einer 
stationdren Auslenkung fahig ist ([5], S. 804), sondern miissen wirklich das dynamische Problem 
in Angriff nehmen. Abb. 7 stellt die Scheibe in der allgemeinen Lage dar, mit den aufern Kraften 
an der Scheibe, die zu einer Dyname St, | : 
Mzim DurchstoBpunkt D zusammengefaBt 
sind. Sie rithren her vom antreibenden 
Moment und der Elastizitat der ausgelenk- 
ten Welle und enthalten méglicherweise 
ebenfalls ein Antrieb- oder Bremsmoment, 
das ‘auf den Scheibenrand wirkt. Dieselbe 


M=Mz+Rxe 


Sy 


\\ 


Abb. 6. Der EinfluB der Abb. 7. Die Bezeichnungen an der Schzibe und die Dyname 
Dampfung. in D oder S. 


Belastung ist auch noch nach dem Schwerpunkt S reduziert eingezeichnet, da wir sie fiir den 
Drallsatz in diesem Punkt benétigen. (Die Belastung ist also eine Dyname in D oder in S.) 
Aus welchen Anteilen wir die Belastung im einzelnen zusammengesetzt zu denken haben, wird 
sich im Laufe der Rechnungen ergeben. 

Als Bezugsysteme wahlen wir das raumfeste System (x, y, z) und das mit der Scheibe fest 
verbundene System (&, 7, ¢), dessen -Achse mit der Exzentrizitat e des Scheibenschwerpunktes S 
und dessen ¢-Achse mit dem Tangenten—Einheitsvektor ) an die elastische Linie im DurchstoB- 
punkt D zusammenfallt. Um Vektoren von einem System ins andere transformieren zu kénnen, 
schreiben wir uns die Transformationsformeln ([2], S. 271) \ 


a; == a, cosy + a, sin p — a,(v, cosy + v, sing) , | 

a, = — a, sing + a, cos p + 4,(v, sin p — vy cos G) , Cita) 

a, = + 4, v, + ay vy + a, | 
auf. Der Winkel gist dabei der Winkel, der zwischen der x-Achse und der &-Achse liegt, wenn wir 
die Scheibe aufrichten, d. h. wenn wir die €-Achse in die z-Achse tiberfithren. Ferner wird voraus- 
gesetzt, daB die Komponenten v, und vy des Tangentenvektors so klein sind, daB ihre Quadrate 
gegen Eins vernachlassigt werden kénnen. Mit andern GréfBen derselben GriSenordnung werden 
wir ebenso verfahren, ohne jedesmal besonders darauf hinzuweisen. Die Transformation ist 
orthogonal und normiert, so daB die inverse Transformation sofort angeschrieben werden kénnte. 


Es empfiehlt sich auch hier, in der (x, y)-Ebene und der (€, 7)-Ebene die komplexe Schreib- 
weise einzufiihren, wodurch wir allerdings zur Unterscheidung der beiden Ebenen gezwungen 
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sind, etwas ungebrauchliche Bezeichnungen zu verwenden. Es bedeutet 
fiir einen beliebigen Vektor a: a,y =, +ia,, 


fiir den Ortsvektor tr: r=x-+1y, 


und da auch hier keine Verwechslung méglich ist, v= v,, =v, +iv,. 
Es lassen sich dann die ersten beiden Gleichungen von (11.1) zusammenfassen 2u 
ce, (a,, ta) exp (—i). (11.2) 
Die Bewegung der Scheibe kann in jedem Augenblick beschrieben werden durch die Kinemate 
im Schwerpunkt S, bestehend aus t, und der Winkelgeschwindigkeit u (a, , Wy, @,;). Fir die 
wirkliche zeitliche Ableitung eines Vektors gilt allgemein 
dy da 
vA aE +uxa, 
wo d’q/dt die Ableitung in bezug auf das Fahrzeug (die Scheibe) und 11 die Winkelgeschwindigkeit 
des Fahrzeugs bedeutet. Wenden wir diese Beziehung auf an, fiir das die scheinbare Ableitung 
verschwindet, so erhalten wir 
Vy, = Wy —W, Vy, Vy = Wz, Ve — Wy 
oder 
0 = —1(@, y — @, 0). (13) 
Nun beriicksichtigen wir weiter, daB die Rotation der Scheibe vorwiegend um die ¢-Achse erfolgt, 
daB also auch gilt 
|Oxy| KO. 
und damit nach (11.1) w, = w,; diese beiden Winkelgeschwindigkeitskomponenten, die @ gleich 
sind, sollen fortan einfach mit w bezeichnet werden. Fiir den Gleichlauf, bei dem sich die ganze 
Figur wie starr um die urspriingliche Achse dreht, ist sogar w,— w,—=0. Auch fiir die als Gegen- 
lauf bezeichnete Bewegung ist unsere Aussage jedenfalls erfiillt, indem dort 


0,=20—o=—0 wd ,,=2 00 
gilt ([5], S. 806). 
Aus (11.2) und (11.3) folgen zwei Beziehungen, die wir im nachsten Abschnitt brauchen 
werden. Fiir die Winkelgeschwindigkeit ergeben die beiden Gleichungen 


Me, = (Wry — 0) exp(—i—) =tvexp(—ip), (11.4) 


und daraus folgt weiter 
De, = (Lv + 0G) exp (— ig). (11.5) 


12. Die Kriterien fiir kritische Drehzahlen mit Kreiselwirkung. Um die Bewegung der Scheibe 
zu untersuchen, stehen uns der Impulssatz und der Drallsatz zur Verfiigung. Der Impulssatz 
lautet mit m als Scheibenmasse 

mt,=8, 
davon lassen sich die x- und y-Komponente zu 
mr, = Kyy (ize1) 


zusammenfassen. 
Den Drallsatz schreiben wir im mitrotierenden Bezugsystem, d. h. als Eulersche Gleichungen 


an. Wir setzen eine so flache Scheibe voraus, da zwischen den Tragheitsmomenten die Be- 


ziehungen gelten 
: Cer 2B 


wodurch sich die Eulerschen Gleichungen zu 


Ba,=— Bow,+M,;, 
Bo, =+Bow,+M,, 


ergeben. St stellt das Moment der Dyname im Scheibenschwerpunkt dar. Die beiden ersten 
Gleichungen lassen sich zusammenfassen als 
B Wey = 1Bwm;,+ Ms,. (12.2) 
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Man ersieht iibrigens aus unserem Vorgehen, daf wir von einer schiefen Aufkeilung der Scheibe 
auf die Welle abgesehen haben. 

Wir setzen nun fir die weitere Untersuchung stationaren Betrieb voraus; die ¢-Komponente 
des Momentes ty (ein méglicherweise an der Scheibe auftretendes Bremsmoment ist darin schon 
beriicksichtigt) soll also verschwinden. M, besteht demnach nur aus der €-Komponente des 
Momentes von & beziiglich dem Schwerpunkt. Ist die GréBe e?/k? sehr klein gegen Eins, so 
schwankt — kleine Auslenkungen des Schwerpunktes vorausgesetzt — die Winkelgeschwindigkeit 
so wenig um den Mittelwert jw, daB wir sie als unabhangig von der Zeit betrachten diirfen. Wir 
nehmen also an, daf stationdrer Betrieb mit stationarem Drehzustand (@ = 0) zusammenfalle 
([5], S. 777) und g = wt wird. 

Da das resultierende antreibende Moment verschwinden soll, bleiben in der Dyname in D 
nur noch die Beitrage, die von der Elastizitat der Welle infolge der Auslenkung geliefert werden. 
Sowohl die Kraft wie auch das Moment fallen in die (€, 7)-Ebene, und es gilt 


Ke KO exp (101) 
und 
M;, = M,,exp(—tot), (12.4) 


wobei, wegen K, = 0, M;, in D und S gleich gro8 ist. 

Den Zusammenhang zwischen der Auslenkung, der Neigung der Tangente und der Belastung 
der Welle haben wir im Abschnitt 8 hergeleitet. Dieser Zusammenhang ist fiir kleine Aus- 
lenkungen linear. Wir schreiben ihn hier in allgemeiner Form, um verschiedene Lagerungen 
gleichzeitig zu erfassen: 


= Kay MUO T a Le Y, | 
iM,,=Bdrj+Bbv. | 


Die GréBen a, b, c, d sind, abgesehen von m und B, die wir aus Dimensionsgriinden willkirlich 
hinzufiigen, die dualen GréBen zu den Maxwellschen EinfluBzahlen. 
Setzen wir nun (12.5) in (12.1) ein, so erhalten wir 


r,t+argtev=0. 


(12.5) 


Es ist aber 

rq =, —eexp(t wt), (12.6) 
und damit wird die obere Gleichung zu 

7, + ar, +cv=aeexp(iwt). (12.7) 


Wir wollen nun ebenfalls (12.2) ins feste Koordinatensystem tiberfithren und beniitzen dazu 
die Gleichungen (11.4), (11.5) und (12.3): 


(tBu +2 Buaw) exp(—iwt) = M,,exp(—iot). 
Daraus ergibt sich mit (12.4), (12.5) und (12.6) 
v—2iwv+bv+dr,—deexp(iwt). (12.8) 


Diese Gleichung stellt zusammen mit (12.7) das Differentialgleichungssystem dar, das unser 
Problem beherrscht. Man kénnte es leicht auf eine einzige Differentialgleichung 4. Ordnung 
in r, zuriickfithren. Eine partikulare Lésung liefert uns der Ansatz 


TT, CXP(COl)e und, v = 9 exp (i wt) . 
Er fiihrt auf das inhomogene lineare Gleichungssystem 
r,(—w*+a)+v,c¢c=ae, 
rp5d+v,(w?+ 6) =de 
fiir ry und vg, und daraus erhalten wir 
2 =e 
ie oe eS a cd’ * (12.9) 


Die kritischen Drehzahlen treten genau dann auf, wenn die Determinante verschwindet, naim- 


lich bei 


T9 


(of — a) (wi tb) Fed=0. (12.10) 


XX. Band 1952. Trésch: Stabilitatsprobleme bei tordierten Staben und Wellen. 213 


Die Lésung des homogenen Differentialgleichungssystems erhalten wir ebenfalls durch einen 
Exponentialansatz 


T,=Toexp(tAw@ot) und v=v,exp(id ab) (12.11) 


Darin verwenden wir die dimensionslose Gri®e A, vervielfacht mit der kritischen Drehzahl Wop 
die ohne Kreiselwirkung und beim Fehlen des Torsionsmomentes auftritt. Damit wir nicht- 
triviale Lésungen erhalten, mu die Determinante des Gleichungssystems verschwinden : 


(A? we — a) (A? w2 — 21 0,0 —b)—ced=0. (12.12) 


Weil «, eine positive Zah] darstellt, wachst die Schwingung unbegrenzt, wenn der Imaginarteil 
von A kleiner als Null wird. 

Die Gleichungen (12.10) und (12.12) erméglichen uns, fiir einen vorgegebenen Stab die 
kritischen Drehzahlen festzustellen, sofern wir den Zusammenhang zwischen Belastung und Aus- 
Ienkung kennen. Wir untersuchen nun die beiden Faille, den einseitig eingespannten Stab und 
den hbeidseitig gelenkig gelagerten Stab, fiir die wir in Abschnitt 8 diesen Zusammenhang be- 
rechnet haben. 


13. Kritische Drehzahlen fiir die fliegende Welle. Fiir die fliegende Welle liefert ein Vergleich 
von (12.5) mit (8.3) fiir die Einflu8gréBen 


a=—4a?, c= — 0-3 (6 —iwl), 

: (13.1) 
_ : er ie Ek 
d gO Piwls b= abr (4 ). 


Die darin neu verwendeten Abkiirzungen bedeuten 9 =k/l, wo der Tragheitshalbmesser 
Be 3 6 
k Ale ist, Una! = i ([2], S. 288). 


Wir bemerken hier nur nebenbei, daB durch das Torsionsmoment die Reziprozitatsbeziehung 
mc = Bd, die sonst gilt, zerstért wird ([5], S. 95). 

Bei der Frage nach den kritischen Drehzahlen stellt man fest, da8 sich aus der partikularen 
Lésung der Differentialgleichung keine kritische Drehzahl ergibt. Um dies nachzuweisen, spaltet 
man die Gleichung (12.10), die jetzt 


3 owt + af (4—12 9 Fat Lot (— 442i 01) =0 (13.2) 


lautet, in Real- und Imaginarteil auf. Wir werden am SchluB dieses Abschnittes darauf zuriick- 
kommen. 

Die Untersuchung der allgemeinen Lisung wird uns zeigen, daf fiir jede Drehzahl eine Lésung 
mit negativem Imaginarteil vorhanden ist, die Beriicksichtigung der Kreiselwirkung also grund- 
satzlich nichts an der Tatsache andert, daB durch das Torsionsmoment jede Drehzahl kritisch 
wird. 

Die charakteristische Gleichung (12.12) der homogenen Differentialgleichung lautet fiir die 
fliegende Welle mit der Abkiirzung u = w/w, 


1 iwl 4, 
(#4) (2 2Au za\4 =) z=0. 


Um uns einen Uberblick iiber die Verhaltnisse zu verschaffen, entwickeln wir die Wurzeln A 
nach der kleinen Grobe iwl: 

A=A, tiwlsa,ters. 
A, geniigt dabei der Gleichung 


(22 — 4) (i —2Au— 5) =9 (13.3) 
oder 
(28 — 4) A, (2 w) = 55 (4 —1)- 
und fiir 
pene ey 


2 a(iw 1) 
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erhalt man 


1 42 ; 

rae 12 Ga — a | 9(A,— u)|. (13.4) 
Gerade fiir dieses A, interessieren wir uns; denn nach Gleichung (12.11) ist die Schwingungs- 
amplitude nach der Zeiteinheit auf r,(1) = r,(0) exp(—A/, wl wo) gestiegen. 

Um A, in Abhangigkeit der Winkelgeschwindigkeit und des Traégheitsmomentes darstellen 
zu kénnen, miiBte 2, aus der Gleichung (13.3) bestimmt und in (13.4) eingesetzt werden. Wir 
ermitteln den Zusammenhang jedoch einfacher graphisch. Den beiden Gleichungen (13.3) und 
(13.4) entnehmen wir ohne Rechnung, daf 2, und ebenfalls /, lediglich das Vorzeichen wechseln, 
wenn man u durch — wersetzt. Wir wahlen u > 0, was keine Kinschrankung der Allgemeinheit 
bedeutet. Um A, in Abhangigkeit von u mit 9? als Parameter der Kurvenschar auftragen zu 
kénnen, berechnet man mit (13.3) aus A, die GréBe 
dat tee ARES) (13.5) 
2 30,3 —4) 
und daraus das zu wu gehorige J, mit (13.4). 

Abb. 8 zeigt, da® bei gré®eren Tragheitsradien der Scheibe die Aufschaukelung langsamer 
erfolgt, daB jedoch die Kreiselwirkung die Welle nie vollstandig stabilisieren kann. AuBerdem 
wachst die Amplitude umso lang- 
samer, je gréBer die Drehzahl ist. 
Abb. 8 umfaBt, nach links und 
nach rechts aufgetragen, gleich- 
zeitig die Falle w< 0 (Antrieb 
an der Welle, Bremsung an der 
Scheibe) und w > 0 (Antrieb an 
der Scheibe, Abnahme des Mo- 
mentes an der Welle). Da Abb.8 
nicht symmetrisch ist, besteht 
zwischen diesen Fallen ein Unter- 
schied. Es sei gleich hier erwahnt, 
daB bei der gelenkig gelagerten 
Welle, wie man aus (16.3) ent- 
nimmt, diese beiden Méglich- 
keiten sich nicht unterscheiden. 
Auf eine Eigentiimlichkeit, die 
sich fiir kleine Tragheitsradien er- 


u= 


tu 
450 +40 ~=-+302=«+20~—=«tC«iS CHOC SD gibt, werden wir bei unserem 
Abb. 8. Das Anwachsen der Amplitude bei der fliegenden Welle. zweiten Beispiel (Abschnitt 16) 
zuruckkommen. 


In unseren Ergebnissen sind die Sonderfalle 9 = 0 und w — 0 ebenfalls enthalten. Wenn wir 
die Kreiselwirkung vernachlassigen, d. h. 9 = 0 setzen, so miissen wir zu den Ergebnissen von 
Hi. Ziegler ([2], S. 288) gelangen, und dies ist tatsichlich der Fall. Andererseits miissen die Er- 
gebnisse in diejenigen von R. Grammel iibergehen, wenn wir die Wirkung des Torsionsmomentes 
vernachlassigen, d. h. w = 0 setzen. Die Gleichung (13.2) gibt nach w/w, aufgelést die bekannten 
kritischen Drehzahlen des Gleichlaufs ([5], S. 806). Dagegen erhalten wir die kritischen Dreh- 
zahlen des Gegenlaufs nicht, obschon unsere Voraussetzungen sie keineswegs ausschlieBen 
(vgl. S. 271). Der Gegenlauf laBt sich jedoch nur deuten als Resonanzerscheinung, erzeugt durch 
auBere Schwingungen. 

Bei der Beurteilung, ob eine Drehzahl als kritisch’ angesprochen werden muB oder nicht, 
darf ein Punkt nicht auBer acht gelassen werden. Als kritisch haben wir solche Drehzahlen be- 
zeichnet, bei denen sich die Scheibe mit der Zeit beliebig weit von der Ruhelage entfernt. In 
Wirklichkeit geniigt es jedoch schon, daf die Scheibe um einen gewissen endlichen Betrag aus 
der Ruhelage ausweicht, um den Rotor zu gefahrden. Ohne Beriicksichtigung des auBeren 
Momentes braucht man auf diese Unterscheidung nicht zu achten, weil sich diese gefahrlichen 
Drehzahlen unmittelbar an die eigentlich kritischen Drehzahlen anschlieBen. In unserem Fall 
kénnen aber sehr wohl gefahrliche Drehzahlen auftreten, obschon eigentliche kritische, wie wir 
gezeigt haben, bei der partikularen Lésung ausgeschlossen sind. Die GréBe der Auslenkung be- 
rechnet sich aus den Gleichungen (12.9) und (13.1). Den gréBten Wert nimmt |ro| offenbar 
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an fiir 
3 9? wt + we(4 — 12 0%) w —404 = 0. 

Setzen wir noch w?= wq, wobei 1 < q (0?) <2 ([5], S. 807), so wird 
8 (1+ 3074) 


Holmes = Ta @—q ~ 


und [Fo|max 14Bt sich in die Grenzen einschlieBen 


33 
8 ay +9 
eh Sg ace =< ag 


Die untere Grenze gilt fiir den Fall, daB die Kreiselwirkung vernachlassigt werden darf. Die 
Gefahrdung der Welle wird also hier durch die Kreiselwirkung erhéht. Eine Abschatzung von 


Tolmax zeigt, da& man hei fliegenden Wellen normaler Lange auch im giinstigsten Fall mit 
Boksax a 200 e rechnen muf. 


14. Der Einflu8 der Dampfung. Das unerwartete Ergebnis, daB bei der fliegenden Welle 
jede Drehzahl kritisch ist, wird bewirkt durch das Torsionsmoment W oder also durch die gegen 
Kins kleine Gré8e wl. Es ist deshalb zu erwarten, da8 auch schon eine geriuge Daimpfung die 
Verhaltnisse grundlegend verandert. Wir wahlen wiederum als einfachsten Ansatz die der 
Geschwindigkeit proportionale Dampfung. Unter dieser Annahme schreibt sich dann (12.1) 


mr, -myt, = K,,. 


Die Rechnungen bleiben nun Schritt fiir Schritt dieselben wie in Abschnitt 12. Die partikulare 
Lésung liefert entsprechend (12.10) eine kritische Drehzahl, wenn die Bedingung 


(Wj; — ty @, — a) (wap +b) +ced=0 (14.1) 
erfiillt ist. 


Die charakteristische Gleichung des homogenen Differentialgleichungssystems lautet an 
Stelle von (12.12) 


(A? we —idwyy — a) (22 w2 —2Aw,w —b)—cd=0, (14.2) 


und die Schwingung wachst wiederum exponentiell an, wenn der Imaginarteil von ( kleiner als 


Null ist. 


15. Die Dampfung bei der fliegenden Welle. Um die Wirkung der Dampfung bei der fliegenden 
Welle zu iibersehen, untersuchen wir die Gleichung (14.2). Wenn wir namlich eine Dampfung 
,einschalten‘‘ und die Dampfungszahl immer gréBer werden lassen, so bleibt vorerst immer 
noch jede Drehzahl kritisch, das Aufschaukeln erfolgt jedoch langsamer, wahrend gleichzeitig 
die partikulare Lésung immer gefahrlicher wird. Bei einer bestimmten Dampfungszahl y,, die 
wir als Grenzdampfungszahl bezeichnen wollen, besteht die Lésung der homogenen Gleichung 
aus ungedampften harmonischen Schwingungen; dagegen wird jetzt die gefahrliche Drehzahl 
zu einer eigentlich kritischen. Fir eine Dampfung, die gréBer als die Grenzdampfung ist, stellt 
die Lésung der homogenen Gleichung abklingende Schwingungen dar, und die eigentlich kritische 
Drehzahl ist wiederum nur noch mehr oder weniger gefahrlich. 

Die Grenzdampfung yj, bestimmen wir durch die Forderung, daf in (14.2) A reell sein soll. 
Fiir die fliegende Welle fiihrt die Aufspaltung in Real- und Imaginarteil auf die beiden Glei- 
chungen 

4 4 en, 
(8 — 4) (# —2du 55) — a Flu! 
und 


l 4A 
(e— 4), ate (# —2Au —sa|= 0. 


Die Dampfungszahl] wird, wie vermutet, proportional zu w 1. Dadurch fallt in der obern Gleichung 
die rechte Seite weg, und ein Vergleich mit (13.3) zeigt, da A zu A, wird. Fir die Grenzdiémpfung 
ergibt sich 
ta Ai 4) 
Yo 944, 


Wywl. 


20 
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Fiir jede Drehzahl erhalten wir allerdings vier verschiedene Werte von A,, also auch vier még- 
liche yo; die wirkliche Grenzdampfung wird durch das gr6Ste unter ihnen beschrieben und kann 
mit Hilfe von (13.5) berechnet werden. f 

Man kann zeigen, da unser — wie wir schon betont haben — sehr unvollkommene Dampfungs- 
ansatz auch bei der Annahme von groBen Dampfungszahlen nicht imstande ist, die Schwingung 
vollstandig zu dampfen. Es 1a48t sich bei jeder festen Dampfungszahl eine Drehzahl angeben 
derart, da® alle héheren Drehzahlen kritisch sind. Diese Grenzdrehzahl hangt noch vom 
Tragheitsradius ab. 


16. Kritische Drehzahlen bei der gelenkig gelagerten Welle. Als zweites Beispiel fiir die Unter- 
suchung von kritischen Drehzahlen wahlen wir die beidseitig gelenkig gelagerte Welle. Dabei 
soll das Moment W nicht am andern Wellenende, sondern an der Scheibe abgenommen werden, 
so daB nur die eine Halfte der Welle unter Torsion steht (Abb. 4). 

Die EinfluBzahlen haben wir in Abschnitt 6 berechnet. Die kritische Drehzahl ohne Torsion 
und Kreiselwirkung betragt bei dieser Lagerung 


9 484 
Se ery 
Damit ergibt sich aus (8.4) und (12.5) 
a= a4 (1 +3q 502), = oa iwl, 
one 1 are 
2 at 2 
Ae 0 cane b= wheal gpl). 


Das weitere Vorgehen ist gleich wie bei der fliegenden Welle. Die Gleichung (12.10) lautet 
3 Sue 1 ae 
oe a(t +3g¢ 01) c og a(} givl)/=0, (16.1) 


und daraus ergibt sich keine eigentlich kritische Drehzahl. 
Als charakteristische Gleichung fiir das homogene Gleichungssystem erhalten wir aus (12.12) 


a—( +3qit02)| |#—24u—7 (1 —gpé wl} =O. (16.2) 


Hine wesentliche Anderung gegeniiber der fliegenden Welle stellt hier die Tatsache dar, daB — 


das Glied cd wegen seiner Gré®enordnung wegfallt und deshalb die charakteristische Gleichung 
aufspaltet. Die Bewegung des Schwerpunktes ist mit der Bewegung des Tangentenvektors nur 
durch das Torsionsmoment, in unserer Naherung also sehr schwach, gekoppelt. Fiir die Unter- 
suchung der Stabilitit der Welle greift man am besten auf die Gleichungen (12.7) und (12.8) 


zuriick. Von Interesse sind wiederum besonders die Imaginarteile A, der Wurzeln 


3) 6 1 twl 
— (1 iwi), j= 8 /u2 | a 
=\IT 6G Jai Zot tog (aut 
da sie das Aufschaukeln der Schwingung je Zeiteinheit bestimmen: 
3 1 
<= i= =e 
128 9 V4 u?9?+1 


Wir bemerken, dafs der eine Wert von der Kreiselwirkung und der Drehzahl unabhangig ist. 


Ay = aIyarns 


(16.3) 


Der andere Wert zeigt uns eine Kigentiimlichkeit, die schon bei der fliegenden Welle (Abb. 8) | 


auftritt. Fir kleine Tragheitsradien (u? 9? <1) ist 


Il 


A= + 199 o(1-+2 uo)’ 


woraus wir schliefen miissen, daB das Ausweichen aus der Ruhelage fiir sehr kleine 0 auBer- 
ordentlich rasch erfolgt. Diese Lésungen miinden fiir 9 > 0 demnach nicht in diejenigen ein, 


die wir ohne Kreiselwirkung aus (8.1) erhalten, und man muB bei der Belastung durch ein | 


Moment W die Kreiselwirkung grundsatzlich immer mitberiicksichtigen. Ist kein Torsions- 
moment vorhanden, so spielen diese Schwingungen sehr hoher Frequenz allerdings keine Rolle. 

Wir streifen noch kurz die Ergebnisse, die wir ohne Kreiselwirkung erhalten. Die Glei- 
chung (16.2) zeigt, da natiirlich auch in diesem Falle jede Drehzahl kritisch ist. Vernach- 
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lassigen wir die Torsion der Welle, so ergibt (16.1) die bekannten kritischen Drehzahlen; die 
Kreiselwirkung spielt keine Rolle, wie aus Symmetriegriinden zu erwarten ist. 

Auch bei dieser Lagerung kann eine gefahrliche Drehzahl auftreten an der Stelle, wo sich 
die urspriingliche kritische Drehzahl befindet. Aus (12.9) folgt fiir den gréBten Ausschlag 


[roles = 577° 
O|max 3 |wl| ° 


wobei wir beriicksichtigt haben, daB die gréBte Amplitude fiir w = w, auftritt. 

Bei dieser Lagerung sind also grundsatzlich die Verhiltnisse genau dieselben wie bei der 
fliegenden Welle. Die kritische Drehzahl, die ohne Torsionsmoment vorhanden ist, wird zwar 
immer noch gefahrlich, aber nicht mehr eigentlich kritisch. Dagegen wird jede Drehzahl kritisch ; 
wenn die Kreiselwirkung iiberhaupt einen Einflu8 hat, so wirkt sie in stabilisierendem Sinne. 


17. Die Dimpfung bei der gelenkig gelagerten Welle. Wenn wir nun auch fiir die beidseitig 
gelenkig gelagerte Welle den Einflu8 unseres Dampfungsansatzes untersuchen, so stellen wir 
fest, daB wir auch hier nur die charakteristische Gleichung (14.2) des homogenen Differential- 
gleichungssystems zu betrachten brauchen. Das Aufspalten der Gleichung durch das Wegfallen 
des Gliedes cd bewirkt hier, da®B zwei Wurzeln von der Dampfung vollig unberihrt bleiben, 
namlich die Wurzeln von 

(A? w2 —2AwW,w— b) =0, 


welche deshalb notwendigerweise zu Instabilitat fihren. Eine Dampfungskraft, wie wir sie 
eingefiihrt haben, geniigt in diesem Falle niemals zur Stabilisierung, wie gro8 wir auch die 
Dampfungszahl wahlen. Es ist méglich, diesem Ubelstande abzuhelfen, wenn man ein Moment 
einfiihrt, das das Wippen der Scheibe abddampft. Wir begniigen uns hier mit der Feststellung, 
da8B ein solcher Ansatz wirklich zum Ziele fiihrt. 

Die Grenzdampfung y,, die wenigstens die abdampfbare Schwingung gerade stabilisiert, 
berechnet sich leicht zu 


3 
Yo= 3q/4 Oe 


Die nicht abdampfbare Schwingung wird iibrigens nach (16.3) mit zunehmender Drehzahl immer 
ungefahrlicher. 


(Eingegangen am 12. Februar 1952.) 


Anschrift des Verfassers: Dipl. Phys. A. Trésch, Institute for Mathematics and Mechanics, New York 
University, 45 Fourth Avenue, New York 3, N. Y. 
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Zur Theorie der rechteckigen Fundamentplatten und Pilzdecken *. 
Von W. Miiller. 


1, Einleitung Fiir die Durchbiegung von Platten und Pilzdecken hat bekanntlich V. Lewe 
nach der Fourierschen Methode Doppelreihen aufgestellt, die aber wegen der schlechten Konver- 
genz fiir die praktische Rechnung wenig geeignet sind. Um die Voraussetzungen méglichst all- 
gemein zu halten und méglichst viele Sonderfalle zu erfassen, wollen wir im Folgenden von der 
Annahme einer elastisch nachgiebigen Unterlage der Rechteckplatte ausgehen und die Um- 
wandlung der dafiir giiltigen Leweschen Doppelreihen in einfache gut konvergierenden Reihen 
methodisch und in direkter Weise mit Hilfe einiger Grundformeln ableiten und uns dabei in 
dieser vorlaufigen Mitteilung zunachst auf die Falle der gleichmaBigen und der Punktbelastung 
beschrinken, zur Vorbereitung auf die in Aussicht genommene Behandlung der Falle mit recht- 
eckiger Angriffs- bzw. Stiitzflache, die mathematisch einige Schwierigkeiten bietet}. 


2. Methode der Doppelreihen. Wenn wir als Reaktionswirkung der Unterlage eine Kraft an- 
nehmen, welche der Durchbiegung w proportional ist, so lautet die grundlegende Differential- 
gleichung 


Ll Eh 


Cap ee : phinvgtont> derelbbid 
AAw + L w= mit ie pea (1) 


wobei E den Elastizitatsmodul, h die Plattendicke und y die Poissonsche Konstante bezeichnet. 
Die Gré8e L hat die Dimension einer reziproken Lange; der Faktor ¢ = L* N (kg/cm*) wird im 
Allgemeinen auch als Bettungsziffer? bezeichnet. Im Falle der schwimmenden Platte ist c mit dem 
spezifischen Gewicht des fliissigen Mediums gleichbedeutend. 

Wenn man von der klassischen Lisung von Navier absieht, so hat namentlich V. Lewe® die 
Methode der Fourierschen Doppelrethen ausgebildet. Sie besteht im wesentlichen darin, daB man 
die Druckfunktion, die als stetig oder stiickweise stetig vorausgesetzt wird, in eine trigonome- 


MUX . NIT 


mI x n It 
sin. = bzw. cos yx 


trische Doppelreihe, etwa nach sin cos = mit den Perioden 2a und 


2b entwickelt, wahrend fiir w ein entsprechender Reihenansatz mit unbestimmten Koeffizienten 
gemacht wird. Nach Bildung von AM4w und Einsetzung in (1) ergibt sich dann durch Gleich- 
setzung zweier Reihen die Méglichkeit, die unbestimmten Koeffizienten zu ermitteln. Dabei 
treten vor allem zwei Falle hervor, der Fall der gleichmaBigen Belastung der Gesamtplatte und 
der Fall der gleichmaBigen Belastung langs rechteckiger Druck- oder Stiitzflachen, der in der 
Grenze in den Fall der punktférmigen Belastung oder Stiitzung durch Einzelkrafte ttbergeht. 

a) Fir die gleichmafige Belastung einer drehungsfrei aufgelegten rechteckigen Platte mit den 
Seiten a und b und den Randbedingungen Aw = 0 fiir x = 0,a und y=0,b haben wir mit 
E = x«/a,n = y/a, A = b/a die Druckfunktion 


16 p er. . 
P= > > pn i mas sin ne (m,n =1,3,5...), (2) 


wenn der Ursprung der Koordinaten in die linke untere Ecke des Rechtecks gelegt wird und 
die Achsen den Rechtecksseiten parallel sind. Mit dem Ansatz 


wa D) D, Cn sin marten naz. (3) 


* Die vorliegende Arbeit ist einer Reihe von Untersuchungen entnommen, die in Verbind i 
Institut fiir Beton u. Stahlbeton der Technischen Hochschule Karlsruhe (Prof. Dr. K. Komen 
gefiihrt wurden. 

* Vel. fiir das Folgende: A. Nadai, Die elastischen Platten, Berlin 1925; K. Girk Fila 

2 , 5 5 dS hentrag- 
werke S. 176—178f, Wien 1948; S. Timoshenko, Theory of 248, Newly onlite 
arene an imoshenko, eory of Plates and Shells S. 239, 248, New York a. 

2 F. Schleicher, Beton u. Hisen, 26. (1927) S. 433f.; K. Hayashi, Zur Theorie des Tri i 

Unterlage, Berlin 1921. Weitere Literatur wird an spaterer Stelle anceusiones aes ecg 


* V. Lewe, Bauingenieur, 1 (1920) S. 37; 1 (1920) S. 631; 3 (1922) S. 314 u. 344; 4 (1923) S. 453. 
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findet man dann, wie bekannt, durch Koeffizientenvergleich 
be 16 py b4 
Cn “2m n[(n2x2 + m2 22 gu)? + (L b)*] (4) 
und daraus die Durchbiegung 


, 
sin maz €sinna > 7 


_- 16 po b* 
Maree m n [rd (n? + m®72)2 + (Lb)] sone) ©) 


b) Bei der lokalen Belastung auf einem Rechteck, dessen Mittelpunkt die Koordinaten 

Xq = & @, Yo= No@ hat, und dessen Seiten 2a a und 2a sind (Abb. 1), hat man die Belastungs- 
funktion in der Form 

P = Pof (*) &(y) ( 


anzusetzen, wo fiir fund g die Reihen 


f(a) = >" a, sin maz é& , g(y) = >! b, sin nm 7 


gelten sollen. Dabei erhalten wir fiir die Koeffizienten a,, und b,, nach Fourier 


(Piece gar | f(#) sin mé, b= | 80) sin nov ‘ (8) 
f 
Bei der Integration ist z B. f(x) fiir 
Xy— AA<K< xp aa 


gleich Eins und auBerhalb dieser Grenzen 


ZY 
Di 


Abb. 1. Aufliegende und lokal belastete 
Platte. Abb. 2. Fundamentplatte. 


Null zu setzen. Entsprechendes gilt fiir g(y). Man braucht also nur von %)— @a bis x) + aa, 
bzw. von y,— a8 bis y, +-af/ zu integrieren und erhalt dann durch eine einfache Rechnung 
Ashe , 4, 
a, = —— sin mré&) sin MO, b, = sin na sinnne. (9) 
Damit wird 


pix ¥) = 2B > D sen im ma Ep sin mara sin nar Pain nae sin mack sin no | 


; (10) 
(he eee 
und schlieBlich ergibt sich durch Koeffzientenvergleich 
een sinma &,;sinmaa sinnz ‘o.sin na F sin ma Esin na tL 
__ 19 Po ® : (11 
yeas 2. (11) 


Soye G + m? 22)? 4 (5) 


Wenn a und f gegen Null gehen und 4 p, a? af =P gesetzt wird, so haben wir bei punktartiger 
Belastung 


; 7 
sin ma &) sin nz sin ma sin na 


(nijpre= 1 2530.4) 28 © (42) 


oe 


2 2 A2)2 1 Lb\* 
(may + (F 
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c) Fir die Fundamentplatte (Abb. 2) und Pilzdecke (Abb. 3) wollen wir hier eine unendlich 
ausgedehnte und von einem regelmabigen Sdulengitter mit rechteckigen Grundquerschnitten, 
belastete bzw. gestiitzte Platte voraussetzen. Legen wir den Ursprung in einen der Saulenmittel- 
punkte und sind 2a und 26 die Abstande der benachbarten Gitterpunkte in der x- und y-Richtung | 
oder die Seiten des Periodenrechtecks, so treten in den Reihen die Elementarfunktionen cos maz & 


und cos nz + auf, sodaB die Ableitungen dw/dé dw/dy in den Randpunkten = 0, 2; y= 0, 2A | 


verschwinden. Fiir die Belastungsfunktion 
setzen wir 


P(~y¥) = af(*) 8) (13) | 


Yj : 
A) mit 


f(x) FH + > dy cos mir, | 
- nt 1 (14h 


1 
B(¥) => bo + >>) bn cos ac 


Ca 


Fiir die Auswertung der Fourierschen Inte” 
grale gilt dieselbe Bemerkung wie frither: 
Man erhalt dann durch Rechnung 


igaune 


DS oS 
dn =—— sn MTC, 
mI 


2B (15) | 
bb = sen | 
Abb, 3. Pilzdecke. bee < sin ie . 


Damit wird die Druckfunktion 


2qB wl 
Pay)=o +> 


yi at ae 
sin m7 0 Cos mag 4+ “4% S$’ —sin na ® cos na = | 
i (16) 
| 
J 


m 
4 ire s 
+i > > —sin mz o sin es con MEE cos Hye 
m mn 7 7 
™m n 


Der Konstanten entspricht dabei eine konstante Einsenkung. Im Falle der Pilzdecke (Abb. 3) 
ist q negativ zu nehmen und L — 0 zu setzen. Fiigt man noch die gleichmaBige Belastung py, 
hinzu und beriicksichtigt, daB aus Gleichgewichtsgriinden py'4ab = q@af 3 pp = a B/A ist. 
so heben sich die Konstanten auf. Macht man nun fiir die Durchbiegung w einen entsprechenden 
Ansatz mit zwei einfachen und einer Doppelreihe mit unbestimmten Koeffizienten, so findet man 
durch Vergleich beider Seiten der Differentialgleichung folgenden Ausdruck, in dem wir den Fal] 
der Fundamentplatte und Pilzdecke zusammenfassen: 


; 1 
| ; - sha lecmne adil 
2q bt) B sinmazacosmazé& A A 


Ww = Wy ae matin + CO] | Ds nalenh SLO 


(17) 


: . a 7 
sin marasinna © cos mI > COS nail 


dit Ye mnt mE E LO | (m,n = 1,2,3...). 


Dabei ist das obere positive Vorzeichen fiir den Fall der Fundamentplatte und das untere negative | 
Vorzeichen fiir den Fall der freien und gestiitzten Pilzdecke mit L = 0 zu nehmen. 
Far den Sonderfall der Punkthelastung (1 = 6 > 0; q+ 4a2a 8 = P) ergibt sich 


; | cos na cos mx Ecosna 
Pb cos ma & 


7 
E 2aN| Sim) (Lbe iG ee (De ee me Pace Fm? 2)? + (L by] (18) | 

(int, 1012 eee) 
Die auftretenden Doppelreihen haben eine schlechte Konvergenz, besonders nach einmaliger und | 
zweimaliger Differentiation und sind deshalb namentlich zur Bestimmung der Momente nicht 
sehr geeignet. Es wird daher unsere Aufgabe sein, die Reihen nach einer Koordinate zu sum- 
mieren oder die Doppelreihe in eine einfache Reihe mit mals Laufzahl zu verwandeln. 


W = Wy 4 


| 
| 
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3. Methode der einfachen Reihen im Fall der aufliegenden Platte unter Vollbelastung. Zur 


oar der Doppelreihe (5) in eine einfache Reihe gehen wir aus von der Kneserschen 
orme 


SS aes =: | 


Ue ae 2 Acoma 
n? + @ n? + @2 20? ' 20 Gnon 


set 8 te n (a — 0) 1 a Cojo O 


(Gee O a7). (19) 
Integriert man beide Seiten nach @ im Intervall von 0 bis, so erhalten wir 
sin n(x — @) C) Sin w O 
2 n(n? + w*) 2 w ( ~ Sin va ; Y) 
Aus (20) geht die weitere Formel hervor 


<1 = sin n(x — 0) 1 Lo 1 il 
va n(n? + @?) (n?+@’2) m2 —@”? ox Fares iz 0) é +o? n+ =) | 


me ‘ 


2 w? w/? 
die wir ohne weiteres auf die Reihe (5) anwenden kénnen, wenn wir noch folgende Umformung 
vornehmen. Setzt man namlich 


(21) 


I =: ie y 
Sin wa Gin w’ x 


| il fe GnwoO w Gina’ 1 


w* — qm’? 


7 (n® -- m? A?)? + (L b)* =a! (n? + 2 m? 22) (n? + hk’? m? 2) , (22) 
so erhalt man 
k =i) asia ba! =|i—i(ey. eh? = 1 +(22) . (23) 


Mit k=y,+1 0m, ko = ¥m—t Om wird 
ym = \F(Va+ E42). = V P(r +E). (24) 


L ; ae iat ‘ : Lei Lea 
vt oO, = le =e ee) 5 Va Om = Va One wes (25) 


mI 


Ersetzen wir jetzt in den Formeln (20) und (21) ~— © durch ay/A und » durch kmA, baw. ow’ 
durch k’ m/, so ergibt sich 


sw sin nw : Ae 
n (n? + m2? k2 7?) (n2 + m? hk’? 7?) L a\4) 
a Teh (= ‘) | (26) 


4 


| 
T 


‘Cinkma(A—n) Sink’ ma (A —7) 
ea k? Gink mad k’? Gin k’ mad ] 


A(mA)*i 2 
Setzt man fiir k und n’ die konjugiert komplexen Werte (25) ein und entwickelt die Funktionen, 
so treten gemaf den Formeln 


oj c tid) = oj ccosd +i Ginesind, Gin(c + id) = Sinecos d +i Wo) csind 


hyperbolische und trigonometrische Funktionen auf. Nach einer etwas umstiandlichen, aber 
nicht schwierigen Rechnung ergibt sich schlieBlich bei Benutzung der Beziehungen (22) bis (25), 
wenn man noch mit 16 p, 6¢/Nx® m multipliziert und iiber m summiert, fir die Durchbiegung 
folgende Formel?: 
SECy. | (22) Pm + Qn | 
8 sin m It n 
pan Le > 1 ‘ | 


~ REE pa (aay I Zon le 


(gfe = Leo a) (27) 


1 Vgl. A. Kneser, Die Integralgleichungen und ihre Anwendungen in der mathematischen Physik, S. 148, 


hweig 1911. : 
Be auch S. Jguchi, Eine Lésung fiir die Berechnung der biegsamen Platten S. 46, Berlin 1933. 
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wo die Funktionen P,,, Q,, und Z,, folgende Bedeutung haben: 
P,, = Sit ym m7 sin bm, mI (2A —y) — Sit Ym ma (2A —y) sin 6, man , 
Qn = Cf im MN COS Om ma (2A — 1) — CO} Ym mz (2A — 7) cos bm MIN , (28) 
Zm = C0) 2 Ym mxA—cos26,, mu. 


Die Formel kann noch vereinfacht werden, wenn man unter der Voraussetzung der analytischen 
Fortsetzbarkeit den Mittelwert der positiven Durcbhiegungswerte fiir die Punkte 7,, den an der 
Seite 71 = 0 gespiegelten Punkt 7, — 1, und die beiden um b davon abstehenden Punkte 
N3=A+n, Ny =A— NH bildet (vgl. Abb. 4.), also 


w = [10 (n) + 10 (—thg) — 20 (79) + 20() 


Dabei ergibt sich fiir die beiden ersten Glieder in der Klammer 


L ee M 3 ta]. 1 
i|) oe ak ie ee 2 


Das Glied mit 7 fallt also heraus. Ferner ist, wie eine einfache Rechnung ergibt, 


Pr (1j1) + Pn(—= 12) — Pm (18) + Pn (Ha) = — 4 (CO) 7m mx A — cos 6, mI A) Cyn (7) i 


Qn (111) + Qm(— M42) — Om (773) + Om (is) = — 4 (CO} 7m mA — cos bn, mI 2) Hm (1) 5 (29) 
Zim = 2 (Col ym ma A — cos bm mI A) (Gof Ym MAA + cos 6, MY) . 
Daraus folgt dann 
4 p sin mat & (3) RC hae | 
as 0 : — 
w= HTD, ms0\4 Coj ym mad + cos mma UF eae (30) 
Sai (2) 
\L a/ | 
wo 
Gn (11) = SI Ym moe(A — 7) sin bn my + Sit Ym m7 sind, ma(A—7) » | (31) 
Hn (7) = CO} Ym m(A — 17) €08 dm mae + CO} Ym may cos dn ma(A—n). J 
‘Papp 4, Die Aufteilung der gefundenen Lésung. Zur Bestatigung der 
Rechnung kann man noch ein anderes Verfahren! anwenden, das 
zugleich Gelegenheit gibt, einige bemerkenswerte Sonderfalle zu be- 
trachten. Wir gehen zunachst aus von dem Fall des Plattenstreifens 
3¢-- - (6 = oo), dessen Lisung w, als Partikularlésung der inhomogenen 
Gleichung zu gelten hat. Fiir diesen Fall haben wir also die Differential- 
| ¥ gleichung 
f d4w Pp 
! oa + T4(w— Hin) =0, (32) 
Aaa Ph naan SMiellans = 
| ie sich bei konstantem py mit w —7 74 = u auf die gewohnliche linear- 
aan j= ==- homogene Gleichung | 
H d‘n | 
ee ee Gee ail ealag (33a) 
| 
! reduziert, deren charakteristische Gleichung die Wurzeln | 
| el ie 4 ; 
eo Agel eee 
hat. Wenn man die daraus sich ergebenden Partikularlésungen 
L Gol “ecose Gin Le since 
Abb. 4. Zur Umf V2 V2 V2 
ia ‘Falle depuVollbelaceane® | Coj oe in Be ; Gin Lx one Lx 
/2 2 /2 
benutzt und die Randbedingungen w = 0, w’’ = 0 fiir « =0, a beachtet, so ergibt sich als 


* E. Estanave, Contribution a de I’ étude l’equilibre élastique d’un plaque etc., Theses prés. a la Faculté 
des Sc. de Paris 1900; ferner M. Lévy, C. R. 129 (1899), S. 535. 
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Lésung 


aie 1 


Lx 1 
ears = oss ee x) + Cvj = (a — x) cos = 
he [2 ia 2 an 
E Gof we -L cos ae ; 


la 2 y2 2 


Um diese Lésung in eine sin-Reihe von der Form 


. x 
w= >) AQ sin mT — 


zu entwickeln, haben wir nach Fourier die Integrale 


r= [oat Te 8g (a — x) sin ma— dx, r= [on en) cose sin “2 dx 
auszurechnen. Die beiden Integrale sind gleichwertig mit den Summen 
i =7@itB)=z [cos = (a —x-+1%) sin ma — duty | cose (a— x—1x) sin ma — dx, 
0 
1 L 1 , L 
, , . ° x . 5 x 
(Sy 5) = 5 | cos 7 (x-+i(a—«))sin ma = dx Bee | Ware (x —i(a—a)) sin ma ~ dx ; 
0 


Durch zweimalige Anwendung der partiellen Integration ergibt sich, wie hier nicht naher aus- 
gefiihrt werden kann, 


S, +32 = = gant ie > 
es 
mz\? aL a (35) 
ee 2 ea ai (—1)™ cos a 
eek Te 
Damit wird 
y ; (™) Airs a a z ("2) a a 
+ (814+3,4+ + 9%) = one (60)*2 + cos $=) - Ba (601 22 + e052) 
(m =1,3,5...). (36) 


Da sich auch die Eins reihenmaSig darstellen la8t in der Form 


4 ere x jis 
=— >'— sin ma — Game leseor (37) 


0 ergibt sich fiir die Durchbiegung des Plattenstreifens folgender Ausdruck: 


i sin ma = 
+ —Po = 38 
Wy 4 . (38) 
ee Gel 


Um die vollstandige Lésung fiir die begrenzte rechteckige Platte zu erhalten, haben wir noch 
die Lésung der homogenen Gleichung 
AAw + L4w=0 


hinzuzufiigen. Setzen wir 
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so geniigt Y(y) der Differentialgleichung 


yay) 2 al Y""4 Ls ! E “)] noe 
a a 


Mit der frither eingefithrten Bezeichnung (22) bis (25) haben wir dann partikulare Lésungen von 
der Art 


Co} ym m7 COS Om MIN ; Co) yn may sin bn MIN , 


und in Verbindung mit w, kénnen wir daher die Gesamtlésung n’ = w,-+ w, in der folgenden 
Gestalt ansetzen: 


OE Seager 2 


wo fiir die Funktion F'(7) bei Beriicksichtigung der Randbedingungen w=0, Aw = 0 fiir y = 0,2 
wieder derselbe Ausdruck gewonnen wird, den wir friiher unter (31) cash haben. 


Eine gewisse Bestatigung der Rechnung erhalt man, wenn man den Grenzfall a =oo durch- 
fiihrt. Verschiebt man den Anfangspunkt des Koordinatensystem in den Punkt (a/2, 0), so ist 


m—1 


sin mae(a + :) = (-- 1) 2 cos maiaé 


an Stelle von sin ma & zu setzen. Da hierbei 


m—l1 
(—1) 2. aL! Sa 
> cos mag = (15D sea 

ferner 

Timmy mac = Timm by, mad ae lim re =) 
und 

, or L L = 2 co 

lim Hn(1) aa: (b— y) cos eo OT 85 = (bh —y) 


wird, so erhalt man in Ubereinstimmung mit der Formel ne fiir den Plattenstreifen parallel 
der x-Richtung 


; iy La ¢ 36 
fe Sol ze — 9) e057 + Co} TE gO 


pa NED? Go bL + bL 
=+te 
airy OS == 2 fl 


Das ist gleichzeitig die Barchhicsdite fiir den ange Balken von konstanter Biegungssteifigkeit 
N=EJ, dessen beide Enden auf festen Stiitzen frei aufliegen, wahrend seine Unterflache 
einen elastischen Untergrund beriihrt. 


lim 2 = 


(40) 


5. Die auf elastischem Untergrund gelagerte Platte mit Punktlast. Zur Umformung des Aus- 
drucks (12) gehen wir aus von der Formel 


5 Tome 7) 1 = Tt M4 
sin na sin UA See 1) Joos 5 — +7) toss - (A—n —%)|- (41) 
Dann wird entsprechend der friiher besprochenen Zerlegung des Nenners 


sin nw sin na) 


1 
Dies Dare im? 22) (n® + k’2m? 12) — (made — hk) * 


il - nwt nat : 1 il 
x > o (11) eos oie (A—7)+) —cos | (A— nN n| a eta he 72 [n® + ke cl ; 


(42) 
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Nach zweimaliger Anwendung der Kneserschen Formel (20) ergibt sich nun 


a v 
4i(ma A)? ees 
mm 
Coj k’ ma(A— ny + n) — Gof k’ ma (A — —n) Coj k ma(A— ny + n) —Coj k ma (A —1 —n) 43 
2k’ maxi Cink’ mad 2Zkmani Cinkmad | if) 
= il | Sin k’ mx (A —m) Cink’mnay  CGinkma(A— mm) Gink may 
Dine Ay za; | 2k’ mand Sink’ mad 2kmand Cinkmad | 
mI 


Bei Einfithrung von k = y,, + i6,, k’ = Ym —id, und Entwickelung der hyperbolischen Funk- 
tionen erhalt man schlieBlich 


4 PR 4P a @ sin maz &)sin mx E[ym Xm(7) + dm Ym(n)] 
Ne Net er oe ha AVF (44) 
ental ie | | (Coj 2 ym m x A — cos 2 dm mad) 
m7 mo 


wo die auftretenden Funktionen folgende Bedeutung haben: 
Xn (4) = Sit ym mu (A — No) cos Om MX(A — 7H) + F (ny) 
— Gof yn mx (A — Ho) sin bn mu(A — Ho) G(x) » 
Yn (y) = Sin y,, ma (A — No) cos Om mx (A — Yo)-G(n) 
+ Co} ym mx (A — no) sin 6, mz (A — Yo) F (7); 
mit (45) 


F(n) = + (Gin Ym m(A+ 7) sin dp, mz (A —n) —Sin ym m(A —7) sin bn ma (A-+ n)] » 
G(7) = 5 [Go Ym MaI(A+- 7) COS Om ma(A — 7) — CO} Ym ma(A — 7) cos J, ma(A+ 7) - 
Die Formel (44) gilt zunachst fiir 7< 7), wahrend die entsprechende Reihe fiir 74> ny nach dem 


Reziprozitatsgesetz aus (44) durch Vertauschung von 7 und 7) hervorgeht. Die Ableitung nach 7 
erleidet dabei keine Unstetigkeit. Man hat, wie mau leicht feststellt, 


Pa sin ma &,sin ma E[ Gin 2 ym ma(A—7p) sin 2 dm mI Hq + Sin 2 ymmaA Hy sin 2 dmma(A— 1) ] 
(ma)8 ei (Coj 2 ym mad — ocs 26mm) 
mar 


(46) 
Wenn die Last in der Mitte der Platte angreift, erhalten wir mit €) = 1/2; 7%) = 1/2 


Pes sin 5 marsin mat &(ym Xm+ Om Ym) 
«= Na (Ee 
m m? os 


\ JU 


—— (ie ale: Stes (A) 
i+ ey (Gof ym ma A + cos dm ma A) 

m It 
wo die Funktionen X,, und Y,, iibergehen in 


Xe CU, maa COS Om mas: F (y) — @of Yn mas sin Sq me G(n) , | 


AS: A 
Yn = Cit, mas, COS Om mas: G (4) + Go} Yn mot ->- sin On mo + F(n) . | 
Die Durchbiegung im Lastpunkt wird dann 


sint > mz (dm Gin ym mA — ymsin Om ma A) 


wp = eee ez fate (49) 


™ JU. 


2 
Fir die gewéhnliche freie Platte (L = 0) erhalt man mit y, = 1, 0m =9, 2 O= ie) die fiir 


Ingenieur-Archiy 
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4 XN giltige Formel 


Pa CDosinmazé,sinmaé x; , 
w= po a ee {Gin ma(A — no) (1+ mad Cotg mad) Sin may 


m= 


— mx Co} may] — ma(A — no) Col mx(A — 9) Sit may} , (50) 


wahrend fiir 7>7) die Formel 


ee 3k = sin mz &)sin ma & A Bf diss ¥ >: 
w= Ta, aie ee {Gin ma(A — y) (1+ mad Cotg mad) Sin ma ro 


m—L 


— mx no Co} ma yo] — ma (A — 1») Co} ma(A—7) Sit max yo} (51) 
zu benutzen ist!. Wir haben in Abb. 5 die Hoéhenschichtlinien der Biegeflache fiir den Fall 
E5=1/2 mo=1/2, A = 1,5 gezeichnet. Bemerkenswert ist, daf die Stelle der gréBten Durchbiegung 
nicht mit dem Lastpunkt zusammenfiallt. 

Ist der Lastpunkt in der Plattenmitte (£5 = 1/2 , 79 = //2), so wird 


° mim. a 
Pa on a ma & 
v= —- - ———. -——-.— 
IN 2s n® Sof mae 
[Sin max (1 + Ue: rq rales mz Wo) may}. (50a) 
Im Lastpunkt erhalt man dann 
: 1 ase 
sin? —_ ma(Gin mad —mad) 
2 
eee Pa Dee 2 
ave m3 Go}? mah 
D 
~ mah (52) 
= P@ 8 oo l mah 
273 .N m3 ( ved) 
Wiebe acen beh Ris Hd, 


Bei der quadrischen Platte (A= 1) hat man z. B. bei Beschrankung 
Abb. 5. Héhenschichtlinien der Biege- 


fliiche einer frei aufliegenden Platte mit auf das erste Glied der gut konvergierenden Reihe 
Punktbelastung 


1 1 
4A=1,53 Ep= 1 y= 1. 
2 2) 


twp = 0,0108 =~ , (53) 


6. Die Durchbiegung einer durchlaufenden Platte mit Punktbelastung auf elastisch nachgiebiger 


Unterlage. Die unter (17) aufgestellte Doppelreihe kann unmittelbar durch die Knesersche Formel 
(19) umgestaltet werden, wenn man von der Darstellung 


cos nal 
A — ss: 
Dem + k2 m2 72) (n? + km? 22) | 


= battad 5 >, cos na -- : 4 
(mad)? (k? — k’?) — d\n? + kh’? (mad)? nn? x? + hk? (mad)? 


ausgeht. Man findet dann mit den frither eingefiihrten Bezeichnungen nach Trennung des Reellen 
und Imaginaren 


(54) 


Yi 


cos nm 
A 


ri vm Rm(n) + Om Tm(7) 1 
2mari(Lb) 1 i Size 2¢ma aye ( 
™m It 


2a (n? + m?2?)? + (L b)4 L *) | : (s3) 
m I 2 
wo folgende Abkiirzungen eingefiihrt wurden: 

Rin(]) = Of ym may sin b, mI (2A —n) + Cod) ym ma(24— n) sin b,, man , 

T,, (4) = Sit ym mH COS Om mI (24 —n) + Sin yp» ma(2A — 7) cos dm man , (56) 

Zn =o) 2y_,maA—cos26, mn. 


1 Vegi. wach K. Grein, Pilzdecken, S. 34, 2. Auflage, Berlin 1948. 


| 
| 
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Man sieht daher, daB sich die erste Lewesche m-Reihe heraushebt und die Formel geht daher im 
Falle der Fundamentplatte iiber in 


Hier ist w)eine Konstante, w, die partikulare, von x unabhangige Lésung, die einem der x-Rich- 
tung parallelen Plattenstreifen von der Breite 2b mit einer Einsparung in y = 0 und 26 und einer 
gleichmafig verteilten Last entspricht. Die GréBe w, gibt als allgemeine Lésung der homogenen 
Gleichung die Durchbiegung einer in x = 0 u. 2a; y = 0 u, 2 b eingespannten Einlastplatte. 


7. Der Fall der Pilzdecke mit Punktstiitzen. Zur Bestatigung der Rechnung wollen wir noch 

den Grenzfall L = 0, a = 6 = 0 betrachten. Man erhalt dann z. B. 
- YmRm(y) +-6mTm(n) oj ma(A—n) [1 + mad Cotg mx A] —ma(A—n) Sin ma(A — 7) 

lim 2 ; (60) 
as Z Sin ma A 
Wenn man den Anfangspunkt in die Mitte des Stiitzenrechtecks legt, also an Stelle von € und 7 
die Summe 1-+ § und 7+ A setzt, so geht wy iiber in 
__. PF a* (—1)™ cosma& 
— 4Nz3 m Cin.ma, 


. 


[Co] man (1+ mad oj mada) — may Sin may]. (61) 


We 
Fir w, haben wir dagegen mit L > 0 


PB 
274 Na n4 


lim w, = 
L-0 
Nun ist 


wo (—1)" cos no 


2\2 
4 = konst. — 7 (1 Fr] : 


“at nf 48 
a 


Wir erhalten also, von einer Konstanten abgesehen, 
ia ee 22 _ 1 aia ae 3 
w= Pe — yy =x ptt —(2)] (63) 


Die Durchfithrung des Grenziiberganges an dem unter (58) stehenden, bereits summierten Aus- 
druck fiir w, ist etwas umstandlich und soll nur kurz skizziert werden. Denn man mu8 hierbei im 
Ziahler die drei ersten Glieder in den Reihenentwicklungen fiir die trigonometrischen und hyperbo- 


U/] 
en de cos ma & m Rm (7) + 6 Te Me ens 
w= wW a UA 2 (7) m m() | {P Ah 
; INTE 2 Peegy Zm '2Niab ao) sy (57) 
ees eS 
Fiir die Summe im letzten Glied erhalt man aus der allgemeinen Formel, wenn man k? = i, 
OD . 1 : / 1 : ; 
k’? — —i, also k = 5 ee Dae es ie (1 —i) und fiir die als konstant zu bewertende GréBe 
mz Aden Wert L b einsetzt, die Reihe 
cos nat 
ee 
aay. 
1+ (F5) 
Ly. L eect, Ly L Li L iB 
SoS = SS (AD SO) Set PA S11 Si J — Sin —= = a eee 
3 j a se ( y) +o} 3 (2b — y) sin a + Sin 5 cos a (2b —y)+ 5 (2 b—y) cos 2 ‘ | 
Coj /2 Lb —cos 2 Lb Li 
(58) 
SLs SO | 
mae iH? Z(b) Lb’ 
Es ergibt sich daher im Ganzen folgender Ausdruck fiir die Durchbiegung : 
P K(y) P cosmmE = ymRm(y) + dmTm(n) 
= ! ee 
w= wot awibe ZO) TINT Oo Zm pe 
memati) 
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lischen Funktionen, d.h. die Glieder bis zur fiinften Potenz des Argumentes beriicksichtigen, 
wihrend im Nenner die ersten beiden Glieder ausreichend sind. Setzen wir also 

fia 2 zh Sinz = | 28 eee 

eee ee oe Bieta pee 


cos 2 = sin z =| 


so erhalten wir 


: 9) L& ((b—y)(b+y!t , (b+y)(b—y)* (b+)? (6—y)* 
lim Ky) +208) = p18) =a ae ROP Bch = | (64) 
(+ yt (by)? Oye Oe Oe =i konst i 
24 L240 240 oH 48 


77 | 
| 
| 

ee 
AS = 
Se awe t= eae 

Boe Jn55 Pik ee on, ~195 719 
¢ ¢ ¢ . SLM = 
i.—o.4 ee --—_- 

Ss . x 20 2198156 7 

2S Re ee Wi ee ya ae cael 


Abb. 6. Héhenschichtlinien fiir die Biegefliche einer durchlaufenden Pilzdecke 
mit Punktstiitzen (A=1,5). 


Damit wird in der Tat 


: Ee, For ae 
Setzen wir w, und w, in unsere Reihe (54) ein, so erhalten wir fiir die gewéhnliche, gleichmabig 


belastete und punktartig gestiitzte Pilzdecke, bezogen auf den Mittelpunkt des Stiitzenrechtecks, 
folgenden Ausdruck fiir die Durchbiegung: 


co 


bt 2\2 3h <2, {care £ 
w= atopy} rr) + pa ) sie ODivis (1+ mx Cotg mx A) Coj man | 


m> Gin mad 


(65) 
—man Sin man). 
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Fiir die Mitte (& = 7 = 0) haben wir 


ay po | pab ome 
enc | Nx 2D, wiGi ey a Tal Cote mm A), (66) 
oder bei Beschrankung auf das erste Glied der Reihe mit CotgaA x1 


pb | 24 1l+aA 
aN Day Sines le (66a) 


We Wy 


Die Konstante wy kann etwa dadurch bestimmt werden, dai man die Stiitzen als starr, also die 


Durchbiegung hier gleich: Null setzt. 


Wir haben die Héhenschichtlinien der Biegungsflache nach (65) fiir die Falle A = 1,5 und A=1in 
den Abb. 6 und 7 zur Darstelllung gebracht. 


Abb. 7. Héhenschichtlinien fiir die Biegefliche einer durchlaufenden Pilzdecke 
mit Punktstiitzen (4 = 1). 


8. Sonderfall einer auf elastischer Unterlage ruhenden Fundamentalplatte mit einer Einzelkraft- 


reihe. Wenn man in den Ausdriicken (59) die Seite bins Unendliche wachsen l]aBt, so erhalt man 


Bes ee Rm(q) + 6m Tm(q) __ eo mm an (Ym Sin On m7 -+ Om COS Om MIN) , | 


b— co Zm 
_by (ep 
erry) Vee net Ley, Ly 
(ir 2" ag Sit COS ||. 
ZG) a) 


Damit ergibt sich, wenn man gleichzeitig das Vorzeichen der Krafte umkehrt, die Durchbiegung 
einer durchlaufenden Platte mit einer einfachen Kraftreihe in Ubereinstimmung mit Nadai! und 
Westergaard? in folgender Form: 


Sia eee E 
eee P e (sin 2 + cos 72) 


4/2 L3aN /2 V : 
| P cos mm &e~?m™™" (ym sin Om mx + Om COS Om mm N) (68) 
ae ae T a : 

sea | 1 +(=) 


1 4, Nadai, a. a. O. S. 188. 
2 H. M. Westergaard, Ingenioren 32 (1923), S. 517f. 
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Die gréBte Durchbiegung im Punkte x = y = 0 wird, von einer Konstanten abgesehen, 


es 
wr = ae +4 reps) ES ‘ (69) 
™ TU. 


Interessant ist ferner der Grenzfall der unendlich ausgedehnten Platte, die eine Kinzelkraft tragt. 
Der Biegungspfeil fiir diesen Fall ergibt sich aus (69),wenn man a ins Unendliche wachsen aft. 
Der erste von w, herrithrende Teil wird dabei Null. Der andere Teil geht mit as ite = == duin 
ein Integral iiber, Man erhalt also 


fo) 
. 


/VTA u4 — u? 
We = Bane / eee dn. (70) 
2/2n NL, 4+ u 
0 
Wenn man nach Westergaard die Substitution 
a ate oii 
2vV1+ v2 


anwendet, so erhalt man in Ubereinstimmung mit einer von H. Hertz! bei der Untersuchung der 
in einem Punkt belasteten schwimmenden Platte abgeleiteten Formel 


eee (11) 


“tf TeENIen la | 1+ 88ND 
0 
Der gréBte Bodendruck unterhalb der Platte wird 
p=NI'‘u,, (72) 
also in dem oben betrachteten Fall, wenn man noch die Nachgiebigkeitsziffer ¢ = L* N der 
Bettung einfiihrt, 


_ PE P)2cd—*) 73 
Nea gtr eof 8077, es 


9. SchluBbemerkung. Die weiteren Folgerungen aus den gegebenen Ansatzen, namentlich die 
Verfolgung des Verlaufs der Momente mit Hilfe der Thetafunktionen, insbesondere fiir den bemer- 
kenswerten Fall der Pilzdecke, sowie eine Methode zur Ableitung der Durchbiegung aus der 
Momentenflache durch eine reine Integration soll spater gegeben werden. SchlieBlich hofft der 
Verfasser auf die mathematisch nicht ganz einfache Behandlung der Fundamentplatten und 
Pilzdecken mit endlich-rechteckigen Angriffs- bzw. Stiitzflachen ausfihrlich zuriickzukommen. 


(Eingegangen am 25, Februar 1952.) 
Anschrift des Verfassers: Prof. Dr. W. Miiller, Miinchen 9, Taubenstr. 4 I. 


1 Vgl. H. Hertz, Wiedemanns Annalen, N, F. 22 (1984), S. 449; fe 
an tee, ( )s ; ferner H. Hertz, Ges. Werke, I. Bd. 


me) A sre : ae ‘lh has ee i Reine.) st ‘ 4 
bcatant Herat pret *) ie = eck bene he Siege 
eed aie Meme, Seca maar ay at Re 
rper. ‘Vou Professor Dr.-Ing. Friedrich Téike, __ 
: So meee Sears ag een eta tg TM he eS tReet eG 
Mit 339 Abbildungen im Text. VIII, 388 Seiten, . | ~ 
ete eee ak --- Ganzleinen DM 45,— = 
| as Sy Test Ke . ia ae / 5 es, 23 NG a 
: Vorwort zu dieser 2, Auflage den Wunsch duBert, daB 3 
5 der Mathematik und der Technik zu schlagen‘, so wird _ ; 
x Technik bekannt ist, was es heiBt, dic Kolle eines Vermitt- = | 
cht nur gerne bestitigen, dafS dieses Werk einen wertvollen = | 
’ auch seine Bewunderung dariiber zum Ausdruck bringen, = =9- |__ 
\ufgabe hier gelést wurd ores a nae es ies Viet ack ee a: = 
Zeitschrift des Vereins deutscher Ingenieure'* 
4 Li, * ott et ee ae, ay oe ta Ae Nie bas , 
mks e. aah we igs 


at: sged iickt. Wahrend der Mathematiker schon lange _ 
den meisten Ing; en noch fremd. Dieses. Buch will nun den Ingenieu- — 
durch eine lebendige und verstandliche Darstellung. Die Probleme _ 
Leser wird in den Stand gesetzt, damit zu operieren. Die vielen, 
1an mit Matrizen praktisch rechnet Das Buch wird dazu bei- 
| moderner Ingenieurmathematik wird. Die Darstel- 

uch das Buch in erster 


en Mathematikern sehr zu empfeblen. 
2 De a tett gc ctans A et 53 


Ry nH, > 


ie ae ee 
sor an der Universit 
PN, fo be eee ong Sse eeN ‘ 


it langem erwartete Buch ist aus Vorlesungen entstanden, die der 


en 
n Anschauung aus, ohne 
PSAs 


er Ebenen.ansc 


7 


LAG/BERLIN. GUTTINGEN - HEIDELBERG 


tal pap 


ER 
ees 


+ 


a 


Einfihrung in die Baustatik. Von Professor Dipl.-Ing. Dr. techn. E. Metan, 
Wien. Mit 242 Textabbildungen. X, 328 Seiten. 1950. DM 28,50; gebunden DM 31,50 


Aus den Besprechungen:... Das vorliegende Werk des bekannten Verfassers befaBt sich mit der Theorie 
der ebenen Tragwerke und setzt dabei die elementaren Begriffe der Statik als bekannt voraus, Auf Grund der aus der 
Kinematik der Tragwerkscheiben hergeleiteten Exkenntnisse wird die Verteilung der @uBeren und ‘inneren Krifte 


bestimmt. Behandelt werden weiter das Prinzip der virtuellen Verschiebungen mit seinen Anwendungen sowie die - 


Formanderungen und Biegelinien der statisch bestimmten und unbestimmten Systeme. Ausfihrlich entwickelt Melan 


das Verfahren der EinfluBlinien und erldutert die Ermittlung dieser Linien nach verschiedenen Verfahren. Die ge- 


wonnenen Erkenntnisse werden in ihrer Anwendung auf einige der am haufigsten angewandten Fachwerksysteme — 
den frei aufliegenden Triger, den Bogen- und Rahmentriger, den mit einem Stabbogen verstirkten Trager und den 
Durchlauftriiger — gezeigt. Das Buch erfreut durch seine klare und iibersichtliche Gliederung und durch die einprag- 
Same Behandlung des Stoffes. Es wird jedem, der sich mit dem behandelten Stoff erstmals vertraut machen will, gute 
Dienste leisten. Aber auch dem in der Praxis stehenden Ingenieur, hezouders dem Stakibaucr, wird es bei der Wieder- 


auffrischung mancher Kenntnisse yon Nutzen ‘sein. ; * ,,VDI-Zeitschrift® _ 
: = : e J ” 
Ebene und raumliche Rahmentragwerke. von p:.-1ng. V. Kupfersehmid. 
Mit 252 Abbildungen. VII, 196 Seiten. 1952. Ganzleinen DM 35,70 


Das Buch behandelt die Berechnung ebener und réumlicher Rahmentragwerke auf einheitlicher ES Ps Der Ver- 


fasser geht dabei von dem Nachweis aus, daB der durch die Belastung eines Rahmens verursachte Formanderungs- | 


zustand aus bekannten Formanderungszustinden abgeleitet werden kann, Als Ausgangssystem wird der ,,Rahmen 
mit verschieblichen, aber unverdrehbaren Knoten“ eingefiihrt, ein System, das sich fiir jeden Verschiebungszustand 
einfach; berechnen laBt, da in ihm alle Stabe in den Knoten fest eingespannt sind. Die praktische Anwendung der 
Methode wird durch eine Reihe von Zahlenbeispielen erleichtert, die so gewahlt sind, daB in Anlehnung an sie jede 
Rahmenform bei beliebiger Belastang berechnet werden kann, Die grundlegenden Gedanken der Behandlung der 
Rahmentragwerke in der vorliegenden Art folgen aus der ,,Deformationsmethode* von Ostenfeld (Verlag Julius 
Springer, Berlin 1926) und dem Momentenausgleichsverfahren von Cross. t 


Dynamik des Bogentrigers und Kreisringes. von Professor Dr.K. Feder- 
hofer, Graz. Mit 35 Textabbildungen und 26 Zahlentafeln. XII, 179 Seiten. 1950. DM 23,— 


Aus den Besprechungen:... Das vorliegende Werk des bekannten Autors stellt im wesentlichen eine Er- 
weiterung unserer Kenntnisse in der Schwingungsforschung dar. Wird doch im Gegensatz zu den bekannten Standard- 
werken tiber technische Schwingungsprobleme der Fall gekriimmter Stabe einer eingehenden Behandlung unterzogen. 
Der Verfasser beschrankt sich aber nicht nur auf die allgemeine Behandlung der Probleme, sondern er ist bemiht 
bis zur zahlenmiBigen Auswertung jedes Problems vorzudringen. So findet man auch die Niherungsverfahren in yoll- 
kommener Art durchgefiihrt und es ist dem Leser miéglich, mittels der zahlreichen Tabellen und Darstellungen jedes 
einschlagige Problem auch rasch tiberschlagsweise zu beurteilen. apa »Acta Physica Austriaca’ 


SummeneinfluBwerte fiir den einfachen Balken und den sym- 


metrischen Zweifeldtriiger fir StraBenbriicken, — von Ds. 10¢. 


Friedrich Schweda, Wien. Mit 46 Abbildungenim Text undin 10 Zahlentafeln. VI, 79 Seiten. 
1952. -. . $teif geheftet.DM 17,70 


Ausgehend vom einfachen Balken werden fiir den symmetrischen Zweifeldtriger fiir die Belastungsschemen von 
StraBenbriicken, die aus einem Einzelfahrzeug und einer vor und hinter diesem angeordneten Gleichlast bestehen, 
Tabellenwerte angegeben, mit deren Hilfe sich die gréSten Momente und Querkrafte solcher Trager in ahnlich rascher 
Weise wie fiir eine felderweise gleichmaBige Vollbelastung ermitteln lassen. Auf den EinfluB der Veriinderlichkeit des 
Traigheitsmomentes in des bei Stahlbetontragern tiblichen Form ist Ricksicht- ‘genommen, 
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